Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,238,384 G→T 6.2% M510I (ATG→ATT uvrB → uvrABC excinuclease protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,238,3840GT 6.2% 1037.6 / 40.7 308M510I (ATG→ATTuvrBuvrABC excinuclease protein B
Reads supporting (aligned to +/- strand):  ref base G (147/142);  new base T (14/5);  total (161/147)
Fisher's exact test for biased strand distribution p-value = 6.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 308 total aligned reads displayed.

TTACCTGGCCGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATG  >  NC_002947/2238240‑2238531
                                                                                                                                                |                                                                                                                                                    
 tACCTGGCCGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAg                                                                                                                                                <  1:5143511/150‑1 (MQ=255)
 tACCTGGCCGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAg                                                                                                                                                <  2:4687978/150‑1 (MQ=255)
 tACCTGGCCGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAg                                                                                                                                                <  1:4081651/150‑1 (MQ=255)
 tACCTGGCCGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAg                                                                                                                                                <  2:1537132/150‑1 (MQ=255)
     tGGCCGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGgtgt                                                                                                                                            <  1:1403952/150‑1 (MQ=255)
         cGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCt                                                                                                                                       >  2:7232592/1‑150 (MQ=255)
             cACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCtggtg                                                                                                                                   >  2:8741775/1‑150 (MQ=255)
             cACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCtggtg                                                                                                                                   >  1:3703066/1‑150 (MQ=255)
             cACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCtggtg                                                                                                                                   >  1:6778725/1‑150 (MQ=255)
                gacgTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGcc                                                                                                                                >  1:9058071/1‑150 (MQ=255)
                   gTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATc                                                                                                                             >  1:5254235/1‑150 (MQ=255)
                   gTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATc                                                                                                                             >  1:4671561/1‑150 (MQ=255)
                   gTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATc                                                                                                                             >  1:4333228/1‑150 (MQ=255)
                       gcgTACGTTACATGCACTCGGACATCGACACGCTGGAACGCATCGAAATCATCCGCGATCTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGaa                                                                                                                                                <  2:6792856/126‑2 (MQ=255)
                        ggTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGa                                                                                                                        <  1:5397703/150‑1 (MQ=255)
                         gTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGAt                                                                                                                       <  2:5713903/150‑1 (MQ=255)
                                aCCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACa                                                                                                                <  1:1819131/150‑1 (MQ=255)
                                 ccTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACaa                                                                                                               <  2:1855224/150‑1 (MQ=255)
                                 ccTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACaa                                                                                                               <  2:4238380/150‑1 (MQ=255)
                                   tGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCGCCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACaagg                                                                                                             >  1:3958774/1‑150 (MQ=255)
                                   tGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACaagg                                                                                                             <  1:1059971/150‑1 (MQ=255)
                                      aCTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACaaggaag                                                                                                          <  1:8097919/150‑1 (MQ=255)
                                      aCTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACaaggaag                                                                                                          <  1:1219621/150‑1 (MQ=255)
                                       cTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACaaggaagg                                                                                                         <  2:26774/150‑1 (MQ=255)
                                       cTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACaaggaagg                                                                                                         <  2:7450414/150‑1 (MQ=255)
                                        tCGGACATCGACACGCTGGAACGCATCGAAATCATCCGCGATCTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACaaggaaggc                                                                                                        >  1:3456907/1‑149 (MQ=255)
                                        tCGGACATCGACACGCTGGAACGCATCGAAATCATCCGCGATCTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACaaggaaggc                                                                                                        >  1:4251976/1‑149 (MQ=255)
                                         cGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTTGGCACCTTCGTCGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGtt                                                                                                       <  1:7072131/150‑1 (MQ=255)
                                           gACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGtttt                                                                                                     <  1:8300264/150‑1 (MQ=255)
                                            aCATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTc                                                                                                    <  1:6452846/150‑1 (MQ=255)
                                            aCATCGACACCGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTc                                                                                                    <  2:5133325/150‑1 (MQ=255)
                                              aTCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTg                                                                                                  >  2:173512/1‑150 (MQ=255)
                                                 gACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGt                                                                                               >  2:689842/1‑150 (MQ=255)
                                                  acacGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGtt                                                                                              <  1:878389/150‑1 (MQ=255)
                                                  acacGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGtt                                                                                              >  2:8386810/1‑150 (MQ=255)
                                                   cacGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCGCCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTc                                                                                             <  2:3958411/150‑1 (MQ=255)
                                                   cacGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTc                                                                                             <  1:2431261/150‑1 (MQ=255)
                                                    acGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTcc                                                                                            <  1:3897326/150‑1 (MQ=255)
                                                      ggTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATTCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGa                                                                                          <  2:8599345/150‑1 (MQ=255)
                                                      ggTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGa                                                                                          <  1:338856/150‑1 (MQ=255)
                                                      ggTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGa                                                                                          <  2:5901008/150‑1 (MQ=255)
                                                       gTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAg                                                                                         <  1:5037145/150‑1 (MQ=255)
                                                           aGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCt                                                                                     <  2:7963291/150‑1 (MQ=255)
                                                            gCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTc                                                                                    >  2:6526168/1‑150 (MQ=255)
                                                                  ggAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGAt                                                                              <  1:7238048/150‑1 (MQ=255)
                                                                   gagaTTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATc                                                                             <  2:3443815/150‑1 (MQ=255)
                                                                   gagaTTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATc                                                                             >  1:7458405/1‑150 (MQ=255)
                                                                     gaTTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCa                                                                           >  1:8533092/1‑150 (MQ=255)
                                                                       ttATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGa                                                                         <  1:8624507/150‑1 (MQ=255)
                                                                       ttATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGa                                                                         <  2:3144320/150‑1 (MQ=255)
                                                                       ttATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGa                                                                         >  2:1676321/1‑150 (MQ=255)
                                                                          tCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCa                                                                      <  2:148734/150‑1 (MQ=255)
                                                                            cgcgATCTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACAAGGAAGGCTTCCTACGTTCGGAAACTTCGCTGGTCCAGACCAtc                                                                    >  1:7858952/1‑149 (MQ=255)
                                                                            cgcgACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCAtt                                                                    <  1:1676321/150‑1 (MQ=255)
                                                                            cgcgACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCAtt                                                                    <  1:7367412/150‑1 (MQ=255)
                                                                             gcgACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTg                                                                   >  2:5494882/1‑150 (MQ=255)
                                                                             gcgACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTg                                                                   >  2:4465091/1‑150 (MQ=255)
                                                                             gcgACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTg                                                                   <  1:6055733/150‑1 (MQ=255)
                                                                             gcgACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTg                                                                   <  2:5079687/150‑1 (MQ=255)
                                                                              cgACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTgg                                                                  <  2:3301431/150‑1 (MQ=255)
                                                                                atcTGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACAAGGAAGGCTTCCTACGTTCGGAAACTTCGCTGGTCCAGACCATCggtc                                                                <  1:2612732/148‑3 (MQ=255)
                                                                                  cTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCcgc                                                              >  2:1627920/1‑150 (MQ=255)
                                                                                   tGCGCCTTGGCGCTTTCGATGTATTGGTGGGTATCAACCTGCTGCGCGAGGGTCTCGATATTCCGGAATGCGGTTTCGTCGCCATTCTGGATGCCGACAAGGAAGGCTTCCTACGTTCGGAAACTTCGCTGGTCCAGACCATCGGTcgtg                                                             >  1:2583932/1‑148 (MQ=255)
                                                                                   tGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCt                                                                                                   >  2:123957/1‑112 (MQ=255)
                                                                                   tGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCcgcg                                                             >  2:2292323/1‑150 (MQ=255)
                                                                                    gcgcCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCcgcgc                                                            >  2:3538131/1‑150 (MQ=255)
                                                                                        cTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCa                                                        <  2:7102586/150‑1 (MQ=255)
                                                                                            gCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTa                                                    >  1:6855939/1‑150 (MQ=255)
                                                                                            gCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTa                                                    >  1:2180892/1‑150 (MQ=255)
                                                                                            gCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTa                                                    >  2:5073645/1‑150 (MQ=255)
                                                                                              aCCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAAc                                                  >  1:6314104/1‑150 (MQ=255)
                                                                                              aCCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAAc                                                  >  1:1714222/1‑150 (MQ=255)
                                                                                                cTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCt                                                >  1:5924857/1‑150 (MQ=255)
                                                                                                cTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCt                                                >  1:8190310/1‑150 (MQ=255)
                                                                                                 ttCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTc                                               >  1:3320075/1‑150 (MQ=255)
                                                                                                 ttCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTc                                               >  1:8332600/1‑150 (MQ=255)
                                                                                                 ttCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTc                                               >  2:205337/1‑150 (MQ=255)
                                                                                                   cGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCaa                                             <  1:3809238/150‑1 (MQ=255)
                                                                                                        tGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGcc                                        >  1:2325966/1‑150 (MQ=255)
                                                                                                          cTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCgg                                      <  2:4156751/150‑1 (MQ=255)
                                                                                                          cTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCgg                                      <  1:2769916/150‑1 (MQ=255)
                                                                                                          cTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCgg                                      >  2:2769916/1‑150 (MQ=255)
                                                                                                           tggtggGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCggg                                     <  1:5680381/150‑1 (MQ=255)
                                                                                                              tggGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCa                                  >  1:700160/1‑150 (MQ=255)
                                                                                                                ggCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCAtt                                <  2:8044237/150‑1 (MQ=255)
                                                                                                                 gCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTc                               <  2:4093434/150‑1 (MQ=255)
                                                                                                                  cATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCt                              <  2:222637/150‑1 (MQ=255)
                                                                                                                   aTCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTg                             <  2:5591545/150‑1 (MQ=255)
                                                                                                                   aTCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTg                             >  1:7664902/1‑150 (MQ=255)
                                                                                                                   aTCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTg                             >  2:4570577/1‑150 (MQ=255)
                                                                                                                   aTCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCCGGCCAttgag                             <  2:7711581/150‑4 (MQ=255)
                                                                                                                       aCCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATAGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACg                         >  2:4843190/1‑150 (MQ=255)
                                                                                                                           gctgcGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTGTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGa                     <  2:3871140/150‑1 (MQ=255)
                                                                                                                                     ggCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACa           <  1:4392428/150‑1 (MQ=255)
                                                                                                                                     ggCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACa           <  1:5585997/150‑1 (MQ=255)
                                                                                                                                        cTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTATGCCGATAACGTAACAGGc        >  1:5782075/1‑150 (MQ=255)
                                                                                                                                        cTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGc        >  2:6384080/1‑150 (MQ=255)
                                                                                                                                         tGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCt       >  2:6067123/1‑150 (MQ=255)
                                                                                                                                          ggACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTc      <  1:2990086/150‑1 (MQ=255)
                                                                                                                                              aTGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATg  >  1:1136969/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTACCTGGCCGACCACGACGTGCGGGTACGCTACCTGCACTCGGACATCGACACGGTGGAGCGGGTGGAGATTATCCGCGACCTGCGCCTGGGCACCTTCGACGTGCTGGTGGGCATCAACCTGCTGCGCGAGGGCCTGGACATGCCCGAGGTGTCGCTGGTGGCCATCCTCGATGCGGACAAGGAAGGTTTTCTGCGTTCCGAGCGCTCCCTGATCCAGACCATTGGCCGCGCGGCACGTAACCTCAATGGCCGGGCCATTCTGTACGCCGATAACGTAACAGGCTCCATG  >  NC_002947/2238240‑2238531

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