Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,897,829 G→A 11.1% intergenic (‑173/+61) 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,897,8290GA11.1% 212.4 / 11.2 79intergenic (‑173/+61)PP_5172/PP_5173hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base G (37/33);  new base A (6/3);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGA  >  NC_002947/5897680‑5897978
                                                                                                                                                     |                                                                                                                                                     
cgcCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACgg                                                                                                                                                       <  2:2671206/150‑1 (MQ=255)
cgcCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACgg                                                                                                                                                       <  2:1907860/150‑1 (MQ=255)
 gcCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggt                                                                                                                                                      >  2:2258406/1‑150 (MQ=255)
  cccggCAGAGGTTTGCAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGAGTTTACCCGCGAATACggtc                                                                                                                                                     <  1:442008/147‑2 (MQ=255)
   tgattggcgggatcatcgtggcgaccttgctgacgctgctgttcctgccggcgctgtacgtggcctggtacaagatccgcgagccgcagaaacacgcctgaaggcccatatggcccCTGGTAGGAGCGGGTTTACCCGCGAATACggtgg                                                                                                                                                    <  2:1436848‑M1/34‑1 (MQ=255)
        agagGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGAt                                                                                                                                               <  2:2643614/150‑1 (MQ=255)
            gTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCAc                                                                                                                                           >  2:910768/1‑150 (MQ=255)
            gTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCAc                                                                                                                                           <  2:1519358/150‑1 (MQ=255)
               tcatcgtggcgaccttgctgacgctgctgttcctgccggcgctgtacgtggcctggtacaagatccgcgagccgcagaaacacgcctgaaggcccatatggcccCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCAt                                                                                                                                        <  2:2654971‑M1/46‑1 (MQ=255)
                gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                       <  1:2988310/150‑1 (MQ=255)
                 aaGGGGTTTTGCGCCTGGCTTGAGGGCACTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCACGAATACGGCGGCCGATGCACCATcc                                                                                                                                      >  2:19390/1‑149 (MQ=255)
                        cgaccttgctgacgctgctgttcctgccggcgctgtacgtggcctggtacaagatccgcgagccgcagaaacacgcctgaaggcccatatggcctCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                               >  1:750711‑M1/96‑150 (MQ=255)
                        cgaccttgctgacgctgctgttcctgccggcgctgtacgtggcctggtacaagatccgcgagccgcagaaacacgcctgaaggcccatatggcccCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                               >  1:2477830‑M1/96‑150 (MQ=255)
                        cgaccttgctgacgctgctgttcctgccggcgctgtacgtggcctggtacaagatccgcgagccgcagaaacacgcctgaaggcccatatggcacCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                               >  1:2522693‑M1/96‑150 (MQ=255)
                         ttGCCCCTGGCTTGTGTGCAATAGGAGTGAGGCGGGATTATTCTTAACAATACTGTTTAAGAATTTTTATATAACCATTTAAGCGAACCCCCCTCAGGTACGATCATCTTTTCCCGCGACTGCTATGGTGCATCGACCACCGTATTcgcg                                                                                                                              <  1:2681304/150‑1 (MQ=255)
                          tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  1:1680486/150‑1 (MQ=255)
                          tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  1:19390/150‑1 (MQ=255)
                          tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  1:910768/150‑1 (MQ=255)
                          accttgctgacgctgctgttcctgccggcgctgtacgtggcctggtacaagatccgcgagccgcagaaacacgcctgaaggcccatatggcccCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             >  1:2654971‑M1/94‑150 (MQ=255)
                              ccTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggtg                                                                                                                                                     >  1:2147688/1‑122 (MQ=255)
                              ccTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggtg                                                                                                                                                     <  2:2147688/122‑1 (MQ=255)
                                    ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                             >  1:157408/1‑140 (MQ=25)
                                    ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                             >  1:1244345/1‑140 (MQ=25)
                                    ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGc                                                                                                                   <  1:2456189/150‑1 (MQ=21)
                                    ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGc                                                                                                                   <  1:553968/150‑1 (MQ=255)
                                    ttGAGGGCAATAGTAGTGAGGCGGGATCATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                             >  1:825324/1‑140 (MQ=21)
                                     gctgctgttcctgccggcgctgtacgtggcctggtacaagatccgcgagccgcagaaacacgcctgaaggcccatatggcccCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             >  1:1709535‑M1/83‑139 (MQ=255)
                                      gAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGCGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTc                                                                                                                 <  2:1806810/150‑1 (MQ=255)
                                                gtagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAgggg                                                                                                       <  1:2258406/150‑1 (MQ=255)
                                                       ggCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATATAAGAGAAGACGCCTCTGGTAGGAGCGGGTTTACACGCGAATACGGTGGTCGAGGCACCATCGCATGCGCGGGTAAACCCGCTCCTACCAGGGGCCATATg                                                                                                <  2:180128/150‑1 (MQ=38)
                                                        gCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggt                                                                                                                                                      <  2:596332/95‑1 (MQ=255)
                                                        gCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggt                                                                                                                                                      >  1:596332/1‑95 (MQ=255)
                                                                 ttCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                       <  2:2942380/101‑1 (MQ=255)
                                                                 ttCTTAAAAATACTGGTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                       >  1:2942380/1‑101 (MQ=255)
                                                                    ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCgg                                                                                                                             >  1:1157524/1‑108 (MQ=17)
                                                                    ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTAc                                                                                                             <  2:606056/124‑1 (MQ=255)
                                                                    ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTAc                                                                                                             >  1:606056/1‑124 (MQ=255)
                                                                    ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCgt                                                                                   >  2:1077017/1‑150 (MQ=255)
                                                                    ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  1:2917686/108‑1 (MQ=255)
                                                                    ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             >  2:2917686/1‑108 (MQ=255)
                                                                      tgacaagatccgcgagccgcagaaacacgcctgaaggcccatatggcccCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             >  2:2638250‑M1/50‑106 (MQ=255)
                                                                      gtacaagatccgcgagccgcagaaacacgcctgaaggcccatatggcccCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  1:2638250‑M1/57‑1 (MQ=255)
                                                                         caagatccgcgagccgcagaaacacgcctgaaggcccatatggcccCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  1:453323‑M1/57‑1 (MQ=255)
                                                                         caagatccgcgagccgcagaaacacgcctgaaggcccatatggcccCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             >  2:453323‑M1/47‑103 (MQ=255)
                                                                                 gcgagccgcagaaacacgcctgaaggcccatatggcccCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             >  1:1647622‑M1/39‑95 (MQ=255)
                                                                                          tttATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCtt                                                             <  1:46280/150‑1 (MQ=255)
                                                                                               ataACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTAcc                                                        >  2:2272505/1‑150 (MQ=255)
                                                                                                 aaCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggtggt                                                                                                                                                   <  1:323588/57‑1 (MQ=255)
                                                                                                 aaCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggtggt                                                                                                                                                   >  2:323588/1‑57 (MQ=255)
                                                                                                 aaCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                      >  1:2772254/1‑150 (MQ=255)
                                                                                                 aaCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                      >  2:1147089/1‑150 (MQ=255)
                                                                                                  aCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAgg                                                     <  1:1147089/150‑1 (MQ=255)
                                                                                                    cATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGcc                                                   >  2:1545991/1‑150 (MQ=255)
                                                                                                    cATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGcc                                                   >  1:3019463/1‑150 (MQ=255)
                                                                                                       ttAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACg                                                >  1:223877/1‑150 (MQ=255)
                                                                                                           gagaACACGCCTCTGGTAGGAGCGGATTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGCTCTTGTACCAGGCCACGTACa                                            >  2:1548664/1‑150 (MQ=255)
                                                                                                                 acGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCg                                      >  1:2908286/1‑150 (MQ=255)
                                                                                                                   gCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGc                                    >  2:2454970/1‑150 (MQ=255)
                                                                                                                    ccTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCa                                                                                                           <  2:1556962/78‑1 (MQ=255)
                                                                                                                    ccTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCa                                                                                                           >  1:1556962/1‑78 (MQ=255)
                                                                                                                     ctctGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                  <  2:223877/150‑1 (MQ=255)
                                                                                                                      tctGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAgg                                 >  1:2156407/1‑150 (MQ=255)
                                                                                                                      tctGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                  <  2:2372664/149‑1 (MQ=255)
                                                                                                                      tctGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                  >  1:2372664/1‑149 (MQ=255)
                                                                                                                           TAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                            <  1:1436848‑M1/150‑1 (MQ=255)
                                                                                                                             GGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCtt                                                             >  1:903497‑M1/1‑115 (MQ=255)
                                                                                                                             GGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCtt                                                             <  2:903497‑M1/115‑1 (MQ=255)
                                                                                                                               AGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                        >  2:2772254‑M1/1‑150 (MQ=255)
                                                                                                                                         CCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg              <  2:327093‑M1/150‑1 (MQ=255)
                                                                                                                                          CCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgc                                                                    >  2:1647622‑M1/1‑95 (MQ=255)
                                                                                                                                          CCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                        >  2:1709535‑M1/1‑139 (MQ=255)
                                                                                                                                          CCGCGAATACGGTGGTCGATGCACCGT‑‑‑ATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgc                                                                    >  2:2243690‑M1/1‑92 (MQ=255)
                                                                                                                                          CCGCGAATACGGTGGTCGATGCACCGT‑‑‑ATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgc                                                                    <  1:2243690‑M1/92‑1 (MQ=255)
                                                                                                                                          CCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACa                                            <  1:245174‑M1/119‑1 (MQ=255)
                                                                                                                                          CCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACa                                            >  2:245174‑M1/1‑119 (MQ=255)
                                                                                                                                          CCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt             >  1:2454970‑M1/1‑150 (MQ=255)
                                                                                                                                          CCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt             <  1:1316410‑M1/150‑1 (MQ=255)
                                                                                                                                          CCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt             <  2:1146866‑M1/150‑1 (MQ=255)
                                                                                                                                          CCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt             <  2:60163‑M1/150‑1 (MQ=255)
                                                                                                                                          CCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt             <  1:2935653‑M1/150‑1 (MQ=255)
                                                                                                                                           CGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc            <  2:2457579‑M1/150‑1 (MQ=255)
                                                                                                                                           CGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc            <  1:1980241‑M1/150‑1 (MQ=255)
                                                                                                                                           CGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc            >  2:2935653‑M1/1‑150 (MQ=255)
                                                                                                                                            GCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACcagag                                                                                                        >  2:2210207‑M1/1‑55 (MQ=255)
                                                                                                                                            GCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACcagag                                                                                                        <  1:2210207‑M1/57‑3 (MQ=255)
                                                                                                                                             CGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGc          >  1:1077017‑M1/1‑150 (MQ=255)
                                                                                                                                               AATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAgg                                 >  2:68970‑M1/1‑125 (MQ=255)
                                                                                                                                               AATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAgg                                 <  1:68970‑M1/125‑1 (MQ=255)
                                                                                                                                                ATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCAc       >  2:2156407‑M1/1‑150 (MQ=255)
                                                                                                                                                  aCGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACga     >  1:437420/1‑150 (MQ=255)
                                                                                                                                                    ggtggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgatg   >  1:629557/1‑150 (MQ=255)
                                                                                                                                                     gtggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGAATAAACCGCGCCGCGACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgatga  >  2:1951125/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGA  >  NC_002947/5897680‑5897978

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