Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,802,655 G→C 10.4% E436D (GAG→GAC) ‡ 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,6550GC10.4% 299.5 / 22.5 96E436D (GAG→GAC) ‡nrdEribonucleoside‑diphosphate reductase 2, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (48/38);  new base C (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 = 9.32e-01

CTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCTTATGG  >  NC_000913/2802506‑2802804
                                                                                                                                                     |                                                                                                                                                     
ctctTCTCAGAAATTTTGCACGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTagag                                                                                                                                                       <  1:826931/150‑1 (MQ=255)
  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                                                                                                                   <  2:295241/150‑1 (MQ=255)
                                    gagTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATg                                                                                                                   >  1:1118078/1‑150 (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:1712684/1‑150 (MQ=255)
                                         tGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCa                                                                                                              <  2:774576/150‑1 (MQ=255)
                                         tGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCa                                                                                                              <  1:522904/150‑1 (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:392419/1‑150 (MQ=255)
                                                                                                      aaTATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTgggg                                                 >  1:1629433/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:1121572/5‑150 (MQ=255)
                                                                                                                                      accgcGCACGGCAGTCTCTACCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCgaga                 >  2:1580198/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)
                                                                                                                                                     |                                                                                                                                                     
CTCTGCTCAGAAATTTTGCAGGTTAACAGCGCCTCAGAGTATGACGAGAATCTCGACTATACCCGCACAGGCCATGATATTTCCTGCAATTTAGGTTCGTTGAATATTGCGCACACCATGGATTCCCCCGATTTTGCCCGCACGGTAGAGACTGCCGTGCGCGGTTTAACGGCAGTATCAGATATGAGTCATATCCGCAGCGTGCCGTCCATCGAAGCCGGAAATGCCGCCTCGCACGCCATCGGACTGGGGCAGATGAATTTACACGGCTATCTGGCGCGAGAAGGCATCGCTTATGG  >  NC_000913/2802506‑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: