Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,381,112 G→C 5.2% A116P (GCT→CCT) 
P65A (CCT→GCT) 
PP_5639 →
PP_3855 ←
hypothetical protein
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,381,1120GC 5.2% 517.3 / 13.4 154A116P (GCT→CCT) 
P65A (CCT→GCT) 
PP_5639
PP_3855
hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (91/55);  new base C (4/4);  total (95/59)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 159 total aligned reads displayed.

ATTCAGACCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTTTT  >  NC_002947/4380963‑4381261
                                                                                                                                                     |                                                                                                                                                     
aTTCAGACCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAgg                                                                                                                                                       >  1:1527559/1‑150 (MQ=255)
 ttCAGACCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGc                                                                                                                                                      >  1:1447520/1‑150 (MQ=255)
          aTACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGa                                                                                                                                             <  1:4350467/150‑1 (MQ=255)
           tACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAg                                                                                                                                            >  2:2839154/1‑150 (MQ=255)
           tACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAg                                                                                                                                            >  2:3231213/1‑150 (MQ=255)
              gagaTTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCACTATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCt                                                                                                                                         >  2:1208950/1‑150 (MQ=255)
               agaTTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCtt                                                                                                                                        >  2:195567/1‑150 (MQ=255)
               agaTTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCtt                                                                                                                                        >  2:195675/1‑150 (MQ=255)
                 aTTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGAAGGCTCCGGCAGAGCCTGCg                                                                                                                                      <  1:3849566/150‑1 (MQ=255)
                 aTTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCAGAGCCTGCGCCGGAGCCtgag                                                                                                                                      <  1:4612913/150‑4 (MQ=255)
                  ttCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGc                                                                                                                                     <  1:5301286/150‑1 (MQ=255)
                   tCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATCTGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCa                                                                                                                                    >  1:5376525/1‑150 (MQ=255)
                      gATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATg                                                                                                                                 >  1:2107682/1‑150 (MQ=255)
                         cGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGAGGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATc                                                                                                                              <  2:1909697/150‑1 (MQ=255)
                            cATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGac                                                                                                                           <  1:2647193/150‑1 (MQ=255)
                                        tGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTg                                                                                                               >  2:2484532/1‑150 (MQ=255)
                                         ggCGTTGCTGACACAATTGTTCGGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGt                                                                                                              <  1:3647525/150‑1 (MQ=255)
                                           cGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTaa                                                                                                            <  1:4701335/150‑1 (MQ=255)
                                               gCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCt                                                                                                        <  2:5349387/150‑1 (MQ=255)
                                                cTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCtt                                                                                                       <  1:2508210/150‑1 (MQ=255)
                                                 tGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTc                                                                                                      <  1:1647996/150‑1 (MQ=255)
                                                  gACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCa                                                                                                     <  2:3468762/150‑1 (MQ=255)
                                                  gACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCa                                                                                                     >  1:3585775/1‑150 (MQ=255)
                                                     acaATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCt                                                                                                                >  2:1879094/1‑136 (MQ=255)
                                                       aaTTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTAcc                                                                                                >  2:2734868/1‑150 (MQ=255)
                                                       aaTTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTAcc                                                                                                >  2:2716100/1‑150 (MQ=255)
                                                       aaTTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTAcc                                                                                                >  1:1154886/1‑150 (MQ=255)
                                                         ttGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTg                                                                                              >  1:3024233/1‑150 (MQ=255)
                                                                 cccTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTa                                                                                      <  2:3378741/150‑1 (MQ=255)
                                                                  ccTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGAAGGCTCCGGCAGAGc                                                                                                                                           <  1:4813796/96‑1 (MQ=255)
                                                                   cTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAAt                                                                                    <  2:2702674/150‑1 (MQ=255)
                                                                     acacCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCaa                                                                                                                                   >  1:2652938/1‑101 (MQ=255)
                                                                     acacCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCaa                                                                                                                                   <  2:2652938/101‑1 (MQ=255)
                                                                     acacCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCaa                                                                                                                                   >  1:3796806/1‑101 (MQ=255)
                                                                     acacCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCtt                                                                                                       >  2:216703/1‑129 (MQ=255)
                                                                     acacCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCtt                                                                                                       <  1:216703/129‑1 (MQ=255)
                                                                     acacCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATAc                                                                                  >  1:371370/1‑150 (MQ=255)
                                                                     acacAAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTGAATAc                                                                                  >  1:3206617/1‑150 (MQ=255)
                                                                         cAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTgc                                                                              >  1:4030162/1‑150 (MQ=255)
                                                                         cAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTgc                                                                              >  1:4030463/1‑150 (MQ=255)
                                                                         cAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTgc                                                                              >  1:809657/1‑150 (MQ=255)
                                                                         cAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTgc                                                                              <  1:3436403/150‑1 (MQ=255)
                                                                         cAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTgc                                                                              >  1:4775543/1‑150 (MQ=255)
                                                                          aGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTgcg                                                                             >  1:4239509/1‑150 (MQ=255)
                                                                               tCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTc                                                                        >  1:1934303/1‑150 (MQ=255)
                                                                                 aGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                      >  2:1438198/1‑150 (MQ=255)
                                                                                 aGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                      >  2:4185120/1‑150 (MQ=255)
                                                                                 aGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                      >  1:1881279/1‑150 (MQ=255)
                                                                                 aGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                      >  1:4621387/1‑150 (MQ=255)
                                                                                 aGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                      >  2:2991121/1‑150 (MQ=255)
                                                                                 aGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                      >  1:3548397/1‑150 (MQ=255)
                                                                                  gCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCggg                                                                     >  1:101225/1‑150 (MQ=255)
                                                                                  gCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCggg                                                                     >  1:1445199/1‑150 (MQ=255)
                                                                                  gCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCggg                                                                     >  2:1273259/1‑150 (MQ=255)
                                                                                   cAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGa                                                                    <  2:3958745/150‑1 (MQ=255)
                                                                                     aTACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAAc                                                                  >  1:2224963/1‑150 (MQ=255)
                                                                                      tACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACt                                                                 <  2:521379/150‑1 (MQ=255)
                                                                                      tACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACt                                                                 >  2:35151/1‑150 (MQ=255)
                                                                                         aaGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGAc                                                              >  1:2148479/1‑150 (MQ=255)
                                                                                          aGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACt                                                             >  2:4625646/1‑150 (MQ=255)
                                                                                           gCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTc                                                            <  2:1509658/150‑1 (MQ=255)
                                                                                           gCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGGCGGGAACTGACTc                                                            <  2:3862766/150‑1 (MQ=255)
                                                                                           gCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCATTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGAATc                                                            >  1:1874298/1‑150 (MQ=255)
                                                                                            cTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCa                                                           <  2:3791628/150‑1 (MQ=255)
                                                                                              aGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAgg                                                         <  1:3525232/150‑1 (MQ=255)
                                                                                                         ttaTACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                      >  1:1689682/3‑126 (MQ=255)
                                                                                                         ttaTACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                      <  2:1689682/124‑1 (MQ=255)
                                                                                                         tGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGc                                              <  2:2928658/150‑1 (MQ=255)
                                                                                                          gATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCa                                             >  2:4894380/1‑150 (MQ=255)
                                                                                                             aCCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACgg                                          <  2:1234679/150‑1 (MQ=255)
                                                                                                               cAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGt                                        >  1:4119021/1‑150 (MQ=255)
                                                                                                               cAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGt                                        >  2:5335288/1‑150 (MQ=255)
                                                                                                               cAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGt                                        >  2:5013579/1‑150 (MQ=255)
                                                                                                                aaGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTa                                       >  2:5154133/1‑150 (MQ=255)
                                                                                                                aaGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTa                                       >  1:1168285/1‑150 (MQ=255)
                                                                                                                    ttCAGCCCGCAGATGCCTCAGGCTCCGGCGTAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAg                                   >  2:3746521/1‑150 (MQ=255)
                                                                                                                          ccGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCcct                             >  1:629268/1‑150 (MQ=255)
                                                                                                                           cGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCcctc                            >  2:3356359/1‑150 (MQ=255)
                                                                                                                                aTGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCtgtg                       <  1:1043784/150‑1 (MQ=255)
                                                                                                                                  gCCTCAGGCTCCGGCGCGGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTgga                     >  1:4306771/1‑150 (MQ=255)
                                                                                                                                    cgaaGGCTCCGGCAGAGCCTGCGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCGGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagc                   >  1:271211/4‑150 (MQ=255)
                                                                                                                                    cTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagc                   <  1:1791949/150‑1 (MQ=255)
                                                                                                                                    cTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagc                   <  1:1919545/150‑1 (MQ=255)
                                                                                                                                     tCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGt                                                                                        <  1:1973601/80‑1 (MQ=255)
                                                                                                                                     tCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGt                                                                                        >  2:1973601/1‑80 (MQ=255)
                                                                                                                                      cAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcgg                 >  2:3458991/1‑150 (MQ=255)
                                                                                                                                       aGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcgga                >  2:5404806/1‑150 (MQ=255)
                                                                                                                                        ggCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGACCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggag               <  2:2224963/150‑1 (MQ=255)
                                                                                                                                         gCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggagc              <  1:5015595/150‑1 (MQ=255)
                                                                                                                                         gCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggagc              >  1:4904603/1‑150 (MQ=255)
                                                                                                                                         gCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggagc              <  1:3231213/150‑1 (MQ=255)
                                                                                                                                         gCTCCGGCAGAGCCTGCGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggagc              >  1:3089123/1‑150 (MQ=255)
                                                                                                                                              ggCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCtg         <  2:1527559/150‑1 (MQ=255)
                                                                                                                                                cgcAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCtgtt       <  2:1447520/150‑1 (MQ=255)
                                                                                                                                                cgcAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCtgtt       <  1:195675/150‑1 (MQ=255)
                                                                                                                                                cgaaGGCTCCGGCAGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCtgtt       >  1:5078080/4‑150 (MQ=255)
                                                                                                                                                  cAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCtgttgt     >  1:4505582/1‑150 (MQ=255)
                                                                                                                                                    ggCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGttt   >  2:4831217/1‑150 (MQ=255)
                                                                                                                                                     gCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGtttt  >  2:2059765/1‑150 (MQ=255)
                                                                                                                                                     gCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGtgtt  >  2:4797138/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTCAGACCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTTTT  >  NC_002947/4380963‑4381261

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: