Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,869,711:1 +C 8.1% intergenic (‑28/+24) thyA ← / ← lgt thymidylate synthase/prolipoprotein diacylglyceryl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,869,7111.C 8.1% 662.1 / 52.3 160intergenic (‑28/+24)thyA/lgtthymidylate synthase/prolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base . (67/80);  new base C (5/8);  total (72/88)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 167 total aligned reads displayed.

CCTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTT  >  NC_002947/5869563‑5869859
                                                                                                                                                     |                                                                                                                                                    
ccTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAATCGGCAggggg                                                                                                                                                     >  2:3720638/1‑150 (MQ=255)
ccTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAggggg                                                                                                                                                     >  2:1073299/1‑150 (MQ=255)
 cTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGa                                                                                                                                                    >  2:3914893/1‑150 (MQ=255)
 cTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGa                                                                                                                                                    >  2:3510805/1‑150 (MQ=255)
    cTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAATACGGCAGGGGGATCg                                                                                                                                                 >  2:4976461/1‑150 (MQ=255)
     tGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGc                                                                                                                                                >  2:398953/1‑150 (MQ=255)
        aGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAACGGCAGGGGCGATCcccct                                                                                                                                             <  1:4247507/150‑2 (MQ=255)
        aGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCAGGGCTTACCTGTAGCAAAAACGGCAGGGGCGATCcccc                                                                                                                                              <  2:1426229/150‑1 (MQ=255)
           tCGAAACGCAGCATCGCGCCCGGCTGGCTGATGGGGCGGATGCCGGTGCGTTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCTCCTGc                                                                                                                                          <  2:1453755/150‑1 (MQ=255)
           tCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATCCCCCTGc                                                                                                                                          >  1:3062654/1‑150 (MQ=255)
           tCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGAGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGc                                                                                                                                          >  2:3100379/1‑150 (MQ=255)
            cGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGGTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATCCCCCTGcc                                                                                                                                         <  2:4675477/150‑1 (MQ=255)
              aaaCGCCGCATCGCGCCGGGCAGGGTGTTGGTGCGTATGCCGGTGCGGTTTCCCTGCACCTTGCCGTTATAAATGACTGCGCGGACCTGGTCCAGATACTTTTTCTGGGCTTACCTTTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGt                                                                                                                                       <  1:2053730/150‑1 (MQ=255)
              aaaCGCAGCCTCGCGCCGGGCAGGCTGTTGGTGCGGATGCGGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGt                                                                                                                                       <  1:2053538/150‑1 (MQ=255)
              aaaCGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGt                                                                                                                                       <  2:3279573/150‑1 (MQ=255)
              aaaCGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGt                                                                                                                                       >  2:3591020/1‑150 (MQ=255)
               aaCGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTg                                                                                                                                      <  2:869310/150‑1 (MQ=255)
                aCGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGTTTTCCCTGAAACGTGCCTTTATCAATGACGTCGCGGACCAGTTCAATATACTGTTTCATGGCTTACCTTTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTgg                                                                                                                                     <  2:2266514/150‑1 (MQ=255)
                aCGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTgg                                                                                                                                     <  2:4063942/150‑1 (MQ=255)
                aCGCAGCACCGCGCCCGGCAGGCTGATGGTGCGGATGCGGGTGCGGTTTCCCTGCACCGTGCGGTTACCAATGACGTCGCGCACCTGGTCCCGAGCCTGTTTCTTGGCTTACCTGTAGCAAACACGGCAGGGGGATCGCCCCTGCCGTgg                                                                                                                                     <  1:1531252/150‑1 (MQ=255)
                  gccgccgcgcgCCCGGCCGGCGGTGGGTGCGGAGGCGGGTGCGGTTCCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCCGGCCCAGATACTGTTTCCTGGCTTACCGGTCGCAAACACGGCAGGGGGATCGCCCCTGCCGTGGGa                                                                                                                                   <  1:543579/143‑1 (MQ=255)
                  gcagcaTCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGa                                                                                                                                   <  1:5449364/150‑1 (MQ=255)
                  gcagcaCCGCCCCGGGCTGGCTGTTGGTGGGGTTGCGGGTGGTGTTTCCCTGCACCGTGCCTTTATAAATGAGGGCGCGGACCCGGTCCAGATACTGTTTCAGGGCTAACCTGTCGCAAACACGGCAGGGGGATCGCCCCTGCCGTGGGa                                                                                                                                   <  1:1416487/150‑1 (MQ=255)
                  gcaccaTCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTCCCCTGCAACGCGCGGTTATCAATGACGTCGCGGACCTGGACCAGATACTGTTTCATGGCTTACCTGTAGCAAACACGGCAGGGGGATCGCCCCTGCCGTGGGa                                                                                                                                   <  1:970193/150‑1 (MQ=255)
                   cagcaTCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATCCCCCTGCCGtttttg                                                                                                                                  <  1:2248524/150‑6 (MQ=255)
                     gcaTCGCGCCCGGCAGGCTGATGGTGCGGATGCGGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCCGGTCCATATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTa                                                                                                                                <  2:2472823/150‑1 (MQ=255)
                        tCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATc                                                                                                                             >  1:2475309/1‑150 (MQ=255)
                        tCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATc                                                                                                                             >  2:1150781/1‑150 (MQ=255)
                        tCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAAGGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATc                                                                                                                             >  2:1151107/1‑150 (MQ=255)
                            gcCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCTTGGCTTACCTGTAGAAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGc                                                                                                                         <  2:2524893/150‑1 (MQ=255)
                            gcCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGc                                                                                                                         <  2:5250881/150‑1 (MQ=255)
                                ggCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGt                                                                                                                     >  1:790108/1‑150 (MQ=255)
                                 gCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGtt                                                                                                                    >  1:2497062/1‑150 (MQ=255)
                                     gCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGc                                                                                                                >  2:5318027/1‑150 (MQ=255)
                                        gATGGTGCGTATGCCGGGGCGGTTGCCCTGCAACGTGCCTTTATCAATGACGTGGCGGACCTGGTCCAGATACTGTTTCTTGGCTTACCTTTAGAAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGt                                                                                                             <  1:571456/150‑1 (MQ=255)
                                           ggTGCGGATGCGGGTGCGGTTTCCCTGCAACGTGCCTTTATCATTGACGTCGCGGACCTGGTCCAGATACTGTTTCTGGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAgg                                                                                                          <  1:4287789/150‑1 (MQ=255)
                                            gTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGTCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGc                                                                                                         >  1:63442/1‑150 (MQ=255)
                                             tGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCt                                                                                                        >  2:3378281/1‑150 (MQ=255)
                                             tGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCt                                                                                                        >  2:1904047/1‑150 (MQ=255)
                                                ggATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCCGGTCCAGATACTGTTTCATGGCTTACCTGTATCAAAAACGGCAGGGGGGAGGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGc                                                                                                     <  1:3510805/150‑1 (MQ=255)
                                                 gATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCg                                                                                                    >  2:5386699/1‑150 (MQ=255)
                                                      cGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATtat                                                                                               >  2:1800351/1‑150 (MQ=255)
                                                      cGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACGGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATGAATCAGGCGGGTTTCGCCGTAGGCTTGCGATtat                                                                                               <  2:762080/150‑1 (MQ=255)
                                                        gTGCGGTTTCCCTGCAACGTGCCTTTATAAATGACTTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACATGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                             <  1:3805971/150‑1 (MQ=255)
                                                        gTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCTTGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                             <  2:595157/150‑1 (MQ=255)
                                                        gTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                             <  1:1338947/150‑1 (MQ=255)
                                                        gTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCAGGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                             <  2:279800/150‑1 (MQ=255)
                                                        gTGCGGTTGCCCTGCAACGTGCCGTTATCACTGCCGTCGCGGACCAGGGCCAGATCGTGTTTCAGGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                             <  1:2585835/150‑1 (MQ=255)
                                                             gTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCcca                                                                                        >  1:3299901/1‑150 (MQ=255)
                                                                     gCAACGTCCCGTTATCAATGACGTCGCGGACCAGCTCCAGATACTGTTTCTTGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCa                                                                                <  2:169497/150‑1 (MQ=255)
                                                                      cAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCAAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCaa                                                                               <  1:2320952/150‑1 (MQ=255)
                                                                           tGCGGTTATCAATGCCGTCGGGGACCCGGTCCTGCTCCGGTTTCCGGGCTTACGTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGa                                                                          <  1:1073299/150‑1 (MQ=255)
                                                                           tGCCGTTATCATTGACGTCGCGTACCTGGTCCAGATACTTTTTCTGGGCTTACCTTTAGAAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGa                                                                          <  1:215706/150‑1 (MQ=255)
                                                                           tGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCTTGGCTTACCTGTAGAAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGa                                                                          <  2:517130/150‑1 (MQ=255)
                                                                            gCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAg                                                                         <  2:3466047/150‑1 (MQ=255)
                                                                             cctttATAAATGACGTCGCGGACCCGGTCCTGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAgg                                                                        <  1:2395970/147‑1 (MQ=255)
                                                                                    cAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTGGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGa                                                                 >  1:2324436/1‑150 (MQ=255)
                                                                                      aTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGcc                                                                                                                                         >  2:2250298/1‑76 (MQ=255)
                                                                                        gACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTat                                                             >  2:3570412/1‑150 (MQ=255)
                                                                                           gTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcat                                                          >  1:593475/1‑150 (MQ=255)
                                                                                           gTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcat                                                          >  1:3348246/1‑150 (MQ=255)
                                                                                            tCGCGGACCAGTTCAATATACTTTTTAATGGCTTACCTGTAGAAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcatc                                                         <  1:3077787/150‑1 (MQ=255)
                                                                                            tCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcatc                                                         >  1:1551540/1‑150 (MQ=255)
                                                                                            tCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcatc                                                         <  1:2136176/150‑1 (MQ=255)
                                                                                            tCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcatc                                                         >  1:909308/1‑150 (MQ=255)
                                                                                             cgcgGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCg                                                        >  1:3141219/1‑150 (MQ=255)
                                                                                                 gACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAACACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAAt                                                    <  1:3756659/150‑1 (MQ=255)
                                                                                                    cAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGca                                                 <  1:1163006/150‑1 (MQ=255)
                                                                                                       gTCCAGATACTGTTTCTTGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAg                                              <  2:4152156/150‑1 (MQ=255)
                                                                                                        tCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGc                                             <  2:2993445/150‑1 (MQ=255)
                                                                                                          cAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCaa                                           <  1:2887367/150‑1 (MQ=255)
                                                                                                           aGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAAc                                          >  1:3481165/1‑150 (MQ=255)
                                                                                                            gATACTGTTTCATGGCTTACCTGTAGCAAAATCGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACt                                         >  2:2405658/1‑150 (MQ=255)
                                                                                                              tACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTgg                                       <  1:398953/150‑1 (MQ=255)
                                                                                                               aCTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGc                                      <  1:399277/150‑1 (MQ=255)
                                                                                                               aCTGTTTCATGGATTACGTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGc                                      <  1:5018127/150‑1 (MQ=255)
                                                                                                                  gTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCtt                                   <  1:3958784/150‑1 (MQ=255)
                                                                                                                   tttCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtg                                  <  1:1668149/150‑1 (MQ=255)
                                                                                                                   tttCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtg                                  >  1:370170/1‑150 (MQ=255)
                                                                                                                   tttCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtg                                  >  2:18111/1‑150 (MQ=255)
                                                                                                                    ttCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtga                                 <  1:4516747/150‑1 (MQ=255)
                                                                                                                    ttCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtga                                 <  1:1788813/150‑1 (MQ=255)
                                                                                                                    ttCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtga                                 <  2:4278099/150‑1 (MQ=255)
                                                                                                                         ggCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtgagtgag                            >  1:1886739/1‑150 (MQ=255)
                                                                                                                               ccTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGc                                                                                           >  2:4925305/1‑81 (MQ=255)
                                                                                                                               ccTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGc                                                                                           <  1:4925305/81‑1 (MQ=255)
                                                                                                                               ccTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGagccagcc                      <  1:2637108/150‑1 (MQ=255)
                                                                                                                                  gTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGaa                   >  2:1754743/1‑150 (MQ=255)
                                                                                                                                   tAGCAAACACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAc                  <  1:1040298/150‑1 (MQ=255)
                                                                                                                                       acacaCGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa              <  1:5268905/146‑1 (MQ=255)
                                                                                                                                          caCGGCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAcgc                >  2:5265558/2‑145 (MQ=255)
                                                                                                                                           aCGGCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTa          >  2:3400831/1‑150 (MQ=255)
                                                                                                                                            cGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTaa         <  1:1800351/150‑1 (MQ=255)
                                                                                                                                             ggCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAAc        >  1:4136439/1‑150 (MQ=255)
                                                                                                                                                aGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGTGGCCACCGACTATCATTGGAATGCACAGCAACTGGCATTGAGTGAGCCAGCCGAACGCGATGTAACCGa     >  1:1340658/1‑150 (MQ=255)
                                                                                                                                                 gggggATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAg    <  1:2479260/150‑1 (MQ=255)
                                                                                                                                                  ggggATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGt   >  1:3190318/1‑150 (MQ=255)
                                                                                                                                                   gggATCGCCCCTGCCGTGGGATTAATCAGGCGGTTTTCGCCGTAGGATTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGtt  <  2:3849369/150‑1 (MQ=255)
                                                                                                                                                   gggATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCaa                                           <  2:263228/109‑1 (MQ=255)
                                                                                                                                                   gggATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCaa                                           >  1:263228/1‑109 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTT  >  NC_002947/5869563‑5869859

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: