Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,185,009 C→G 8.3% intergenic (‑240/‑43) PP_4557 ← / → PP_4558 hypothetical protein/GNAT family acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,185,0090CG 8.3% 874.1 / 52.5 264intergenic (‑240/‑43)PP_4557/PP_4558hypothetical protein/GNAT family acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (141/101);  new base G (10/12);  total (151/113)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 273 total aligned reads displayed.

ACGAAGTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTGGGTAGC  >  NC_002947/5184860‑5185158
                                                                                                                                                     |                                                                                                                                                     
aCGAAGTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGc                                                                                                                                                       >  1:4363067/1‑150 (MQ=255)
  gAAGTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGCCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCct                                                                                                                                                     >  1:3628504/1‑150 (MQ=255)
     gTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTg                                                                                                                                                  <  1:2502805/150‑1 (MQ=255)
     gTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTg                                                                                                                                                  <  1:1932514/150‑1 (MQ=255)
       aCCGCTCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGc                                                                                                                                                <  2:4392661/150‑1 (MQ=255)
        ccGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGcc                                                                                                                                               >  2:1074959/1‑150 (MQ=255)
           aTCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGCCtga                                                                                                                                            <  1:3769671/150‑3 (MQ=255)
           aTCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCg                                                                                                                                            >  2:3331932/1‑150 (MQ=255)
           aTCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCg                                                                                                                                            >  2:1434846/1‑150 (MQ=255)
            tCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGCCtgat                                                                                                                                           <  2:3020154/150‑4 (MQ=255)
            tCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGCCtgat                                                                                                                                           <  2:3591834/150‑4 (MQ=255)
            tCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGt                                                                                                                                           >  1:1321616/1‑150 (MQ=255)
            tCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTCGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGt                                                                                                                                           >  2:1182889/1‑150 (MQ=255)
             cccGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGtt                                                                                                                                          <  1:2271373/150‑1 (MQ=255)
               cGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTc                                                                                                                                        >  1:4122749/1‑150 (MQ=255)
                gCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCa                                                                                                                                       >  1:758764/1‑150 (MQ=255)
                gCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCa                                                                                                                                       >  1:3309967/1‑150 (MQ=255)
                gCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCa                                                                                                                                       >  2:4575703/1‑150 (MQ=255)
                gCTTCGTACTGCCAGGGGGGATGAACAAAAGTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGGTTTGCTACAATCTGCGACCAAGAGGCCGCAAATCATCAGGCCTGAGCGGCTGGCcgcgggggg                                                                                                                                       >  1:4508837/1‑145 (MQ=255)
                 cTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAg                                                                                                                                      <  2:5472830/150‑1 (MQ=255)
                    cGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGt                                                                                                                                   >  2:39579/1‑150 (MQ=255)
                        cTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAggg                                                                                                                               <  1:5620261/150‑1 (MQ=255)
                           ccAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTCGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTc                                                                                                                            <  2:1658958/150‑1 (MQ=255)
                             aGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCAc                                                                                                                          >  1:4112543/1‑150 (MQ=255)
                                   gATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTc                                                                                                                    <  2:5040012/150‑1 (MQ=255)
                                    aTGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCt                                                                                                                   <  1:2445123/150‑1 (MQ=255)
                                    aTGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCt                                                                                                                   <  2:170372/150‑1 (MQ=255)
                                    aTGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTACAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCt                                                                                                                   <  1:2796815/150‑1 (MQ=255)
                                        aCAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAAc                                                                                                               >  1:4908949/1‑150 (MQ=255)
                                         cAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACa                                                                                                              >  1:3712489/1‑150 (MQ=255)
                                         cAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACa                                                                                                              >  1:3107906/1‑150 (MQ=255)
                                          aaaaTTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACaa                                                                                                             <  2:2740981/150‑1 (MQ=255)
                                          aaaaTTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACaa                                                                                                             <  1:3902618/150‑1 (MQ=255)
                                            aaTTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCCCCATTTCTGAACAAtt                                                                                                           >  1:1611231/1‑150 (MQ=255)
                                                   cTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGt                                                                                                    <  2:883884/150‑1 (MQ=255)
                                                    tCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCGTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGTACAATTGTACTGtt                                                                                                   >  2:4308403/1‑150 (MQ=255)
                                                     cccGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTc                                                                                                  <  1:3669257/150‑1 (MQ=255)
                                                     cccGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTc                                                                                                  <  2:4536444/150‑1 (MQ=255)
                                                      ccGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCAGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCg                                                                                                 >  1:448137/1‑150 (MQ=255)
                                                        ggCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGAc                                                                                               <  1:788394/150‑1 (MQ=255)
                                                              cAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGGCTCGTTTGAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTc                                                                                         >  1:1737203/1‑150 (MQ=255)
                                                              cAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTc                                                                                         <  1:5548962/150‑1 (MQ=255)
                                                               aaTAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACCATTGTACTGTTCGACGCTTTCg                                                                                        >  1:4292540/1‑150 (MQ=255)
                                                               aaTAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCg                                                                                        >  1:5477069/1‑150 (MQ=255)
                                                               aaTAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCAGTTCTGAACAATTGTACTGTTCGACGCTTTCg                                                                                        >  1:351609/1‑150 (MQ=255)
                                                                   gTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGaca                                                                                    <  2:2493353/150‑1 (MQ=255)
                                                                    ttGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAAGTGTACTGTTCGACGCTTTCGGacac                                                                                   >  2:5152715/1‑150 (MQ=255)
                                                                      gAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCACAGGCCTCGTATCAGGGTAGGGCTCACCATTTCTGAACAATAGTACTGTTCGACGCATTCGGAcacaa                                                                                 >  2:3274484/1‑149 (MQ=255)
                                                                           ccGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGc                                                                            <  1:1845374/150‑1 (MQ=255)
                                                                           ccGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGc                                                                            <  2:2958840/150‑1 (MQ=255)
                                                                                   ggATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCt                                                                    <  1:501915/150‑1 (MQ=255)
                                                                                       atTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGAc                                                                <  2:4683767/150‑1 (MQ=255)
                                                                                       atTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGAc                                                                <  1:1056962/150‑1 (MQ=255)
                                                                                         tCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCa                                                              >  1:1445132/1‑150 (MQ=255)
                                                                                           aCCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACt                                                            >  2:2316902/1‑150 (MQ=255)
                                                                                           aCCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACt                                                            >  2:1908832/1‑150 (MQ=255)
                                                                                           aCCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGATTGCTGGACCACt                                                            <  2:707192/150‑1 (MQ=255)
                                                                                             ccTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTc                                                          >  1:1216753/1‑150 (MQ=255)
                                                                                             ccTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTc                                                          >  1:2596022/1‑150 (MQ=255)
                                                                                              cTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCt                                                         >  2:2005500/1‑150 (MQ=255)
                                                                                               ttGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTAGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTa                                                        >  1:2941485/1‑150 (MQ=255)
                                                                                                tGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTAc                                                       <  2:2089609/150‑1 (MQ=255)
                                                                                                tGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTAc                                                       <  2:4920866/150‑1 (MQ=255)
                                                                                                     gTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCa                                                  >  2:5611957/1‑150 (MQ=255)
                                                                                                     gTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCa                                                  >  2:528538/1‑150 (MQ=255)
                                                                                                          gCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGcaca                                             <  2:2446451/150‑1 (MQ=255)
                                                                                                            aacaaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGc                                           >  2:1459755/1‑150 (MQ=255)
                                                                                                            aacaaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGc                                           >  1:4017427/1‑150 (MQ=255)
                                                                                                            aacaaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGc                                           >  2:3874614/1‑150 (MQ=255)
                                                                                                            aacaaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACCTCTACCGCCAGCACAGc                                           >  1:1872007/1‑150 (MQ=255)
                                                                                                                aTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCAc                                       <  1:2019843/150‑1 (MQ=255)
                                                                                                                aTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCAc                                       <  2:415538/150‑1 (MQ=255)
                                                                                                                 tCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGGACAATTGTACTGTTCGGCGCTTTCGGACACAGAACGCCGTGTGCTGGACCACTTCTACCGCCAGCACAGCTCACg                                      >  1:689177/1‑150 (MQ=255)
                                                                                                                   tGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCCGCACAGCTCACGCa                                    >  1:4694684/1‑150 (MQ=255)
                                                                                                                    gCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCAt                                   >  1:5137156/1‑150 (MQ=255)
                                                                                                                    gCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCAt                                   >  2:3768572/1‑150 (MQ=255)
                                                                                                                     cGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATg                                  >  2:4043010/1‑150 (MQ=255)
                                                                                                                      gACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATgc                                 >  2:2864487/1‑150 (MQ=255)
                                                                                                                         aaaGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATgcgcg                              >  1:5349818/1‑150 (MQ=255)
                                                                                                                          aaGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTgg                                                                  <  1:4175780/113‑1 (MQ=255)
                                                                                                                               gCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCaa                        <  1:4812324/150‑1 (MQ=255)
                                                                                                                               gCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCaa                        <  2:3628504/150‑1 (MQ=255)
                                                                                                                               gCCGCAAATCATCAGACCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCaa                        <  1:3702421/150‑1 (MQ=255)
                                                                                                                                 cGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGc                      <  2:2221282/150‑1 (MQ=255)
                                                                                                                                  gCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCg                     >  2:5077736/1‑150 (MQ=255)
                                                                                                                                     tgaaacgaGGCCAGAGGCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACg                  >  2:3396333/8‑150 (MQ=255)
                                                                                                                                        aacgaGGCCAGAGGCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCg               >  2:4899704/5‑150 (MQ=255)
                                                                                                                                         atcaGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGa              >  1:5417549/1‑150 (MQ=255)
                                                                                                                                         atcaGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGa              >  1:2084134/1‑150 (MQ=255)
                                                                                                                                           gaGGCCAGAGGCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAAc            >  1:3320077/2‑150 (MQ=255)
                                                                                                                                           gaGGCCAGAGGCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAAc            >  2:5184926/2‑150 (MQ=255)
                                                                                                                                              gCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACtgt         >  1:3557394/1‑150 (MQ=255)
                                                                                                                                              gCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACtgt         <  2:4122749/150‑1 (MQ=255)
                                                                                                                                               ccAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACtgtg        >  2:5348185/1‑150 (MQ=255)
                                                                                                                                                cAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTgg       <  2:3309967/150‑1 (MQ=255)
                                                                                                                                                cAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTgg       >  1:943481/1‑150 (MQ=255)
                                                                                                                                                 agagCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTggg      <  2:353837/150‑1 (MQ=255)
                                                                                                                                                 agagCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTggg      <  2:1559760/150‑1 (MQ=255)
                                                                                                                                                  gagCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTGGGt     <  2:1206603/150‑1 (MQ=255)
                                                                                                                                                     ccTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTGGGTAGc  >  1:519188/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGAAGTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTGGGTAGC  >  NC_002947/5184860‑5185158

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