Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,802,656 A→T 10.4% T437S (ACT→TCT)  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,6560AT10.4% 247.7 / 16.5 96T437S (ACT→TCT) nrdEribonucleoside‑diphosphate reductase 2, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (48/38);  new base T (4/6);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

CTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCTTATGG  >  NC_000913/2802508‑2802804
                                                                                                                                                    |                                                                                                                                                    
ctGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGAc                                                                                                                                                     <  1:464346/150‑1 (MQ=255)
   cTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGc                                                                                                                                                  >  2:512560/1‑150 (MQ=255)
     cAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACAATGGATTCCCCCGATTTTGCCAGCACGGTAGAGACTGCCg                                                                                                                                                >  1:66689/1‑150 (MQ=255)
       gAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTg                                                                                                                                              <  1:1446707/150‑1 (MQ=255)
       gAAATTTTGCAGGTTAACAGCGCCTCAGAGAATGACGAGAATCGCGACTATACCCGCACAGGCCATGATATTTCGTGCAATTTAGGTTGGTTGAATATTGCGCACACCATGGATTCCCGCGATTTTGCCCGCACGGTAGAGACTGCCGTg                                                                                                                                              >  2:1625776/1‑150 (MQ=255)
          aTTTTGCAGGTTCACTGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTCAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGCAGTCTCTACCGTgcgg                                                                                                                                           <  1:325491/150‑2 (MQ=255)
          aTTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGCAGTCTCTACCGTgcgg                                                                                                                                           <  2:430441/150‑2 (MQ=255)
          aTTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGCAGTCTCTACCGTgcgg                                                                                                                                           <  1:1410269/150‑2 (MQ=255)
           ttttGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGCAGTCTCTACCGTgcggg                                                                                                                                          <  2:877809/150‑3 (MQ=255)
           ttttGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGCAGTCTCTACCGTgcggg                                                                                                                                          <  1:716657/150‑3 (MQ=255)
            tttGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCgg                                                                                                                                         >  2:478316/1‑150 (MQ=255)
             ttGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGt                                                                                                                                        <  2:731995/150‑1 (MQ=255)
             ttGCAGGTTAACAGCGCCTCAGAGTATGACGAGAAGCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGt                                                                                                                                        <  1:997544/150‑1 (MQ=255)
              tGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGCAGTCTCTACCGTgcgggcaa                                                                                                                                       <  1:34680/150‑6 (MQ=255)
               gCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGttt                                                                                                                                      <  1:17773/150‑1 (MQ=255)
                cAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTa                                                                                                                                     >  2:160912/1‑150 (MQ=255)
                  tgTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCACGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAAc                                                                                                                                   >  1:983503/2‑150 (MQ=255)
                  ggTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAAc                                                                                                                                   <  1:685741/150‑1 (MQ=255)
                  ggTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCACCCGATTTTGCACGCACGGTAGAGACTGCCGTGCGCGGTTTAAc                                                                                                                                   >  2:819805/1‑150 (MQ=255)
                     tAACAGCGCCTCAGAGTATTACGAGTATCTCGAGTATACCGGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGc                                                                                                                                <  1:98502/150‑1 (MQ=255)
                     tAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATTTTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGc                                                                                                                                >  2:1029550/1‑150 (MQ=255)
                       aCAGCGCCTCAGAGTTTGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAg                                                                                                                              <  1:192539/150‑1 (MQ=255)
                        cAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGt                                                                                                                             >  2:1658816/1‑150 (MQ=255)
                            gcCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCTCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCa                                                                                                                         <  2:983503/150‑1 (MQ=255)
                            gcCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCa                                                                                                                         >  1:757111/1‑150 (MQ=255)
                             ccTCAGAGTATGACGAGAATCTCGACTATACACGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAg                                                                                                                        >  1:1713679/1‑150 (MQ=255)
                                cAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGata                                                                                                                     <  2:1648746/150‑1 (MQ=255)
                                  gagTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATg                                                                                                                   >  1:1001244/1‑150 (MQ=255)
                                  gagTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATg                                                                                                                   >  1:1118078/1‑150 (MQ=255)
                                  gagTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATg                                                                                                                   <  2:295241/150‑1 (MQ=255)
                                    gTATGTCGAGAATCACGACTCTACACGCAAAGGCCATGCTATTTCCTTCAATTTAGGTGCGGTGAATATTGCGCACCCCATTTATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGGGCGCGGTTTAACGGCAGTATCATATATGAg                                                                                                                 <  2:126085/150‑1 (MQ=255)
                                    gTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCAGTTTAACGGCAGTATCAGATATGAg                                                                                                                 <  2:1571164/150‑1 (MQ=255)
                                    gTAAGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAg                                                                                                                 <  2:1341808/150‑1 (MQ=255)
                                      aTGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTc                                                                                                               >  2:352963/1‑150 (MQ=255)
                                       tGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATCTTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGGTTAACGGCAGTATCAGATATGAGTCa                                                                                                              <  2:229510/150‑1 (MQ=255)
                                       tGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGTATTCCCCCGATTTTGCCCGCACGGGAGAGACTGCCGTGCGCGGTTTAACGGCAGTAGCAGATATGCGTCa                                                                                                              >  1:329249/1‑150 (MQ=255)
                                       tGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCa                                                                                                              <  2:774576/150‑1 (MQ=255)
                                       tGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCa                                                                                                              <  1:522904/150‑1 (MQ=255)
                                       tGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCa                                                                                                              >  2:1712684/1‑150 (MQ=255)
                                        gACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCat                                                                                                             >  2:158867/1‑150 (MQ=255)
                                          cGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATTTGAGTCatat                                                                                                           >  2:1082523/1‑150 (MQ=255)
                                            agaATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATcc                                                                                                         >  2:10985/1‑150 (MQ=255)
                                             taaTCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTCGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCg                                                                                                        <  1:1333240/149‑1 (MQ=255)
                                                    gACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTg                                                                                                 <  2:65209/150‑1 (MQ=255)
                                                       tataCCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCg                                                                                              >  2:60461/1‑150 (MQ=255)
                                                        ataCCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGt                                                                                             >  2:626999/1‑150 (MQ=255)
                                                          aCCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTcc                                                                                           >  2:1416418/1‑150 (MQ=255)
                                                               cacaGGCCATGATATTTCGTGCTATTTAGGTTCGTTGAATATTGCGCACACCTTGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGa                                                                                      <  2:1161776/150‑1 (MQ=255)
                                                               cacaGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGa                                                                                      >  1:67368/1‑150 (MQ=255)
                                                                acaGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGAATTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGaa                                                                                     <  1:718847/150‑1 (MQ=255)
                                                                   ggCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGcc                                                                                  >  1:1349607/1‑150 (MQ=255)
                                                                   ggCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGcc                                                                                  <  2:1094130/150‑1 (MQ=255)
                                                                     ccATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCgg                                                                                >  1:1324606/1‑150 (MQ=255)
                                                                            atTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATgcc                                                                         >  1:474626/1‑150 (MQ=255)
                                                                               tCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGCCCATCGAAGCCGGAAATgccgcc                                                                      <  2:1032065/150‑1 (MQ=255)
                                                                               tCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATATTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATgccgcc                                                                      >  1:424227/1‑150 (MQ=255)
                                                                                ccTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCt                                                                     >  1:833600/1‑150 (MQ=255)
                                                                                  tGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCg                                                                   <  1:158867/150‑1 (MQ=255)
                                                                                    cAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCCCGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCa                                                                 <  1:1184931/150‑1 (MQ=255)
                                                                                                 ttGAATATTCCGCACACCATGGATTCACCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTg                                                    <  1:819805/150‑1 (MQ=255)
                                                                                                    aaTATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTgggg                                                 >  1:1629433/1‑150 (MQ=255)
                                                                                                    aaTATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTgggg                                                 >  1:392419/1‑150 (MQ=255)
                                                                                                       atTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAg                                              >  1:571243/1‑150 (MQ=255)
                                                                                                           cgcACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGa                                          <  1:1658816/150‑1 (MQ=255)
                                                                                                            gcACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGaa                                         <  2:67368/150‑1 (MQ=255)
                                                                                                             cacaCCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAAt                                        >  1:1666/1‑150 (MQ=255)
                                                                                                             cacaCCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAAt                                        >  2:543415/1‑150 (MQ=255)
                                                                                                               cacCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCCCGCCATCGGACTGGGGCAGATGAAttt                                      >  2:1049549/1‑150 (MQ=255)
                                                                                                               cacCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAAttt                                      <  1:1029550/150‑1 (MQ=255)
                                                                                                                 ccATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTac                                    >  2:582198/1‑150 (MQ=255)
                                                                                                                   aTGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTacac                                  >  2:1540648/1‑150 (MQ=255)
                                                                                                                    tGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACg                                 >  1:1567894/1‑150 (MQ=255)
                                                                                                                    tGGATTCCCCCGATTATGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCCGCGTGCCGTCCAGCGCAGCCGGAAATGACGACTCGCACGCCGTCGGAGTGGGGCTGATGAATTTacaag                                 >  2:355991/1‑148 (MQ=255)
                                                                                                                      gATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGc                               >  1:1240131/1‑150 (MQ=255)
                                                                                                                       aTTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCt                              <  2:1324606/150‑1 (MQ=255)
                                                                                                                       aTTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCt                              <  2:1001244/150‑1 (MQ=255)
                                                                                                                          cccccGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATc                           >  1:1228717/1‑150 (MQ=255)
                                                                                                                            cccGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCCCGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTg                         >  2:1592694/1‑150 (MQ=255)
                                                                                                                              cGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGgc                       <  1:478316/150‑1 (MQ=255)
                                                                                                                              cGAATTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCTATGCCGGAAATGCCGCATCGCACGCCTTCGGACTGGGGCAGATGAATTTACACGGCTATCTGgc                       <  1:1668092/150‑1 (MQ=255)
                                                                                                                               gATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGgcg                      >  1:403194/1‑150 (MQ=255)
                                                                                                                                aTTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGgcgc                     >  2:1362336/1‑150 (MQ=255)
                                                                                                                                 ttttGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGcgcg                    <  2:1349607/150‑1 (MQ=255)
                                                                                                                                   ttGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCgag                  >  2:1495432/1‑150 (MQ=255)
                                                                                                                                    accgcGCACGGCAGTCTCTACCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGATTTTACACGGCTATCTGGCGCgaga                 >  1:5881/5‑150 (MQ=255)
                                                                                                                                    accgcGCACGGCAGTCTCTACCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCgaga                 >  2:1580198/5‑150 (MQ=255)
                                                                                                                                    accgcGCACGGCAGTCTCTACCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCgaga                 >  2:1121572/5‑150 (MQ=255)
                                                                                                                                    accgcGCACGGCAGTCTCTACCGTGCGCGGTTTAACGGCAGTATCAGATAAGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCgaga                 >  1:478782/5‑150 (MQ=255)
                                                                                                                                      cccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAg               >  1:371815/1‑150 (MQ=255)
                                                                                                                                      cccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAg               <  1:160912/150‑1 (MQ=255)
                                                                                                                                      cccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAg               >  1:856241/1‑150 (MQ=255)
                                                                                                                                      cccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAg               >  1:850584/1‑150 (MQ=255)
                                                                                                                                      cccGCACGGTAGAGACTGCCGTGAGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGTACTGGGGCAGATGAATTTACACGGCTATCTGCCGCGAGCAg               >  1:1703807/1‑150 (MQ=255)
                                                                                                                                       ccGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAgg              <  2:403194/150‑1 (MQ=255)
                                                                                                                                        cGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGc             <  1:1335915/150‑1 (MQ=255)
                                                                                                                                           aCGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATc          <  1:352963/150‑1 (MQ=255)
                                                                                                                                                  agaCTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCTTATg   <  2:1629433/150‑1 (MQ=255)
                                                                                                                                                   gaCTGCCGGGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCTTATgg  <  2:534550/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCTTATGG  >  NC_000913/2802508‑2802804

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: