Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,908,805 T→A 7.3% intergenic (+56/+51) PP_5176 → / ← spuH hypothetical protein/spermidine/putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,908,8050TA 7.3% 644.3 / 5.9 239intergenic (+56/+51)PP_5176/spuHhypothetical protein/spermidine/putrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (114/85);  new base A (15/6);  total (141/98)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

TCATGTAGCCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACG  >  NC_002947/5908656‑5908954
                                                                                                                                                      |                                                                                                                                                     
   tGTAGCCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGTGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGaa                                                                                                                                                    <  2:2914738/150‑1 (MQ=255)
   tGTAGCCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGaa                                                                                                                                                    <  1:3906510/150‑1 (MQ=255)
   tGTAGCCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGaa                                                                                                                                                    <  2:6509/150‑1 (MQ=255)
     gcaaggcgattcagcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAgg                                                                                                                                                  >  2:5509367‑M1/116‑150 (MQ=255)
     gcaaggcgattcagcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAgg                                                                                                                                                  >  2:3994466‑M1/116‑150 (MQ=255)
        cccGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCGGAAGGcgg                                                                                                                                               >  2:701312/1‑148 (MQ=255)
           cgattcagcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGa                                                                                                                                            >  2:52755‑M1/110‑150 (MQ=255)
            gattcagcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAg                                                                                                                                           >  1:5534641‑M1/109‑150 (MQ=255)
            gTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAg                                                                                                                                           <  2:2835692/150‑1 (MQ=255)
             tACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAgg                                                                                                                                          <  1:3642370/150‑1 (MQ=255)
              ttcagcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGt                                                                                                                                         >  1:5088859‑M1/107‑150 (MQ=255)
               tcagcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTc                                                                                                                                        >  2:1393071‑M1/106‑150 (MQ=255)
                cagcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTcc                                                                                                                                       >  1:2310984‑M1/105‑150 (MQ=255)
                aTTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACGATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTcc                                                                                                                                       >  2:805787‑M2/1‑148 (MQ=255)
                 agcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTccc                                                                                                                                      >  1:4064445‑M1/104‑150 (MQ=255)
                    aggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgt                                                                                                                                   >  1:2316447‑M1/101‑150 (MQ=255)
                    aggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGACCCCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                  >  1:2081651‑M1/101‑150 (MQ=255)
                     tgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTccctgtg                                                                                                                                  >  2:626572‑M2/1‑143 (MQ=255)
                     tgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTccctgtg                                                                                                                                  >  1:1251700‑M2/1‑143 (MQ=255)
                     ggcgttcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                  >  1:3762443‑M1/100‑150 (MQ=255)
                     ggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                  >  2:4759426‑M1/100‑150 (MQ=255)
                     ggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                  >  2:3448333‑M1/100‑150 (MQ=255)
                     aggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggcccttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                  >  2:2171469‑M1/100‑150 (MQ=255)
                      gcgatcgcagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                 >  2:3824388‑M1/99‑150 (MQ=255)
                      gcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                 >  2:1987169‑M1/99‑150 (MQ=255)
                      gCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAACGCTGTGAGGTccctgtgg                                                                                                                                 <  2:4901124‑M2/150‑9 (MQ=255)
                        gatcgacgagggtgctgcagccaacgcttcgcagccgcaggtcacgcgcccggcgcaggtggctgcgtcggcctgacctggccccttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGa                                                                                                                               >  1:2773042‑M1/97‑150 (MQ=255)
                        cAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTccctgt                                                                                                                                   <  2:4153933‑M2/146‑7 (MQ=255)
                        cAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTccctgt                                                                                                                                   >  1:4153933‑M2/1‑140 (MQ=255)
                         aaTCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATTCCCGCTCCCACAGGGCTCACCACAACCCAGAAGGCTGTGAGGTccctgggggcg                                                                                                                              <  2:326552‑M2/150‑12 (MQ=255)
                         aaTCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCCCGCCCACAGGGCACACCACAACACAGAAGGCCGTGTGGTgcctgtgggag                                                                                                                              <  1:531968‑M2/150‑12 (MQ=255)
                           tCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCAAGCCCGCTCCCACCCGGCGCCACCCCAAACAGAAGGCTGTGAGGTcccggggggggcg                                                                                                                            <  1:1419485‑M2/150‑14 (MQ=255)
                             ccTTGAGGGCTGGATTGGTCATAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCCGGCCTGCCCCCACCCACACGGCCAAACACAAACCAACCGGCTGGGCGGGGCCGGTGGGAGCggg                                                                                                                          >  1:1404237/1‑150 (MQ=255)
                             ccTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATTCCCCCTACCACAGGGCTCACCACAACCCTGAAGGCTGTGGGGTcccggtgggagcggg                                                                                                                          <  1:4452460‑M2/150‑16 (MQ=255)
                               ttGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACACGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCg                                                                                                                        >  2:4266871/1‑150 (MQ=255)
                               ttGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCAAGCCCGCGCCCACAGGGCTCACCACAACCCAGAAGGCTGTGAGGTCCCGGTGGGAGCGGGCg                                                                                                                        >  1:2749684/1‑150 (MQ=255)
                                tGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGt                                                                                                                       <  2:608497/150‑1 (MQ=255)
                                     gCTGGTTTGGTCAGAAGTATCAGGGACTATCCACATTGGGTTGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGCGGTCCCTGTGGGAGCGGGCGTGCCcg                                                                                                                  <  1:715644/150‑1 (MQ=255)
                                        ggATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGa                                                                                                               <  2:4175151/150‑1 (MQ=255)
                                              gTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCGCCCACAGGGCTCACCACAACCCGGAAGGCAGTGAGGTCCCGGTGGGGGCCGGCGTGCCCGCGCAAGGGc                                                                                                         >  1:2273272/1‑150 (MQ=255)
                                               tCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGcc                                                                                                        >  2:4102393/1‑150 (MQ=255)
                                               tCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCGGAAGGCGGGGAGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGcc                                                                                                        >  2:2516130/1‑150 (MQ=255)
                                               tCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCAGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCCAAAGGGcc                                                                                                        >  2:4039708/1‑150 (MQ=255)
                                                 ataaGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGATTCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAg                                                                                                      <  1:3676840/148‑1 (MQ=255)
                                                 aGAAGTATCAGGGCCTATCCACCTTGGGTTGCGGTATGGATTTGAATCGGTGCTGTAGCCTTCGCGGGGAGGCCCGCTCCCACCGGGCTCACCCCATCCCTGACGGCTTTGAGTTCCCTGTGTGAGCGGGCGTGCCCGCGAAAGGGCCAg                                                                                                      <  1:847431/150‑1 (MQ=255)
                                                      tATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCGCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAg                                                                                                 >  2:1761414/1‑150 (MQ=255)
                                                           gggACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCGGTGAGGTCCCTGTGGGAGCGGGCGTGCCCCCGAAAGGGCCAGGGCAGGCCGa                                                                                            >  2:1856636/1‑150 (MQ=255)
                                                                    cacaTTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCCTGCCCGCTCCCACCCGGCTCACCCCAACCCGGAAGGCGGTGAGGTCCCGGGGGGGGCGGGCGTGCCCCCGGAAGGGCCAGGGGAGGCCCGCGCCGCCAc                                                                                   >  1:2639502/1‑150 (MQ=255)
                                                                    cacaTTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCCGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCAc                                                                                   <  2:1224819/150‑1 (MQ=255)
                                                                       aTTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCCCGCCCCCAGGGCTCACCACAACCCAGAAGGCAGTGAGGTACCAGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCAccag                                                                                >  1:4497209/1‑148 (MQ=255)
                                                                            gATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCCCGCCCACACGGCTCACCCCAACCCGGAATGCTGTGTGGTCCCGGGGGGAGCGGGCGGGGCCGCGCCAGGGCCCGGGCAGGCCGCCGCCGCCACCCGCGCCg                                                                           >  1:559196/1‑150 (MQ=255)
                                                                               gCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGgc                                                                        <  2:1261932/150‑1 (MQ=255)
                                                                                  gTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCt                                                                     >  1:4668559/1‑150 (MQ=255)
                                                                                   gagggCTTTGAATCGGTGCTGGCGCCTTCGCGGGCATGCCCGCGCCCACCGGGCCCACCCCCCCCCGGACGGCGGTGAGGGCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCtt                                                                    <  2:3056098/147‑1 (MQ=255)
                                                                                   gaTGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCtt                                                                    <  2:4450453/149‑1 (MQ=255)
                                                                                                cGGTGGTGGCGCCTGCGCGGGCAGGCCCGCCCCCACCGGGCTCACCAAAACCCTGAAGGGTGTGCGGTCCCTGTGGGCGGGGGGGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGc                                                       <  1:1459702/150‑1 (MQ=255)
                                                                                                     cgggAGCCTTCGCGGGCAGGCCCGCGCCCACCGGGCTCACCACACCCCTGAAGGCTGGGAGGTCCCTGGGGGCGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGc                                                  <  1:69242/148‑1 (MQ=255)
                                                                                                       ggAGCCTTCGCGGGCAAGCCCGCGCCCACAGGGCTCACCACAACCCAGTAGGCGGTGCGGTCCCGGGGGGCGCGGGCGTGCCCGCGGCAGGGCCCCGTCAGGGCGACGCAGCCCCCAGCGCCCGGCGCGTGGCCTGCGGCTGCGAACCgg                                                >  2:669785/1‑149 (MQ=255)
                                                                                                          gCCTGCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTgg                                             <  2:1191041/150‑1 (MQ=255)
                                                                                                             ttCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTg                                          >  1:5405344/1‑150 (MQ=255)
                                                                                                             ttCGCGGGCATGCCCGCTACCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGGTGGCTg                                          >  2:48347/1‑150 (MQ=255)
                                                                                                               cgcgggcagGCCCGCCCCCACCGGGCTCACCACAACCCTTAAGGCTGTGAGGTCCCTGGGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgca                                        >  1:1914984‑M1/10‑150 (MQ=255)
                                                                                                               cgcgGGCATGCCCGCTCCCACAGGGCTCACCACACCCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgca                                        <  2:4663211/150‑1 (MQ=255)
                                                                                                               cgcgGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTCACGGCGGTGTGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgca                                        <  1:4251023/150‑1 (MQ=255)
                                                                                                                      aTGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCct                                 >  1:3584895/1‑150 (MQ=255)
                                                                                                                         CCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCtt                              <  2:557889‑M1/150‑1 (MQ=255)
                                                                                                                          CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGAc                                                                 >  1:385786‑M1/1‑114 (MQ=255)
                                                                                                                          CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGACAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                             <  1:889786‑M1/150‑1 (MQ=255)
                                                                                                                              cCCCACCGGGCTCACCACAACCGTGACGGCTGTGTGGTCCCTGTGGGAGCGGGGGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATc                         >  1:2516130‑M1/2‑150 (MQ=255)
                                                                                                                              TCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATc                         >  2:3326344‑M1/1‑150 (MQ=255)
                                                                                                                                ccACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTccctgtgggagcgggcgtgcccgcgaaggctccagcaccgattcaaagccataccgcatcccaatgtggatagtccctgatacttctgaccaatccagccctcaagggattggc                       <  1:1509596‑M2/150‑115 (MQ=255)
                                                                                                                                ccACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTccctgtgggagcgggcgtgcccgcgaaagctccagcaccgattcaaagccataccgcatcccaatgtggatagtccctgatacttctgaccaatccagccctcaagggattggc                       <  2:2015201‑M2/150‑115 (MQ=255)
                                                                                                                                ccACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTccctgtgggagcgggcatgcccgcgaaggctccagcaccgattcaaagccataccgcatcccaatgtggatagtccctgatacttctgaccaatccagccctcaagggattggc                       <  2:4815492‑M2/150‑115 (MQ=255)
                                                                                                                                CCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGc                       <  2:1316960‑M1/150‑1 (MQ=255)
                                                                                                                                CCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGc                       <  1:301353‑M1/150‑1 (MQ=255)
                                                                                                                                CCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGc                       <  1:3226998‑M1/150‑1 (MQ=255)
                                                                                                                                 CACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGcc                      <  2:3664767‑M1/150‑1 (MQ=255)
                                                                                                                                  ACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCct                     <  1:4710602‑M1/150‑1 (MQ=255)
                                                                                                                                     GGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgct                  >  1:3840239‑M1/1‑150 (MQ=255)
                                                                                                                                            accacAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCg            >  2:4987215/1‑150 (MQ=255)
                                                                                                                                            accacAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCg            >  1:1036242/1‑150 (MQ=255)
                                                                                                                                             ccacAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGc           >  2:2468945/1‑150 (MQ=255)
                                                                                                                                             ccacAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGc           >  2:1638806/1‑150 (MQ=255)
                                                                                                                                              cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGTGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc          >  1:1726648/1‑150 (MQ=255)
                                                                                                                                              cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc          >  2:1890457/1‑150 (MQ=255)
                                                                                                                                              cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc          >  1:3288440/1‑150 (MQ=255)
                                                                                                                                               acaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCAGCTGAAAAGCCt         >  1:45620/1‑150 (MQ=255)
                                                                                                                                                caACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCtt        <  2:4853900/150‑1 (MQ=255)
                                                                                                                                                caACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCtt        <  1:3394629/150‑1 (MQ=255)
                                                                                                                                                 aaCCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGc                                                                                                         >  2:861329/1‑52 (MQ=255)
                                                                                                                                                  aCCCTTAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCTATTGACCTGCGGCTGCTAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGc      <  1:4992969/150‑1 (MQ=255)
                                                                                                                                                   cccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcg     >  2:3515023/1‑150 (MQ=255)
                                                                                                                                                    ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGc                       <  1:165512/131‑1 (MQ=255)
                                                                                                                                                    ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcga    >  2:1462961/1‑150 (MQ=255)
                                                                                                                                                    ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcga    >  2:1548876/1‑150 (MQ=255)
                                                                                                                                                     cTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgac   >  2:5146914/1‑150 (MQ=255)
                                                                                                                                                      tGAAGGCTGTGAGGTccctgtgggagcgggcgtgcccgcgaaggctccagcaccgattcaaagccataccgcatcccaatgtggatagtccctgatacttctgaccaatccagccctcaagggattggcaacaatgtaccgggctacatg  <  1:5445228‑M2/150‑136 (MQ=255)
                                                                                                                                                      tGAAGGCTGTGAGGTccctgtgggagcgggcatgcccgcgaaggctccagcaccgattcaaagccataccgcatcccaatgtggatagtccctgatacttctgaccaatccagccctcaagggattggcaacaatgtaccgggctacatg  <  2:1163857‑M2/150‑136 (MQ=255)
                                                                                                                                                      tGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacg  >  1:1003786/1‑150 (MQ=255)
                                                                                                                                                      tGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGCCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacg  >  2:2314375/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCATGTAGCCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACG  >  NC_002947/5908656‑5908954

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: