Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,867,733:1 +C 7.3% 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,867,7331.C 7.3% 314.7 / 19.2 82intergenic (‑28/+24)thyA/lgtthymidylate synthase/prolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base . (35/41);  new base C (3/3);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-02

CCTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTT  >  NC_002947/5867585‑5867881
                                                                                                                                                      |                                                                                                                                                    
ccTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGACGTTATCAATGACGTCTCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAgggg                                                                                                                                                       >  1:109958/1‑149 (MQ=255)
 cTTCTGCAGGTCGAAACGCTGCATCGCGCCCGGCTGGCTGATGGTGCGGATGCCGGGGCGGTTTCCCTGAAACGTGCGGTTATAAATGACGTCGCGGACCTGGTCCAGATACTTTTTCATGGCTTACCTGTAGCAAAAACGGCAggggg                                                                                                                                                     <  4:2999/149‑1 (MQ=255)
    cTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGAt                                                                                                                                                   <  3:180541/148‑1 (MQ=255)
       cAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCCTGTAGCAAAAACGGCAGGGGGATCGc                                                                                                                                                <  1:196362/149‑1 (MQ=255)
       cAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGcc                                                                                                                                               >  3:303574/1‑149 (MQ=255)
             gAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTTCCGTTATCAATGACGTCGCGCACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATCCCCCTGc                                                                                                                                          <  2:165943/148‑1 (MQ=255)
             gAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATCCCCCTGcc                                                                                                                                         <  7:353596/149‑1 (MQ=255)
                  gcagcaTCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTgg                                                                                                                                     >  8:288640/1‑148 (MQ=255)
                    agcaTCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGAt                                                                                                                                  >  6:320052/1‑149 (MQ=255)
                    agcaTCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGAt                                                                                                                                  >  1:305411/1‑149 (MQ=255)
                    agcaTCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGAGACTGATTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGAt                                                                                                                                  >  1:310046/1‑149 (MQ=255)
                        tCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAAt                                                                                                                              <  7:208496/149‑1 (MQ=255)
                         cgcgCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAAt                                                                                                                              >  2:88156/1‑148 (MQ=255)
                          gcgcCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAAt                                                                                                                              >  2:188183/1‑147 (MQ=255)
                          gcgcCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCa                                                                                                                            <  7:41691/149‑1 (MQ=255)
                             cccGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGCCGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGc                                                                                                                         <  2:137392/149‑1 (MQ=255)
                              ccGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTCTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCg                                                                                                                        <  1:88156/149‑1 (MQ=255)
                               gggCAGGCTTTTGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGCACCAGGTCAAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCgg                                                                                                                       <  4:5144/148‑1 (MQ=255)
                               cGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCgg                                                                                                                       <  2:163064/149‑1 (MQ=255)
                                 gCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGt                                                                                                                     <  3:169996/149‑1 (MQ=255)
                                     gCTGTTGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCg                                                                                                                 <  6:73428/149‑1 (MQ=255)
                                        gATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATAGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCg                                                                                                              >  3:191626/1‑149 (MQ=255)
                                          tGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTa                                                                                                            <  1:188183/149‑1 (MQ=255)
                                             tGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGc                                                                                                         >  1:294597/1‑149 (MQ=255)
                                               cGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCtt                                                                                                       <  4:303574/149‑1 (MQ=255)
                                                        gTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATtat                                                                                               >  4:169668/1‑148 (MQ=255)
                                                         tGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATaa                                                                                             <  1:190582/149‑1 (MQ=255)
                                                             gTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGccc                                                                                         <  5:320052/149‑1 (MQ=255)
                                                               ttCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCccac                                                                                       >  7:208375/1‑149 (MQ=255)
                                                                  ccTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACATCAAGGCGATCGCCCCTGCCGTGGGATAAATAAGCCGGGTTTCGCCGTAGGATTGCTATTATAAGAccac                                                                                       >  4:198/1‑146 (MQ=255)
                                                                  ccTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGATTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCGGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCccac                                                                                       >  3:272897/1‑146 (MQ=255)
                                                                    tGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGa                                                                                   >  1:300014/1‑148 (MQ=255)
                                                                      cAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATc                                                                                 <  3:253703/148‑1 (MQ=255)
                                                                       aaCGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCaa                                                                               <  3:346391/149‑1 (MQ=255)
                                                                       aaCGTGCCGTTATCAATCACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCa                                                                                >  1:172675/1‑148 (MQ=255)
                                                                             ccGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAg                                                                         <  8:349389/149‑1 (MQ=255)
                                                                              cGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAgg                                                                        >  5:75838/1‑149 (MQ=255)
                                                                                ttATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGAAAGCCCCTGCCGTGGGATTAATAAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCCCAAGATCATCACGAGGc                                                                       >  8:353752/1‑148 (MQ=255)
                                                                                    cAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATtat                                                                                               <  8:82184/120‑1 (MQ=255)
                                                                                    cAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATtat                                                                                               >  7:82184/1‑120 (MQ=255)
                                                                                     aaTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGa                                                                 >  5:180759/1‑149 (MQ=255)
                                                                                     aaTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGa                                                                 <  3:198/149‑1 (MQ=255)
                                                                                     aaTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGa                                                                 >  6:55258/1‑149 (MQ=255)
                                                                                       tGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACt                                                               <  2:300014/149‑1 (MQ=255)
                                                                                        gACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGGGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTa                                                              >  2:72373/1‑149 (MQ=255)
                                                                                         aCGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAAt                                                                                                                              <  4:260933/84‑1 (MQ=255)
                                                                                         aCGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAAt                                                                                                                              >  3:260933/1‑84 (MQ=255)
                                                                                         aCGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTat                                                             <  7:288640/149‑1 (MQ=255)
                                                                                         aCGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTat                                                             <  8:208375/149‑1 (MQ=255)
                                                                                         aCGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTa                                                              <  3:169668/148‑1 (MQ=255)
                                                                                           gTCGCGGACCAGTACCAGATACTGTATCATGGCTTACCTGTAGCAAAAACTGCAGGGGGATCTCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCACATCAAGACGAGGCCACCGACTatca                                                           <  6:377555/149‑1 (MQ=255)
                                                                                            tCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcat                                                          >  5:151120/1‑149 (MQ=255)
                                                                                            tCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcat                                                          >  6:50436/1‑149 (MQ=255)
                                                                                             cgcgGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGAGCGCCCCTGCCGTGGGATTAATCAGGCGGGGTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcatc                                                         <  7:249371/149‑1 (MQ=255)
                                                                                              gcgGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCg                                                        <  2:122357/149‑1 (MQ=255)
                                                                                              gaggACCAGGTCCAGATAATGTTTGATTGCTTACCTTTAGCAAAAACTGCAGGGGGATAGCCCCTGCCGTGGGATTAATCAGGCTGGTTTCGCCGTAGGCTTGCGATTATAAGCACACCAGATCAAGACGAGGCCACCGACTatcatc                                                         <  2:109958/146‑1 (MQ=255)
                                                                                                   ccAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATg                                                   <  2:305411/149‑1 (MQ=255)
                                                                                                   ccAGGTCCAGATACTGATTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATg                                                   <  2:310046/149‑1 (MQ=255)
                                                                                                     aGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGca                                                 <  6:213956/149‑1 (MQ=255)
                                                                                                     aGGTCCAGATAATGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGca                                                 <  3:249251/149‑1 (MQ=255)
                                                                                                            gATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAAc                                          >  1:227199/1‑149 (MQ=255)
                                                                                                              tACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGACCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTg                                        >  5:261475/1‑149 (MQ=255)
                                                                                                                cTGTTTGATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGc                                      <  2:172675/149‑1 (MQ=255)
                                                                                                                 tGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGACACCGACTATCATAGGAATGAACAGCAAATCGcc                                     >  5:221958/1‑149 (MQ=255)
                                                                                                                    ttCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCtt                                   <  5:55258/148‑1 (MQ=255)
                                                                                                                     tCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtg                                  <  2:221724/148‑1 (MQ=255)
                                                                                                                          gCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtgagtga                             <  7:353752/148‑1 (MQ=255)
                                                                                                                            ttACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGagcc                          >  6:104134/1‑149 (MQ=255)
                                                                                                                              aCCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGa                                                                 <  6:44351/108‑1 (MQ=255)
                                                                                                                              aCCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGa                                                                 >  5:44351/1‑108 (MQ=255)
                                                                                                                                 cTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGagccagcc                      <  6:180759/149‑1 (MQ=255)
                                                                                                                                    tAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGaa                   >  7:300265/1‑149 (MQ=255)
                                                                                                                                      gCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAcg                 >  5:166837/1‑149 (MQ=255)
                                                                                                                                         aatactgCAGGGGGATCGCCCCTGCCGTGGGATTAATGAGGCGGTTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa              <  4:191626/143‑1 (MQ=255)
                                                                                                                                         aaaaCGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa              >  1:120539/1‑149 (MQ=255)
                                                                                                                                         aaaaCGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa              >  4:151564/1‑149 (MQ=255)
                                                                                                                                         aaaaCGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCACATCAAGACGAGGCAACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGATCCATCCGAAcgcca              >  8:171360/1‑147 (MQ=255)
                                                                                                                                          ccaCGGCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGACCGCGAt             >  2:106837/3‑149 (MQ=255)
                                                                                                                                            aCGGCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTCGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATg            >  1:124611/1‑148 (MQ=255)
                                                                                                                                               gCATGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGTGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAAc        <  7:138561/149‑1 (MQ=255)
                                                                                                                                               gCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa              <  3:351898/143‑1 (MQ=255)
                                                                                                                                               gCAGGGGCGATCCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa              >  4:351898/1‑143 (MQ=255)
                                                                                                                                               gCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAAc        <  6:166837/149‑1 (MQ=255)
                                                                                                                                                 aGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAAcc       >  8:82153/1‑148 (MQ=255)
                                                                                                                                                  gggggATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGa     <  7:82153/149‑1 (MQ=255)
                                                                                                                                                   ggggATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAAcc       <  7:114794/146‑1 (MQ=255)
                                                                                                                                                   ggggATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAAcc       >  8:114794/1‑146 (MQ=255)
                                                                                                                                                   ggggATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAg    <  2:294597/149‑1 (MQ=255)
                                                                                                                                                     ggATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCGTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGtt  <  5:50436/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCTTCTGCAGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGGATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTT  >  NC_002947/5867585‑5867881

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