Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,899,748 A→G 23.7% intergenic (‑170/+64) PP_5172 ← / ← PP_5173 hypothetical protein/RND family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,899,7480AG23.7% 234.5 / 69.7 124intergenic (‑170/+64)PP_5172/PP_5173hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base A (39/52);  new base G (6/23);  total (45/79)
Fisher's exact test for biased strand distribution p-value = 4.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

Only 100 of 125 total aligned reads displayed.

GGCGCCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGA  >  NC_002947/5899600‑5899897
                                                                                                                                                    |                                                                                                                                                      
    ccGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggtg                                                                                                                                                   <  1:1222910/150‑1 (MQ=255)
    ccGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggtg                                                                                                                                                   <  1:5145053/150‑1 (MQ=255)
    ccGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggtg                                                                                                                                                   <  2:5038032/150‑1 (MQ=255)
     cGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATAAggggg                                                                                                                                                  >  1:3880113/1‑150 (MQ=255)
      gggCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggtggt                                                                                                                                                 <  1:3018662/150‑1 (MQ=255)
       ggCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTc                                                                                                                                                <  2:4955732/150‑1 (MQ=255)
       ggCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTc                                                                                                                                                <  2:1140282/150‑1 (MQ=255)
       ggCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTc                                                                                                                                                <  2:650355/150‑1 (MQ=255)
          agagGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAAAAAGGTTGTCGAt                                                                                                                                             >  2:4295663/1‑150 (MQ=255)
           gagGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCgatggtgcatc                                                                                                                                            <  2:4409255/150‑6 (MQ=255)
           gagGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCgatggtgcatc                                                                                                                                            <  1:4709022/150‑6 (MQ=255)
              gTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCAc                                                                                                                                         <  1:2973443/150‑1 (MQ=255)
                  gtaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                     <  2:3558344/148‑1 (MQ=255)
                  gatggggTTTTGCGCCTGGCTTGAGGGCATTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                     <  1:3447153/147‑1 (MQ=255)
                  gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                     <  1:2114210/150‑1 (MQ=255)
                  gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                     <  2:2070500/150‑1 (MQ=255)
                  gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                     <  2:4322570/150‑1 (MQ=255)
                      gggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTTCCCGCGAATACGGTGGTCGATGCACCATCGCAt                                                                                                                                >  1:4439190/1‑150 (MQ=255)
                          tttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                            <  1:2253970/150‑1 (MQ=37)
                           ttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgcg                                                                                                                           <  1:4734010/150‑1 (MQ=37)
                           ttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgcg                                                                                                                           <  1:1402708/150‑1 (MQ=25)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                          <  2:2673079/150‑1 (MQ=25)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                          <  2:4828466/150‑1 (MQ=25)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                          <  1:1489676/150‑1 (MQ=25)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                          <  2:283479/150‑1 (MQ=25)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                          <  1:175669/150‑1 (MQ=25)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                          <  1:5180645/150‑1 (MQ=25)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                          <  2:1970726/150‑1 (MQ=25)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                          <  2:2956559/150‑1 (MQ=255)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                          <  2:859915/150‑1 (MQ=255)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                          <  2:639016/150‑1 (MQ=37)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                          <  1:3921393/150‑1 (MQ=25)
                            tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                          <  2:3510557/150‑1 (MQ=25)
                            ggcgcCTGGCTTGAGGGCAATAGTAGTGAGGCGCGATTATTCTTAAAAATACTGTTGAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                          <  1:231357/149‑1 (MQ=21)
                             gcgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCggg                                                                                                                         <  2:3830533/150‑1 (MQ=255)
                                  tGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAAACCCCAAAACCGTGGTCGATTCACCATCGCCATCGCGGGTAAAc                                                                                                                    >  1:3199664/1‑150 (MQ=255)
                                  gggCTTGAGGGCAAGAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTCAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAc                                                                                                                    <  2:3139535/149‑1 (MQ=255)
                                    gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAAACCCCAAAACCGTGGTCGATTCACCCACGCATTCGCCGGGAAAacc                                                                                                                  >  1:4634242/1‑150 (MQ=255)
                                      ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGc                                                                                                                >  2:1833496/1‑150 (MQ=255)
                                      ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGc                                                                                                                <  1:4169722/150‑1 (MQ=255)
                                      ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGc                                                                                                                >  1:1263963/1‑150 (MQ=255)
                                        gAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTACATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTc                                                                                                              <  2:3480644/150‑1 (MQ=255)
                                        gAGGGCAATAGGAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAATAATTTTTACATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTc                                                                                                              <  2:3480269/150‑1 (MQ=255)
                                           ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGCATACGGTGGTCGTTGCACCCTCGCATTCGCGGGTAAAACCGCTCCTa                                                                                                           >  1:3456539/1‑150 (MQ=255)
                                           ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCAGAGAATAAGGTGGTCGATGCACCAAAGCATTAACGGGTAAAACCGCtacaa                                                                                                           >  1:3464301/1‑146 (MQ=255)
                                              aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCCCGAAAAAGGTGGTCGATGCACCAACGCAATAGCGGGGAAAACCGCaccaacca                                                                                                        >  1:1587638/1‑150 (MQ=255)
                                              aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGAAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCTCATTCGCGGGTAAACCCGCTCCTAccc                                                                                                        >  2:5225390/1‑149 (MQ=255)
                                               aTAGTAGTGAGGCGGTATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGAGGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAg                                                                                                       <  2:4165101/150‑1 (MQ=255)
                                                tagtagTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATTCACCAACGCATTCGCGGGTAAACCCGCTCCAACCAgg                                                                                                      >  2:1960805/1‑150 (MQ=255)
                                                 agtagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGACCACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAggg                                                                                                     <  2:834567/150‑1 (MQ=255)
                                                  ggtgtGAGGCGGGATTATTCTTAAAAATACTTTTTAAGAATTTTTATATAACATTTTAAGAGACCACGCCTCGGGCAGGAGCGGGTTTCCCCGCGAATACGGTGGTCTATGCACCATCTCCGGCGCGGGTAAACCCGCTCCTACCAgggg                                                                                                    <  1:5225390/147‑1 (MQ=255)
                                                        aGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCAc                                                                                                                                         >  2:4697032/1‑108 (MQ=255)
                                                        aGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCAc                                                                                                                                         <  1:4697032/108‑1 (MQ=255)
                                                          gCGGGATATTTATTAAAAATATTTTTTAATATTTTTTATATAACCATTTAAGATAACACGCCTCTGGTTGGAGCGGTTTTACCCGCGATTACGGTGGTCGATGCACCATCGCTTTCGGGGTTAAACCCGCTCCTACCAGGGGCCATATgg                                                                                            <  1:865191/150‑1 (MQ=255)
                                                               attattCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAcc                                                                                                                   >  2:3525328/1‑122 (MQ=255)
                                                               attattCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAcc                                                                                                                   <  1:3525328/122‑1 (MQ=255)
                                                                 tattatTAAAAATACTGTTTAATAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCa                                                                                     <  2:1897831/150‑1 (MQ=255)
                                                                    tCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGc                                                                                  <  2:3669159/150‑1 (MQ=255)
                                                                      ttAAAAATACTTTTTTAGAATTTTTATATAACCTTTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGGGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCgt                                                                                <  2:4338316/150‑1 (MQ=255)
                                                                      ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATAAGGGCCTTCAGGCgt                                                                                >  2:2136921/1‑150 (MQ=255)
                                                                      ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGGTTAACCCCCAAAAAAGTGGTGGATTCAACAAAGCAATACCCGGGAAAAACCCTACAACAAGCGGCCATAAGGGGCATCAGGCgt                                                                                >  1:2801450/1‑150 (MQ=255)
                                                                       tAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGc                                                                                                                                           <  1:1446671/91‑1 (MQ=255)
                                                                       tAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGc                                                                                                                                           >  2:1446671/1‑91 (MQ=255)
                                                                            aTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGGAAACCCGCTCCTAACAGGGGCCAAATGGGCCTTCAGGCGTGTTTCt                                                                          >  1:5327964/1‑150 (MQ=255)
                                                                               cTGTTTAATATTTTTTATATACCCTTTTAATCCAACACGCCTGTGGTAGGCGGGGGTTACCCCGCGAATAGGGTGGTCGATGCACCATCGCCTGCGGGGGAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCg                                                                       <  1:745576/150‑1 (MQ=255)
                                                                                    tAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcg                                                                  <  2:3151912/150‑1 (MQ=255)
                                                                                              tatataACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGt                                                        <  2:3696585/150‑1 (MQ=255)
                                                                                                  aaacatattAACACAACACGCCTGTGTTCGGGGGGGTTTTCCCGGCTAATAGGGTGGTCGAGCCACCAGCGCAGGCGGGGGAAACCCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCa                                                    <  1:1960805/143‑1 (MQ=255)
                                                                                                   aaCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                   <  2:4803183/150‑1 (MQ=255)
                                                                                                       aTTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCa                                               <  2:4469266/150‑1 (MQ=255)
                                                                                                           aaGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                          >  1:262947/1‑71 (MQ=255)
                                                                                                            agagAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGGTGCACCAACGCATTCGCGGGTAAACCCGCGCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTGTGCGGCTCGCGGATGGTGTACCAGGCCACGTAc                                          >  2:1486453/1‑150 (MQ=38)
                                                                                                             gagaACCCGCCTGTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACa                                         <  2:3379332/150‑1 (MQ=255)
                                                                                                                     gCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGCTGCCCCAAAGCAATCCCCGGTAAAACCCCTCCGACCACGGGCCCTATGGGCCGTCCGGCGGGTGTCTGCGGCGCGCGGCGCGTGTACCAGGCCCCCGGCAGCGCCCGc                                 >  2:11582/1‑150 (MQ=255)
                                                                                                                           ggTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGAGGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAc                           <  2:518491/150‑1 (MQ=255)
                                                                                                                             tAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                         <  2:1873944/150‑1 (MQ=255)
                                                                                                                                gAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATCTGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcagcag                      <  2:825009/150‑1 (MQ=255)
                                                                                                                                gAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcagcag                      <  2:824800/150‑1 (MQ=255)
                                                                                                                                 aGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                     <  2:2846912/150‑1 (MQ=255)
                                                                                                                                  gCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCg                    >  1:1223460/1‑150 (MQ=255)
                                                                                                                                            ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTa                                                       <  2:147372/105‑1 (MQ=17)
                                                                                                                                            ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTa                                                       >  1:147372/1‑105 (MQ=17)
                                                                                                                                            ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt          >  2:1971833/1‑150 (MQ=25)
                                                                                                                                            ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt          >  2:2534436/1‑150 (MQ=25)
                                                                                                                                            ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGAGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGt                                                        <  2:2366600/104‑1 (MQ=17)
                                                                                                                                            ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGAGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGt                                                        <  2:2366065/104‑1 (MQ=17)
                                                                                                                                            ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGAGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGt                                                        >  1:2366600/1‑104 (MQ=17)
                                                                                                                                            ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGAGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGt                                                        >  1:2366065/1‑104 (MQ=17)
                                                                                                                                            ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATTTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt          >  1:4543825/1‑150 (MQ=255)
                                                                                                                                            ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt          >  2:4685509/1‑150 (MQ=255)
                                                                                                                                            ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt          >  2:3821057/1‑150 (MQ=255)
                                                                                                                                            ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt          >  1:4902079/1‑150 (MQ=255)
                                                                                                                                            ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt          >  2:3878545/1‑150 (MQ=255)
                                                                                                                                             cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACa                                         >  2:1804101/1‑118 (MQ=21)
                                                                                                                                             cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACa                                         <  1:1804101/118‑1 (MQ=21)
                                                                                                                                             cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc         <  2:3629867/150‑1 (MQ=255)
                                                                                                                                             cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc         <  2:1263963/150‑1 (MQ=255)
                                                                                                                                              gcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCg        >  1:4625652/1‑150 (MQ=255)
                                                                                                                                                 aaTACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCa     <  2:3880113/150‑1 (MQ=255)
                                                                                                                                                    aCGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACga  <  2:3199664/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GGCGCCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGA  >  NC_002947/5899600‑5899897

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: