Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,639,555 A→T 9.7% intergenic (+25/+35) creD → / ← arcA inner membrane protein/response regulator in two‑component regulatory system with ArcB or CpxA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,639,5550AT 9.7% 334.9 / 20.4 123intergenic (+25/+35)creD/arcAinner membrane protein/response regulator in two‑component regulatory system with ArcB or CpxA
Reads supporting (aligned to +/- strand):  ref base A (62/49);  new base T (7/5);  total (69/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

Only 100 of 123 total aligned reads displayed.

TGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACG  >  NC_000913/4639412‑4639701
                                                                                                                                               |                                                                                                                                                   
tggCTCTGGCCGGAATGATTTTTGTTACCCGTAATATCTACTGTTATGCGTTTTCACTGCCTAAAATGAAAGCCATAAATGAAGTTACAACGGACGATGATTTACGTTTCTGTAAATAAGGTTTAAAAATAAAAACGGCGCTAAAAAgcg                                                                                                                                               <  1:1306163/150‑1 (MQ=255)
tggCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAgcg                                                                                                                                               <  2:455742/150‑1 (MQ=255)
 ggCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAATCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAgcgc                                                                                                                                              <  1:1021081/150‑1 (MQ=255)
   cgcTGGCCGGCATGATTTTTTTTACACTTAATATCGACTTGTTTGCTTTTTAACTGCCTAAAATTAAAGCCAGTAAAGATTTTACAACGGACGATGTGTTACGTATCTGGAAATAAGGTTGAAAAAAAAAAACGGCGCTAAAAAGCGCCg                                                                                                                                            <  1:878693/150‑1 (MQ=255)
     gtGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAATTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAATACGGCGCTTTTTAGCGCCG‑tt                                                                                                                                         <  2:1227178/149‑1 (MQ=255)
     cTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑tt                                                                                                                                         >  2:68289/1‑150 (MQ=255)
      tGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑ttt                                                                                                                                        >  2:1336364/1‑150 (MQ=255)
        gCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTTTTTAGCGCCG‑ttttt                                                                                                                                      <  1:688767/150‑1 (MQ=255)
         ccGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑tttttt                                                                                                                                     >  2:470967/1‑150 (MQ=255)
         ccGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑tttttt                                                                                                                                     >  1:1430832/1‑150 (MQ=255)
         ccGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑tttttt                                                                                                                                     >  1:1424139/1‑150 (MQ=255)
           ggCGTGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑tttttttt                                                                                                                                   >  2:1618655/1‑150 (MQ=255)
           ggCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑tttttttt                                                                                                                                   >  2:607036/1‑150 (MQ=255)
            gCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCGGGAAAGAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTg                                                                                                                                  >  2:326202/1‑150 (MQ=255)
             cATGAGGTTTGTGACCCGAAAAATCGACTGGTATGCGTTTTAACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGa                                                                                                                                 <  2:1386750/150‑1 (MQ=255)
               tgatgTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACg                                                                                                                               >  1:1027639/1‑150 (MQ=255)
               tgatgTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACg                                                                                                                               >  2:873128/1‑150 (MQ=255)
               tgatgTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCAGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACg                                                                                                                               >  1:358265/1‑150 (MQ=255)
                gatgTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACgg                                                                                                                              >  1:632350/1‑150 (MQ=255)
                 atgTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTTTTTAGCGCCG‑tttttatttttcaa                                                                                                                             <  2:333978/150‑7 (MQ=39)
                    tttGTGACCCGTAATATCGTCTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACggtggt                                                                                                                          <  1:1103857/150‑1 (MQ=255)
                      tgtgACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACTGTTGTaa                                                                                                                        >  2:1406921/1‑150 (MQ=255)
                       gtgACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGTAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTaaa                                                                                                                       >  2:903455/1‑150 (MQ=255)
                       gtgACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAACCGGCGCTAAAAAGCTCCG‑GTTTTTTTGACGGTGGTaaa                                                                                                                       >  1:308479/1‑150 (MQ=255)
                         gACCCGTAATATCGACTGGTATGTGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGc                                                                                                                     >  1:660012/1‑150 (MQ=255)
                         gACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACCACGGACGATGAGTTAAGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGc                                                                                                                     >  1:1302607/1‑150 (MQ=255)
                            ccGTAATATCGACTGGTATGCGTTATCACTGCCCAAAATGAAACCCAGTAACGAAGTTACAACGCACGATGAGTTACGTAGCTGGAAATACGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGa                                                                                                                  <  1:1620253/150‑1 (MQ=255)
                                   atCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCt                                                                                                           >  2:1677420/1‑150 (MQ=255)
                                      gACTGGTATGCGTTTTCAATGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTcc                                                                                                        <  1:820436/150‑1 (MQ=255)
                                       aCTGGTATGTGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCa                                                                                                       <  2:660012/150‑1 (MQ=255)
                                          ggTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGAt                                                                                                    >  1:1261420/1‑150 (MQ=255)
                                           ttATGCTTTTTCACTCCCGAAAATTAAAGCCAGAAAATATGTAAAAACGGACGATGATTTACGTATCTGTAAATAAGTTTTAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATc                                                                                                   <  1:607036/149‑1 (MQ=255)
                                           gTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATc                                                                                                   <  1:1336364/150‑1 (MQ=255)
                                               gCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCg                                                                                               <  1:1618655/150‑1 (MQ=255)
                                                cGTTTTCACTGCCTAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGc                                                                                              <  2:1281427/150‑1 (MQ=255)
                                                cGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGc                                                                                              >  2:1330375/1‑150 (MQ=255)
                                                 tttttCACTGCCTAAAATGAAAGCCATTAAAGAAGTTACAACGGACTATGATTAACGTATCTGGAAATAAGGTTTAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCa                                                                                             <  1:1252830/149‑1 (MQ=255)
                                                 gTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCa                                                                                             >  2:1326875/1‑150 (MQ=255)
                                                  ttttCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAg                                                                                            <  1:972653/150‑1 (MQ=255)
                                                      cACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGc                                                                                        >  1:1264592/1‑150 (MQ=255)
                                                        cTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTATTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGa                                                                                      >  2:900540/1‑150 (MQ=255)
                                                          gCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATa                                                                                    >  2:668050/1‑150 (MQ=255)
                                                          gCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATATTCCAGATCACCGCAGAAGCGATa                                                                                    >  1:1662462/1‑150 (MQ=255)
                                                            cGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAAc                                                                                  >  1:112963/1‑150 (MQ=255)
                                                             gAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAAcc                                                                                 <  1:237941/150‑1 (MQ=255)
                                                               aaaTGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACggtgg                                                                                                                           >  1:1485345/1‑106 (MQ=255)
                                                                  tGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCAc                                                                            >  1:1121947/1‑150 (MQ=255)
                                                                   gAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCAcc                                                                           >  1:455363/1‑150 (MQ=255)
                                                                       tccAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGa                                                                       <  1:1038293/149‑1 (MQ=255)
                                                                                      aCAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGAtt                                                        >  1:1667529/1‑150 (MQ=255)
                                                                                       cAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGAttt                                                       <  1:303193/150‑1 (MQ=255)
                                                                                         ccgcaCGATGAGTTACGTATATTGAAATAAGTTTTAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtcc                                                     <  1:1714243/146‑1 (MQ=255)
                                                                                           ggACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccgg                                                   <  1:999972/150‑1 (MQ=255)
                                                                                           ggACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccgg                                                   <  1:1152838/150‑1 (MQ=255)
                                                                                             aCGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcg                                                 <  2:388596/150‑1 (MQ=255)
                                                                                              cGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgt                                                >  2:239061/1‑150 (MQ=255)
                                                                                                  gAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTTACGGTTGTAAATCCGATTAATAATACAGATAACCGCAGAAGCGAATACCTTAACCGTGAATTGTTGCGTTTAATtacggcgtatac                                            >  1:361414/1‑148 (MQ=255)
                                                                                                    gTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccgg                                          <  2:767790/150‑1 (MQ=255)
                                                                                                    gTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccgg                                          >  1:546167/1‑150 (MQ=255)
                                                                                                    gTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccgg                                          <  1:68289/150‑1 (MQ=255)
                                                                                                     ttACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccggc                                         <  2:369620/150‑1 (MQ=255)
                                                                                                      tACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccggcg                                        >  2:77456/1‑150 (MQ=255)
                                                                                                      tACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCAATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccggcg                                        <  1:21571/150‑1 (MQ=255)
                                                                                                        cGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccggcgta                                      >  1:1367437/1‑150 (MQ=255)
                                                                                                          tATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGa                                    >  2:1664922/1‑150 (MQ=255)
                                                                                                           aTCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCCCCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGAt                                   >  1:786650/1‑150 (MQ=255)
                                                                                                           aTCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGAt                                   <  1:640274/150‑1 (MQ=255)
                                                                                                           aTCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGAt                                   >  2:476390/1‑150 (MQ=255)
                                                                                                           aTCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGAt                                   >  1:1693592/1‑150 (MQ=255)
                                                                                                             cTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTGCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTc                                 >  2:1565731/1‑150 (MQ=255)
                                                                                                             cTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTc                                 >  2:896794/1‑150 (MQ=255)
                                                                                                                gAAATAAGGTTGAAAAATTAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCATGCGTATCCTGCGTAGATTCGaa                              <  2:241534/150‑1 (MQ=255)
                                                                                                                 aaaTAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGaaa                             >  2:876949/1‑150 (MQ=255)
                                                                                                                   aTAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAAtt                           >  2:1687743/1‑149 (MQ=255)
                                                                                                                      tggTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGttt                        <  1:756113/149‑1 (MQ=255)
                                                                                                                      aGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTACAACGTGAATGGTGGAGATGATTTCCGGCGTATCCGGGGGAGATGCGAAATGtt                         >  2:213279/1‑150 (MQ=255)
                                                                                                                        tttGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAc                      >  1:31850/2‑150 (MQ=255)
                                                                                                                         ttGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACg                     <  1:1664922/150‑1 (MQ=255)
                                                                                                                         tagAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTACTCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACg                     <  1:115465/148‑1 (MQ=255)
                                                                                                                          tGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGa                    >  2:336571/1‑150 (MQ=255)
                                                                                                                              aaaTAAACACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAc                <  1:668050/150‑1 (MQ=255)
                                                                                                                               aaaaaaaaCGGCGCTTTTTAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcg               >  1:984350/4‑150 (MQ=255)
                                                                                                                               aaTAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcg               <  2:1367437/150‑1 (MQ=255)
                                                                                                                               aaTAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcg               <  2:1416730/150‑1 (MQ=255)
                                                                                                                                aaaaaaaCGGCGCTTTTTAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcgg              >  2:1002250/3‑150 (MQ=255)
                                                                                                                                 tAAAAACGGCGCTAAAAAGCGCGG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTAGAAATGTTTACGAATAcggc             >  2:1216013/1‑150 (MQ=255)
                                                                                                                                 tAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcggc             <  2:1583621/150‑1 (MQ=255)
                                                                                                                                 aaaaaaCGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcggc             <  2:1302607/149‑1 (MQ=255)
                                                                                                                                   aaaaCGGCGCTTTTTAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcggcgg           >  2:1569905/1‑150 (MQ=255)
                                                                                                                                    aaaCGGCGCTTTTTAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGa          >  1:1100512/1‑150 (MQ=255)
                                                                                                                                    aaaCGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGa          >  1:530490/1‑150 (MQ=255)
                                                                                                                                     aaCGGCGCTTTTTAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGAt         >  2:657310/1‑150 (MQ=255)
                                                                                                                                     aaCGGCGCTTTTTAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGAt         >  2:799541/1‑150 (MQ=255)
                                                                                                                                     aaCGGCGCTTTTTAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCAACGTGAATGGTGGCGAGGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGAt         >  2:362907/1‑150 (MQ=255)
                                                                                                                                       cGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCg       <  2:1430832/150‑1 (MQ=255)
                                                                                                                                       cGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGCGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCg       >  1:794962/1‑150 (MQ=255)
                                                                                                                                       cGGCGCTAAAAAGCGACG‑TTTTTTTTGACTGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCg       <  2:1424139/150‑1 (MQ=255)
                                                                                                                                          cgcTAAAAAGCGCGG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCa    <  1:1546709/150‑1 (MQ=255)
                                                                                                                                          cgcTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACAAATACGGCGGATCGTCa    <  2:505095/150‑1 (MQ=255)
                                                                                                                                            cTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACg  >  2:344445/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACG  >  NC_000913/4639412‑4639701

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