Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,130,096 C→A 5.8% intergenic (+73/‑109) PP_4511 → / → PP_4512 AraC family transcriptional regulator/DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,130,0960CA 5.8% 688.0 / 6.9 212intergenic (+73/‑109)PP_4511/PP_4512AraC family transcriptional regulator/DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base C (97/97);  new base A (14/3);  total (112/100)
Fisher's exact test for biased strand distribution p-value = 1.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

GTGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGC  >  NC_002947/5129947‑5130245
                                                                                                                                                     |                                                                                                                                                      
gggctggcgctggtggccagtttgccgaagacgccggtaaggccgaagaacagcgcgccaatgtgcagtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagtgccggccctTTCGCGGGTAAAc                                                                                                                                                        <  2:5280630‑M1/13‑1 (MQ=255)
gTGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAC                                                                                                                                                        <  2:358492‑M1/150‑1 (MQ=255)
 tGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACC                                                                                                                                                       <  1:862702‑M1/150‑1 (MQ=255)
   ctggcgctggtggccagtttgccgaagacgccggtaaggccgaagaacagcgcgccaatgtgcagtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagtgccggccctTTCGCGGGTAAACCCg                                                                                                                                                     >  1:5544752‑M1/135‑150 (MQ=255)
   ctggcgctggtggccagtttgccgaagacgccggtaaggccgaagaacagcgcgccaatgtgcagtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagtgccggccctTTCGCGGGTAAACCCg                                                                                                                                                     >  1:5285567‑M1/135‑150 (MQ=255)
    cTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGC                                                                                                                                                    <  1:121933‑M1/150‑1 (MQ=255)
     ggcgctggtggccagtttgccgaagacgccggtaaggccgaagaacagcgcgccaatgtgcagtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagtgccggccctTTCGCGGGTAAACCCGCt                                                                                                                                                   >  1:3840948‑M1/133‑150 (MQ=255)
      tttCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTC                                                                                                                                                  <  1:281380‑M1/150‑1 (MQ=255)
       ttCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTTACCCGCgaa                                                                                                                                                 <  1:4027884‑M1/150‑4 (MQ=255)
       ttCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCC                                                                                                                                                 <  2:2383779‑M1/150‑1 (MQ=255)
        tCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTTACCCGCgaan                                                                                                                                                <  1:3764194‑M1/150‑5 (MQ=255)
        tCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTTACCCGCgaag                                                                                                                                                <  1:640647‑M1/150‑5 (MQ=255)
        tCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTTACCCGCgaag                                                                                                                                                <  2:4049453‑M1/150‑5 (MQ=255)
           gcggcgATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACA                                                                                                                                             <  2:4547930‑M1/150‑1 (MQ=255)
           gcggcgATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACA                                                                                                                                             >  2:640339‑M1/1‑150 (MQ=255)
              gcgATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGa                                                                                                                                          <  1:1045915/150‑1 (MQ=255)
                 aTGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGc                                                                                                                                       <  2:1895611/150‑1 (MQ=255)
                            gTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGaa                                                                                                                            <  2:5118564/150‑1 (MQ=255)
                             ttCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAt                                                                                                                           <  2:945761/150‑1 (MQ=255)
                             ttCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAt                                                                                                                           <  2:3126360/150‑1 (MQ=255)
                             ttCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAt                                                                                                                           <  2:4183501/150‑1 (MQ=255)
                              tCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtg                                                                                                                          <  2:1500093/150‑1 (MQ=255)
                                 ggTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttg                                                                                                                       <  2:885990/150‑1 (MQ=255)
                                  gtgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                      <  1:1943407/150‑1 (MQ=255)
                                  gtgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgttgt                                                                                                                      <  1:4925546/150‑1 (MQ=255)
                                   tgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTtgtg                                                                                                                    <  1:629378/150‑1 (MQ=255)
                                   tgtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCTACAGGATGCCCCCGCCGGAATGTtgtg                                                                                                                    >  2:4502541/1‑150 (MQ=255)
                                    gtACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGc                                                                                                                   <  1:4228503/150‑1 (MQ=255)
                                        acCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCaaaa                                                                                                               <  2:560184/150‑1 (MQ=255)
                                        acCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCCCGCCCACACCGTGCCCCCCCCAGAAAGTTGTGGaaaa                                                                                                               >  1:3906836/1‑150 (MQ=255)
                                         ccTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAc                                                                                                              >  1:3538270‑M2/1‑149 (MQ=255)
                                          cTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcc                                                                                                             >  1:1365941‑M2/1‑148 (MQ=255)
                                           tttGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcct                                                                                                            >  1:661546‑M2/1‑147 (MQ=255)
                                           tttGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcct                                                                                                            >  2:4341520‑M2/1‑147 (MQ=255)
                                           tttGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCCCCTGAATGTTGTGCAAAAcag                                                                                                            <  2:3993505‑M2/150‑4 (MQ=255)
                                           tttGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCCCGCCTACACGATTCCCCCCCCGGAAAGTTGTGCAAAAcca                                                                                                            <  2:2366686‑M2/150‑4 (MQ=255)
                                             tGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgt                                                                                                          >  1:589380‑M2/1‑145 (MQ=255)
                                             tGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgt                                                                                                          >  1:4627410‑M2/1‑145 (MQ=255)
                                              ggCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgta                                                                                                         >  1:793056‑M2/1‑144 (MQ=255)
                                              ggCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgta                                                                                                         >  1:2320348‑M2/1‑144 (MQ=255)
                                              ggCTTTGCGGCGAGAGTCTGGCGAGAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgta                                                                                                         >  2:4349823‑M2/1‑144 (MQ=255)
                                              agaacagcgcgccaatgtgcagtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagtgccggccccTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTa                                                                                                         >  2:2481129‑M1/92‑150 (MQ=255)
                                               gaacagcgcgccaatgtgcagtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagtgccggccctTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAg                                                                                                        >  2:560769‑M1/91‑150 (MQ=255)
                                                    gcggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCTACAGGATGCCCCCGCCGGAATGTTGTGCAAAAcctgtaggagcg                                                                                                   <  2:4719748‑M2/150‑13 (MQ=255)
                                                     cggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTAAACCGCTCCTACAGGTTTCCCCCGGCGGGATGTTGTGCAAAAagtgtaggttcgg                                                                                                  <  2:125171‑M2/150‑14 (MQ=255)
                                                     cggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGttttgcacaacattcaggcgggggcatcctgtaggagcgg                                                                                                   <  2:141237‑M1/150‑41 (MQ=255)
                                                     cggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgtaggagcgg                                                                                                  >  1:4719295‑M2/1‑137 (MQ=255)
                                                     cggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgtaggagcgg                                                                                                  >  1:4075520‑M2/1‑137 (MQ=255)
                                                     cgcgccaatgtgcagtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagtgccggccccTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCgg                                                                                                  >  1:3345230‑M1/85‑150 (MQ=255)
                                                       gcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGttttgcacaacattcaggcgggggcatcctgtaggagcgggt                                                                                                 <  1:4606134‑M1/150‑43 (MQ=255)
                                                       gcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgtaggagcgggt                                                                                                >  1:2896450‑M2/1‑135 (MQ=255)
                                                        cgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACCGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgtaggagcgggtt                                                                                               >  1:1098224‑M2/1‑134 (MQ=255)
                                                            agTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgtaggagcgggtttacc                                                                                           >  1:521129‑M2/1‑130 (MQ=255)
                                                            agTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgtaggagcgggtttacc                                                                                           >  1:3359535‑M2/1‑130 (MQ=255)
                                                             gTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgtaggagcgggtttaccc                                                                                          >  2:5459805‑M2/1‑129 (MQ=255)
                                                             gTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGATTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAcctgtaggagcgggtttaccc                                                                                          >  2:5459829‑M2/1‑129 (MQ=255)
                                                                 ggCGATAACAGACAGCAGTCTTGTGACAGAATGTGTTAGCAAGAGGATTTAAAATTTGACTATGCCCGCCCCTTCGACGGTAAACCGCTCCTACAGGATGCCCACGCCTGAATGTTGTGCAAAAACAGTAGGAGAGGGTGTACCAGCGaa                                                                                     >  2:2757482/1‑150 (MQ=255)
                                                                  agtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagtgccggccctTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGt                                                                                                          <  1:2188850‑M1/58‑1 (MQ=255)
                                                                    gATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGttttgcacaacattcaggcgggggcatcctgta                                                                                                          >  2:5228889‑M1/1‑95 (MQ=255)
                                                                    gATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATTCCCCCGCCTGGATGTTGTGCAAAACCGGTAGGCGCGGGTTTACCCGCGAAAg                                                                                   >  2:427220/1‑150 (MQ=255)
                                                                      cgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagtgccggccctTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGt                                                                                                          <  2:4634215‑M1/58‑1 (MQ=255)
                                                                         cTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCGGAATGTTGTGCAAAACCGGTAGGGGCGGGGTTACCCGCGGAAGGGCCg                                                                              >  2:2225885/1‑150 (MQ=39)
                                                                         cTGACAGCAGTCTTGCGACAGAAGTTGTTTGCCGGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGAGGCCCCCGCCTGAATGTGGTGCAAACCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCg                                                                              <  2:735348/150‑1 (MQ=38)
                                                                          tGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTC                                                                                                                                                  <  1:3061119‑M1/82‑1 (MQ=255)
                                                                                 cgTCTTGCGACCGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTGTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTCCAGGATGCCCCCGCCTGAATTTTGTGCAAACCCTGTTGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGcc                                                                      <  2:5493104/149‑1 (MQ=38)
                                                                                 agTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCTACACGATGCCCCCCCCGGAATGTTGTGCAAAACCGGTAGGAGCGGGGTTACCCGCGAAAGGGCCGGCCCTGcc                                                                      >  2:5006547/1‑150 (MQ=38)
                                                                                   tCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTAAACGCTGCCCCCCCCTGAATGTTGTGCAAAACCGGTAGGTGCGGGTTTACCCCCCAAAGGGCCGGCACAGccgg                                                                    >  2:3896566/1‑148 (MQ=255)
                                                                                   tCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGGTAACCCCCGCCGACACGCTGCCCCCCCCGGGAAGTTGTGGAAAAACCGGTGGGGCCGGGTTTCCCCCGGAATGGCCCGCCCTggcgg                                                                    >  1:3579608/1‑146 (MQ=255)
                                                                                     ttGCGAAAGAATTTGTTTGACTGAGGCTTTAAACTTTTCCTTTGCCGGCCCCTTCGCGGGTAAAAACCCACCGAAACGATGCCCCCCCCAGAATTTTTTTCAAAAAAAGTATGTGAGGGTTTAAACCCCCAAGGGCCCGCCCTGCCAGTa                                                                  >  2:3837749/1‑150 (MQ=255)
                                                                                         gACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgtt                                                                                                                        >  1:2441468/1‑93 (MQ=255)
                                                                                          aCAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTAAACCGCGCCAACAGGATGCCCCCCCCGGGAAGTTGTGCAAAAACGGGAGGCGCGGGTTTTCCCCCGGAAGGGCCGGCCCAGCCGGTAGatat                                                             >  1:5077822/1‑150 (MQ=255)
                                                                                          aCAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCGCTCCTACAGGCTGCCCCCGCCGGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGCAAGGGGCGGCCCTGCCGGTAGAGAt                                                             >  2:3034392/1‑150 (MQ=38)
                                                                                            aGAATGTTTTTGGCCTATGCATTACCCTTTGGCTTTTCCGGCCCCCTCGGGGTTAAACCCGCCCCAAACGGTTGCCCCCGGCTTTATTTTGGTCAAAACCTTTAGGATGGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGAtct                                                           <  2:2524082/150‑1 (MQ=255)
                                                                                                    tttGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACCGGATTCCCCCGCCTGAATGTTGTGCAAAAACGGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAg                                                   >  2:3679350/1‑150 (MQ=38)
                                                                                                    tttGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCGACAGGATGCCCCCGCCGGAATGTTGTGCCAAACCTGGAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAg                                                   >  2:4092287/1‑150 (MQ=255)
                                                                                                      tGCCTGAGGCTTTACACTTTGCCTCTCCCGGCCCCTACCCGGGTAAACCCGCTCCTACAGCATGCCCCCCCCTGAATGTTGCCCAAAACCTGAAGCAGCGGGTTAACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGa                                                 <  2:3343135/150‑1 (MQ=255)
                                                                                                          tgtggCTTTACACTTTGCTTTTGCGGGCCCCTTCGCGGGTATACCCGTTCCTACTGTATGCCGCCGCCTGTATTTTGTGCAAAACCTTTTGGTGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACcca                                             <  2:2256178/147‑1 (MQ=255)
                                                                                                          tGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCAACAGGCTGCCCCCGCCTGAATGTTGTGCAAAACCTGTGGGAGCGGGGGTAACCCCGAAAGGGCCCGCACTGCCAGTCGAGATATCCATCAAGGAcccc                                             >  2:5456702/1‑149 (MQ=37)
                                                                                                                 ttACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCCCTCCTAAAGGATTCCCCCGCCAGAATGTTGTGCAAAACCAGTAGGCGCGGGTTTAACCCCGAAAGGGCCCGCCCAGCCGGGAGAGAACTCCATCAACGCCCCACCAAGaa                                      >  2:3614782/1‑150 (MQ=255)
                                                                                                                  tACACTTTGCCTGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCCACAGGAGGCCCCCGCCTGATTGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAac                                     <  2:3061562/150‑1 (MQ=38)
                                                                                                                    cacTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCCCGCCTACAGGCTGCCCCCCCCGGAATGTTGTGGCAAAACGGGCGGGGCGGGTTTACCCCCGCAAGGGCCCGCCCGGCCGGGGGAGCGCTCCAGCCAGGGCCCCCCCGGAaccc                                   >  2:295548/1‑149 (MQ=255)
                                                                                                                       tttgcctttgccGGCCCCTTCGCGGGGCAAACCCCGCCGACCCGAAGCCCCCCCCCGGATGTTGTGCCAAACCGGGGGGGGCCGGGTTACCCCCGCAACGGCCGGCCCAGCCGGGAGAGATCTCCCAAAAGGGCCCACCAAGAcccaccc                                >  2:232846/1‑149 (MQ=255)
                                                                                                                         tgcctttgccGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGTTGCCCCCGCATGAATTTTGTGCAAAACCTGTAGGTGCGGTTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACcaccgca                              <  2:1157396/150‑1 (MQ=38)
                                                                                                                             tttgccGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGCATGTTGTGCAAAACCTGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGCACGGCCAGTAGAGATCTCCATCCAGGACCCACCAAGAACCaccgcaccgc                          >  2:5051221/1‑150 (MQ=38)
                                                                                                                                  cGGCCCCTGCGCGGGAAAACCCGCTCCTCCAGGTTGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCg                     <  2:2934794/150‑1 (MQ=255)
                                                                                                                                        cTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAAACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCAc               >  2:770638/1‑150 (MQ=255)
                                                                                                                                        cTTCGCGGGTAAACCCGCACCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCGGGAGGAGCGGGGTTACCCGCGAAAGGGCCGGCACAGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCAc               >  2:1923238/1‑150 (MQ=255)
                                                                                                                                          TCGCGGTTAACCCCGCTCCTACAGGATGCCCCCGCCTGATTGTTGTGCAAATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTg             >  2:2859451‑M1/1‑150 (MQ=255)
                                                                                                                                            GCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGca           <  2:4975206‑M1/150‑1 (MQ=255)
                                                                                                                                              gggTAAACCCGCTCCTACAGGttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaggggccggcaaaggcaaagtgtaaagcctcaggcaaacacattctgtcgcaagactgctgtcagttatcgccagac          >  2:4433073‑M1/1‑21 (MQ=255)
                                                                                                                                              gggTAAACCCGCTCCTACAGGttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaggggccggcaaaggcaaagtgtaaagcctcaggcaaacacattctgtcgcaagactgctgtcagttatcgccagac          >  2:4433007‑M1/1‑21 (MQ=255)
                                                                                                                                               GGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACAt        <  1:2601210‑M1/150‑1 (MQ=255)
                                                                                                                                                  AAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTgg     <  2:2099522‑M1/150‑1 (MQ=255)
                                                                                                                                                  AAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTgg     <  1:4864872‑M1/150‑1 (MQ=255)
                                                                                                                                                   aaCCCGCTCCTACAGGttttgcacaacattcaggcgggggcatcctgtaggagcgggtttacccgcgaaggggccggcaaaggcaaagtgtaaagcctcaggcaaacacattctgtcgcaagactgctgtcagttatcgccagactctcg     >  1:450530‑M1/1‑16 (MQ=255)
                                                                                                                                                   AACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgc    <  2:3514612‑M1/150‑1 (MQ=255)
                                                                                                                                                    ACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcg   <  2:1650303‑M1/150‑1 (MQ=255)
                                                                                                                                                    ACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcg   <  1:3709844‑M1/150‑1 (MQ=255)
                                                                                                                                                     CCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcgc  <  2:4796004‑M1/150‑1 (MQ=255)
                                                                                                                                                     CCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGgcgc  <  1:5196089‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCACATTGGCGC  >  NC_002947/5129947‑5130245

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: