Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,802,654 A→T 5.1% E436V (GAG→GTG)  nrdE → ribonucleoside‑diphosphate reductase 2, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,802,6540AT 5.1% 243.1 / 2.1 85E436V (GAG→GTG) nrdEribonucleoside‑diphosphate reductase 2, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (36/44);  new base T (2/2);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

CTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCAT  >  NC_000913/2802508‑2802795
                                                                                                                                                  |                                                                                                                                             
ctGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAgaga                                                                                                                                             >  2:1378554/1‑149 (MQ=255)
 tGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACACGCACAGGCAATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGAGCACACCATGGATTACCCCGATTTTGCCAGGACGGTAGGGAc                                                                                                                                            >  2:1003187/1‑149 (MQ=255)
      aGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCg                                                                                                                                       >  2:1276011/1‑149 (MQ=255)
        aaaTTTTGCAGGTTAACAGCGCGTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGAAATTTCCTGCAATTTAGGATCGTTGAATATTGCGCACACCATGTATTCCACCGATTTTGCACGAACGGTAGATACAGCAGTg                                                                                                                                     >  1:965426/1‑149 (MQ=255)
         aaatttGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTgc                                                                                                                                    <  2:1203231/146‑1 (MQ=255)
           ttttGCTGGTTAACAGCGCCTCATAGTATGACGAGAATCTCGACTATACCCGAACATGCCATGATATTTCCTGCAATTTAGGTTCGTTGATTATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGCAGTCTCTACCGTgcgg                                                                                                                                  <  2:393664/149‑2 (MQ=255)
           ttttGCAGGTTAACAGCGCCTCAGAGTATGCCGAGAATCTCGACTATAGCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGCAGTCTCTACCGTgcgg                                                                                                                                  <  2:211645/149‑2 (MQ=255)
           ttttGCAGGTTAACAGCGCCGCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTgcgc                                                                                                                                  <  1:1046068/149‑1 (MQ=255)
                cAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGttt                                                                                                                             >  2:837330/1‑149 (MQ=255)
                cAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGGGCGCGGttt                                                                                                                             >  2:474613/1‑149 (MQ=255)
                 tggTTAACAGCCCCTCAGAGTATGACGCGAATCTCGACAATACCCGCACAGGCCAAGATATTTCCTGCAATAAAGATTCGATGAAAATTGCGCACACAAAGGATTCCCCCGAATTTTCCCGCACGGTAGAGACTGCCGTGCGCCGTTTa                                                                                                                            >  2:1076880/2‑149 (MQ=255)
                   gTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAAc                                                                                                                          <  1:556529/149‑1 (MQ=255)
                   gTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCAGCACGGTAGAGACTGCCGTGCGCGGTTTaag                                                                                                                          >  2:1273716/1‑148 (MQ=255)
                     tAACAGCTCATCAGCTGATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACgg                                                                                                                        <  1:1148619/149‑1 (MQ=255)
                     tAACAGCGCCTCAGAGTATGACGAGAATCTCGATTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACgg                                                                                                                        <  2:589062/149‑1 (MQ=255)
                      aaCAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGc                                                                                                                       <  1:546208/149‑1 (MQ=255)
                      aaCAGCGCCTCAGAGTATGACCAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGc                                                                                                                       <  2:249463/149‑1 (MQ=255)
                       aCAGCGCCTGAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCa                                                                                                                      <  1:252859/149‑1 (MQ=255)
                       aCAGCGCCTCAGCGTATGTCGCTAATCGCGACTATTACCGCACAGTCCATGATATTTCCTCCAATTTAGGTTCGTTGAATATTGCGCACTCCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCa                                                                                                                      <  2:577494/149‑1 (MQ=255)
                           ggcCTCAGAGTATGCCGAGAATCTCGACATTACCCGCACAGGCCATGATATTTCCTGCACTTTAGGTTCGTTGAATCTTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTAt                                                                                                                  <  1:1558330/148‑1 (MQ=255)
                             ccTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGACATGATATTTCCTACAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCa                                                                                                                <  2:86792/149‑1 (MQ=255)
                               tCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTCCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGa                                                                                                              <  1:1348967/149‑1 (MQ=255)
                                cAGAGTAAGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGat                                                                                                             <  1:1447598/149‑1 (MQ=255)
                                  gagTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGatat                                                                                                           <  1:1378554/149‑1 (MQ=255)
                                    gTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGa                                                                                                         <  2:814718/149‑1 (MQ=255)
                                    gTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCAGGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGa                                                                                                         >  2:899923/1‑149 (MQ=255)
                                     tATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAg                                                                                                        <  2:1398267/149‑1 (MQ=255)
                                      atgcagtgaATCTCCTATAAACCCGCACATTCCATGATTTTTCCTTCAAATTAGGTTCGTTGAATATTGCGCACACCATGGAATCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGt                                                                                                       <  1:1078080/142‑1 (MQ=255)
                                      aTGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGt                                                                                                       <  2:546636/149‑1 (MQ=255)
                                      aTGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATAATGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGt                                                                                                       <  2:561102/149‑1 (MQ=255)
                                          cGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCata                                                                                                   >  1:377923/1‑149 (MQ=255)
                                             taaTCTCGACTATACCCGCACAGGCAATGATCTTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCTTGGTTAGCCCCGCTTTTGCCCGCACGGTAGACACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATcc                                                                                                <  2:434865/148‑1 (MQ=255)
                                             gagtctcGACTATACCCGCACAGGCCAAGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCCCGGTAGAGCCTGCCGTGCGCGGTTTTACGGCAGTATCAGATATGAGTCAAATcc                                                                                                >  2:383459/4‑149 (MQ=255)
                                                tctcGACTATACCCGCCCAGGCCATGATATTTCCTGCAATTTTGGTTCGTTGAATATTGCGCACACCATGGAGTCCCCCGATTTTGCCCGCCCGGTAGAGTCTGCCGTGCGCGGTGTAACGGCAGTATCAGATTTGAGTCATATCCGcc                                                                                             >  2:493473/1‑148 (MQ=255)
                                                cctcGACTATACCCGCACCGGCCATGATATTTCCAGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGCTTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGGGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCa                                                                                             <  2:930166/148‑1 (MQ=255)
                                                  tcGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATTTGAGTCATATCCTCAGc                                                                                           >  1:1400728/1‑149 (MQ=255)
                                                  tcGACTAAACCCGCACAGGCCATGATATTACCTGCAATTTAGGATCGTTGAATATTGAGCACACCATGCCTTCCCCCGATTTAGCCCGCACGGTAGAGACTGCCATGCGCTGTATAACGGCAGTATCAGATATGAGTCATATCCCCAGc                                                                                           >  2:801456/1‑149 (MQ=255)
                                                    gACTATACCCGCAGAGGCCATGATATTTCCTGCAATATAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGt                                                                                         <  2:486746/149‑1 (MQ=255)
                                                     aCTACACCCGCACGTGCGATGTTCGTTCCTGCAAGTAAGGTTCGTGGAATATTGCGCCCACCATGGAGTCCGCCGATATTGCCGGGACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTg                                                                                        <  1:1086257/149‑1 (MQ=255)
                                                         taCCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTATAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGt                                                                                    <  2:1260973/149‑1 (MQ=255)
                                                           cccGCACAGGTCATGATATTTCCTGCAATTTAGGTTCGTTTAATATTGCGCTCACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTcc                                                                                  <  1:1133242/149‑1 (MQ=255)
                                                            ccGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCCGCGTGCCGTCCa                                                                                 >  1:747203/1‑149 (MQ=255)
                                                            ccGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCa                                                                                 <  2:34197/149‑1 (MQ=255)
                                                            ccGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGGATCAGAGATGAGTCATATCCGCAGCGTGCCGTCCa                                                                                 >  2:1483086/1‑149 (MQ=255)
                                                              gCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATc                                                                               >  2:1168453/1‑149 (MQ=255)
                                                                 caGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGaa                                                                            <  1:444877/149‑1 (MQ=255)
                                                                  aGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCCCGGTAGAGACTGCCGTGCGAGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAg                                                                           >  2:1060189/1‑149 (MQ=255)
                                                                    gCCATGTTATTTCCTGCAATTTATGTTCGTTGAATTTTGCGTACACCATGGTTTCGCCCGATTTTGCCCGAACGGTTGAGACTGCCGTGCGCGGTTTAAAGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGcc                                                                         <  1:1050409/149‑1 (MQ=255)
                                                                    gCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAAGGGCAGTATCAGATATGAGTCATATACGCAGCGTGCCGTCCATCGAAGcc                                                                         >  1:795878/1‑149 (MQ=255)
                                                                    gCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGcc                                                                         >  2:946613/1‑149 (MQ=255)
                                                                          atatTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAAt                                                                   >  1:799599/1‑149 (MQ=255)
                                                                             tttCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATgcc                                                                >  1:1398734/1‑149 (MQ=255)
                                                                              ttCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATgccg                                                               >  1:825062/1‑149 (MQ=255)
                                                                                 cTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCt                                                            >  1:636850/1‑149 (MQ=255)
                                                                                  tGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTc                                                           <  1:1195591/149‑1 (MQ=255)
                                                                                    cAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAGCGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGc                                                         <  1:246634/149‑1 (MQ=255)
                                                                                      aTTTAGGTGAGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCAc                                                       <  1:891576/149‑1 (MQ=255)
                                                                                         tAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCACGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGcc                                                    <  1:566556/149‑1 (MQ=255)
                                                                                            gTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGACGCCGGAAATGCCGCCTCGCACGCCATc                                                 >  2:651737/1‑149 (MQ=255)
                                                                                            gTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATc                                                 <  1:726987/149‑1 (MQ=255)
                                                                                                  tGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTg                                           <  1:709203/149‑1 (MQ=255)
                                                                                                   gAATAGAGCCCACACCATGGATTCCCCCTAGATTGCCCCCCCGGTAGAGACTGCAGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGGGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACggg                                          <  2:1027389/149‑1 (MQ=255)
                                                                                                    aaTATTGCTCACACGATGGATTCGCCCGATTTTGCCCGCAGGGTAGAGACTTCCGTGCGCGGTTTAACGGCAGTATCAGATATGATTGATATAAGCAGAGTGCCGTCCATCGAAGCCGGAAATGCCGCCTAGCACGCGATGGGACTggg                                         >  2:865636/1‑149 (MQ=255)
                                                                                                     atatTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTgggg                                        <  1:521743/149‑1 (MQ=255)
                                                                                                        ttGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAg                                     <  1:946613/149‑1 (MQ=255)
                                                                                                         ggcgcACACCATGGATTCCCCCGATTTAGCCCCCACGGTAGAGACTGCCGGGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGGGCCGGCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGa                                    <  1:77931/148‑1 (MQ=255)
                                                                                                            gcACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGa                                 <  1:1261409/149‑1 (MQ=255)
                                                                                                            gcACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGa                                 >  1:1341670/1‑149 (MQ=255)
                                                                                                                    tGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTacac                         <  1:563925/149‑1 (MQ=255)
                                                                                                                        ttCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCt                     >  2:576883/1‑149 (MQ=255)
                                                                                                                         tCCCCCGATTTTCCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATCTGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATTAATTTACACGGCTa                    <  2:1550884/149‑1 (MQ=255)
                                                                                                                            cccGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCt                 >  2:195461/1‑149 (MQ=255)
                                                                                                                            cccGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCATCGCACGCCATCGGACTGGGGCAGATCAATTTACACGGCTATCt                 >  1:432326/1‑149 (MQ=255)
                                                                                                                             ccGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACTGCAGTATCAGATAGGAGTCATATCCGCAGCGTGCCGTCCAACGAAGCCGGAAATGCCTCATCGCACGCCCTCGGCCTGGGGCAGAGGAATTTACACGGCTATCTg                >  1:1552355/1‑149 (MQ=255)
                                                                                                                             ccGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTg                >  2:807371/1‑149 (MQ=255)
                                                                                                                              cGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTgg               <  2:795878/149‑1 (MQ=255)
                                                                                                                                      cgcGCACGGCAGTCTCTACCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGaa       >  2:1424811/3‑149 (MQ=255)
                                                                                                                                      cgcGCACGGCAGTCTCTACCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCAGCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGaa       >  1:1389670/3‑149 (MQ=255)
                                                                                                                                       ccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAg      >  1:396190/1‑149 (MQ=255)
                                                                                                                                       ccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAg      >  2:76632/1‑149 (MQ=255)
                                                                                                                                       ccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAg      >  1:1188047/1‑149 (MQ=255)
                                                                                                                                       ccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCCACGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTTCACGGCTATCTGGCGCgagacg      >  1:399754/1‑147 (MQ=255)
                                                                                                                                       ccGCACGGTAGAGACTGCCGTGCGCGGGTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGGCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGCAGTTACACGGCTAGCTGGCGCGGGAAg      >  1:1157704/1‑149 (MQ=255)
                                                                                                                                           aCGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCAt  <  2:303612/149‑1 (MQ=255)
                                                                                                                                           aCGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATAGCCGCAGCGTGCCGTCCATCGAAGCCGGAATTGCCGCCGCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGctt  <  1:899923/149‑3 (MQ=255)
                                                                                                                                           aCGGTAGAGACTGCCGTGCGAGGTTTAACGTCAGTATCAGATATGAGTCATATCCGCAGCGAGCCGTGAATCGAAGCCGGAAATGCCGCCTCGCCAGCCATCGGACTTCGGCAGATGCATTTACACCGCTATCTGGCTCGCGCAGGCAt  >  1:875241/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCAT  >  NC_000913/2802508‑2802795

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