Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,898,122 C→T 5.0% A397A (GCC→GCT)  sdaA → L‑serine dehydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,898,1220CT 5.0% 258.5 / 2.0 79A397A (GCC→GCT) sdaAL‑serine dehydratase 1
Reads supporting (aligned to +/- strand):  ref base C (40/35);  new base T (2/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGT  >  NC_000913/1897975‑1898269
                                                                                                                                                   |                                                                                                                                                   
aaTGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCa                                                                                                                                                    >  7:222608/1‑149 (MQ=255)
aaTGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGCAATGCCa                                                                                                                                                    >  4:282557/1‑149 (MQ=255)
 aTGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCAt                                                                                                                                                   >  1:302846/1‑149 (MQ=255)
  tGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCTAtt                                                                                                                                                  >  1:270554/1‑149 (MQ=255)
  tGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGAATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCTAtt                                                                                                                                                  >  1:265760/1‑149 (MQ=255)
    gctgctGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGc                                                                                                                                                <  3:264420/149‑1 (MQ=255)
        ctgcGGGTCTTGCGGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGTAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCctct                                                                                                                                            <  3:300851/149‑1 (MQ=255)
         tgcGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtg                                                                                                                                           >  5:53532/1‑149 (MQ=255)
           cGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACTTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtgt                                                                                                                                          >  2:73959/1‑148 (MQ=255)
               tCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAgg                                                                                                                                     <  8:222608/149‑1 (MQ=255)
                 ttGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGc                                                                                                                                    >  1:306526/1‑148 (MQ=255)
                   gCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGAt                                                                                                                                 >  3:61265/1‑149 (MQ=255)
                     aGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtgtg                                                                                                                                         <  1:41120/139‑1 (MQ=255)
                     aGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCtgtg                                                                                                                                         >  2:41120/1‑139 (MQ=255)
                     aGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGAtt                                                                                                                                >  6:89957/1‑148 (MQ=255)
                          tgctgGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAAcgcc                                                                                                                          >  5:206059/1‑149 (MQ=255)
                            ctgGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAAcgccgc                                                                                                                        <  2:306526/149‑1 (MQ=255)
                                    tAGCCCGGAACAGGTTTGCGTGGGGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCTTTGAACGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATgg                                                                                                                <  2:269795/149‑1 (MQ=255)
                                    tAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATgg                                                                                                                <  6:206059/149‑1 (MQ=255)
                                        ccGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGctct                                                                                                            >  6:292315/1‑149 (MQ=255)
                                            aaCAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTgcgc                                                                                                        >  7:266257/1‑149 (MQ=255)
                                             aCAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTGTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGcgcc                                                                                                       <  3:300928/149‑1 (MQ=255)
                                                 gTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCAc                                                                                                   >  4:65485/1‑149 (MQ=255)
                                                 gTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCAc                                                                                                   >  1:42834/1‑149 (MQ=255)
                                                 gTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCAc                                                                                                   >  3:218498/1‑149 (MQ=255)
                                                    tGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCa                                                                                                 <  5:292315/148‑1 (MQ=255)
                                                     gCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGt                                                                                               >  8:94371/1‑149 (MQ=255)
                                                       gTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTg                                                                                              >  8:227147/1‑148 (MQ=255)
                                                        tGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCa                                                                                            >  1:301364/1‑149 (MQ=255)
                                                            ggcggAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCAcc                                                                                          >  1:241057/1‑147 (MQ=255)
                                                              cggAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACcgcgc                                                                                      <  4:52422/149‑1 (MQ=255)
                                                                        cATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTg                                                                             >  1:232796/1‑148 (MQ=255)
                                                                          tGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGAt                                                                          <  5:119530/149‑1 (MQ=255)
                                                                          tGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGAGCCGTGCATTGAGCGTAATGCTATTGCCTCAGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGAt                                                                          <  2:270554/149‑1 (MQ=255)
                                                                          tGGAAAACAACCTTGTGTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCTATTGCCTCTGTGAAGGAGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGAt                                                                          <  2:265760/149‑1 (MQ=255)
                                                                           ggAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGAt                                                                          >  8:82119/1‑148 (MQ=255)
                                                                            gAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGAt                                                                          >  6:329646/1‑147 (MQ=255)
                                                                            gAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGAt                                                                          >  1:114221/1‑147 (MQ=255)
                                                                              acacAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAgg                                                                      <  5:89957/149‑1 (MQ=255)
                                                                               cacaACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGt                                                                     >  3:97744/1‑149 (MQ=255)
                                                                                acaACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGt                                                                     >  6:37765/1‑148 (MQ=255)
                                                                                acaACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTc                                                                    <  2:42834/149‑1 (MQ=255)
                                                                                   aCCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGCAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATc                                                                 <  3:282557/149‑1 (MQ=255)
                                                                                   aCCTTGCTTTAACCTGCCACCCGGGTGCCGGGCAGGTTCCGGCGCCGTGCATTGAGCGTCCTGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGTCTCTGCGCCGCACCAGTGCACCGCGCGGCTCGCTGGATAAGGTCATc                                                                 <  6:251717/149‑1 (MQ=255)
                                                                                     cTTGGTTTAACCAGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCAATGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGAACCGCGCGTCTCGCTGGATAAGCTCATCg                                                                >  2:52893/1‑148 (MQ=255)
                                                                                        ggTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaaa                                                             <  2:301364/148‑1 (MQ=255)
                                                                                        ggTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaaa                                                             <  6:53532/148‑1 (MQ=255)
                                                                                        ggTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGaaa                                                             <  3:65485/148‑1 (MQ=255)
                                                                                         gTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATACCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACAAGTGCACCGCGCGTCTCGCTGGATAAGGTCATctca                                                              >  6:84971/1‑143 (MQ=255)
                                                                                            tAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGAt                                                         <  4:61265/148‑1 (MQ=255)
                                                                                             aaCCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATg                                                        >  8:31430/1‑148 (MQ=255)
                                                                                               ccTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTAc                                                     <  5:329646/149‑1 (MQ=255)
                                                                                               ccTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCACACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTAc                                                     >  4:93989/1‑149 (MQ=255)
                                                                                                   cGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGaaa                                                 <  5:12632/149‑1 (MQ=255)
                                                                                                    gACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAAc                                                >  8:189872/1‑149 (MQ=255)
                                                                                                    gACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAAc                                                <  7:210581/149‑1 (MQ=255)
                                                                                                    gACCCGGTTGCAGGGCAGGTTCAGGTGCCGAGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAAc                                                <  5:174638/149‑1 (MQ=255)
                                                                                                       ccGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCgg                                             >  3:136322/1‑149 (MQ=255)
                                                                                                       ccGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCgg                                             <  7:189872/149‑1 (MQ=255)
                                                                                                        cGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGt                                            <  4:65428/149‑1 (MQ=255)
                                                                                                          gTTGCAGGTCATGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTTCCTCTGTGAAGGCGATTAACGCCGCGCTGATTGCTCTGCGCCGCACCAGTGCACCTCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGtta                                          <  8:266257/149‑3 (MQ=255)
                                                                                                          gTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTaa                                          <  7:94371/149‑1 (MQ=255)
                                                                                                           ttGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAg                                         <  2:232796/149‑1 (MQ=255)
                                                                                                             gcagggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAgg                                        <  4:136322/148‑1 (MQ=255)
                                                                                                               agggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACa                                     <  7:227147/149‑1 (MQ=255)
                                                                                                                gggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                    <  2:302846/149‑1 (MQ=255)
                                                                                                                 ggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACa                                     >  5:150924/1‑147 (MQ=255)
                                                                                                                 ggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAAGTCATCGAAACGATGTACGAAACCGGTAAGGACAt                                    >  5:167413/1‑148 (MQ=255)
                                                                                                                 ggcaggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGAGTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACa                                     <  6:150924/147‑1 (MQ=255)
                                                                                                                   caggTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGaa                                 <  4:177330/149‑1 (MQ=255)
                                                                                                                      gTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGc                              >  5:7695/1‑149 (MQ=255)
                                                                                                                       ttCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGc                              >  1:229783/1‑148 (MQ=255)
                                                                                                                             tGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATAc                       <  7:82119/149‑1 (MQ=255)
                                                                                                                               ccGTGAATTGAGCGTAATGCCATTGCCTCTGTGAAGGAGATTAACGCCGCGCGGATGGCTCTGCGACGCACGAGTGCACAGCGCTTCTCGCTGGATAAGGTCATCCAAAGGATGTACGAAACAGGTAAGGGCATGAACGCCAAATACcg                     >  4:110683/1‑149 (MQ=255)
                                                                                                                                cGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACcgc                    <  2:241057/149‑1 (MQ=255)
                                                                                                                                 gTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGCACGCCAAATACcgc                    >  1:160631/1‑148 (MQ=255)
                                                                                                                                         aGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGCACGCCAAATACCGCGAACCCTCa           >  4:155929/1‑149 (MQ=255)
                                                                                                                                             tAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAcgcg       >  5:49382/1‑149 (MQ=255)
                                                                                                                                                tGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtg    >  7:295461/1‑149 (MQ=255)
                                                                                                                                                  ccATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCggtggt  >  6:116215/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCCGTGCATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGT  >  NC_000913/1897975‑1898269

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑