Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,995,364 A→T 6.9% intergenic (‑5/‑454) uvrY ← / → yecF response regulator in two‑component regulatory system with BarA/DUF2594 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,995,3640AT 6.9% 291.3 / 10.8 88intergenic (‑5/‑454)uvrY/yecFresponse regulator in two‑component regulatory system with BarA/DUF2594 family protein
Reads supporting (aligned to +/- strand):  ref base A (37/43);  new base T (3/3);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

ATTTGTCCGGCACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTGTCCACAGTTATGGTCACGCC  >  NC_000913/1995222‑1995511
                                                                                                                                              |                                                                                                                                                   
aTTTGTCCGGCACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTcc                                                                                                                                              >  2:1548046/1‑150 (MQ=255)
aTTTGTCCGGCACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTcc                                                                                                                                              >  1:122939/1‑150 (MQ=255)
aTTTGTCCGGCACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTcc                                                                                                                                              >  1:122957/1‑150 (MQ=255)
aTTTGTCCGGCACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTcc                                                                                                                                              >  2:332085/1‑150 (MQ=255)
aTTTGTCCGGCACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTATTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTCGTCATCAACAAGTAGAACGTTGATCAAAGGACTATCTcc                                                                                                                                              >  1:1316778/1‑150 (MQ=255)
        ggCACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAg                                                                                                                                      >  1:1415330/1‑150 (MQ=255)
            ccacTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGAt                                                                                                                                  >  2:1519455/1‑150 (MQ=255)
            ccacTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGAt                                                                                                                                  >  2:264190/1‑150 (MQ=255)
              acTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATaa                                                                                                                                >  1:1457282/1‑150 (MQ=255)
               cTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAAc                                                                                                                               <  2:217771/150‑1 (MQ=255)
                   aCGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTa                                                                                                                           <  1:831869/150‑1 (MQ=255)
                   aCGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTa                                                                                                                           <  1:1397600/150‑1 (MQ=255)
                       cGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGa                                                                                                                       <  2:1681502/150‑1 (MQ=255)
                       cGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGa                                                                                                                       <  1:1299953/150‑1 (MQ=255)
                        gTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGAt                                                                                                                      <  1:1370968/150‑1 (MQ=255)
                         tCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGAta                                                                                                                     <  1:665147/150‑1 (MQ=255)
                             cACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGAtagttg                                                                                                                 <  1:1693019/150‑2 (MQ=255)
                              aCCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGAtagttag                                                                                                                >  1:578240/1‑150 (MQ=255)
                              aCCGCACGATGCCTCACCGACGACCTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGAtagttag                                                                                                                >  1:578278/1‑150 (MQ=255)
                               ccGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATTCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGAtagttagt                                                                                                               >  2:1696266/1‑150 (MQ=255)
                                cGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTc                                                                                                              >  1:1372001/1‑150 (MQ=255)
                                 gCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCa                                                                                                             >  1:394053/1‑150 (MQ=255)
                                     gATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCg                                                                                                         >  1:1628564/1‑150 (MQ=255)
                                      aTGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAAGAGGGATAACGCGACTGATAGTTAGTCATTCGt                                                                                                        <  2:355532/150‑1 (MQ=255)
                                         ccTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGtata                                                                                                     >  1:204800/1‑150 (MQ=255)
                                          cTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTatat                                                                                                    >  2:164447/1‑150 (MQ=255)
                                          cTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTAACCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTatat                                                                                                    >  2:817761/1‑150 (MQ=255)
                                           tCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATAtt                                                                                                   <  1:348667/150‑1 (MQ=255)
                                           tCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATAtt                                                                                                   <  1:348887/150‑1 (MQ=255)
                                              ccGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGc                                                                                                <  1:908217/150‑1 (MQ=255)
                                                    acTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGtttt                                                                                          <  2:1625995/150‑1 (MQ=255)
                                                    acTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGtttt                                                                                          <  2:135516/150‑1 (MQ=255)
                                                      ttttATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTa                                                                                        >  1:699821/1‑150 (MQ=255)
                                                      ttttATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTa                                                                                        <  2:1114659/150‑1 (MQ=255)
                                                      ttttATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTa                                                                                        >  1:487564/1‑150 (MQ=255)
                                                      ttttATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTa                                                                                        <  2:1424937/150‑1 (MQ=255)
                                                      ttttATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTa                                                                                        <  2:224873/150‑1 (MQ=255)
                                                       tttATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTAt                                                                                       <  2:459537/150‑1 (MQ=255)
                                                       tttATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTAt                                                                                       <  1:1042498/150‑1 (MQ=255)
                                                           taCCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCa                                                                                   >  2:946225/1‑150 (MQ=255)
                                                               cTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGATTATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACtt                                                                               >  2:696817/1‑150 (MQ=255)
                                                               cTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGATTATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACtt                                                                               >  2:696795/1‑150 (MQ=255)
                                                               cTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGCGGTCATCAACAAGTAGAACGTTGATCAAAGGATTATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACtt                                                                               >  2:624483/1‑150 (MQ=255)
                                                                    tatCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGaaa                                                                          <  1:312084/150‑1 (MQ=255)
                                                                    tatCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGaaa                                                                          <  1:312126/150‑1 (MQ=255)
                                                                        ttCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGtt                                                                      <  2:863458/150‑1 (MQ=255)
                                                                        ttCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACCtt                                                                      <  2:1195174/150‑1 (MQ=255)
                                                                        ttCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCGCCAGAAATAGGGATAACGCTACTGTTAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGtt                                                                      <  2:926339/150‑1 (MQ=255)
                                                                           cAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTaaa                                                                   <  2:943112/150‑1 (MQ=255)
                                                                           cAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTaaa                                                                   <  2:659035/150‑1 (MQ=255)
                                                                                tGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAAtgg                                                              <  1:1395959/150‑1 (MQ=255)
                                                                                  cgtcgtATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAAtggtg                                                            <  1:534056/150‑1 (MQ=255)
                                                                                             tGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGAtt                                                 >  1:115961/1‑150 (MQ=255)
                                                                                             tGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGAtt                                                 >  1:341484/1‑150 (MQ=255)
                                                                                             tGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGAtt                                                 <  2:710120/150‑1 (MQ=255)
                                                                                               cgcgcACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGAtttt                                               >  1:1545091/1‑150 (MQ=255)
                                                                                                gcgcACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGAttttt                                              <  2:1703199/150‑1 (MQ=255)
                                                                                                   cACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAGTATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTtctc                                           >  2:1379920/1‑150 (MQ=255)
                                                                                                   cACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTtctc                                           >  2:1379914/1‑150 (MQ=255)
                                                                                                   cACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTtctc                                           <  1:946225/150‑1 (MQ=255)
                                                                                                   cACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTGTtctc                                           >  2:1379928/1‑150 (MQ=255)
                                                                                                      cAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGCTAAAAATGGTGGGTATTCGATTTTTCTCTCg                                        <  1:1590675/150‑1 (MQ=255)
                                                                                                         ttCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGttt                                     >  1:331343/1‑150 (MQ=255)
                                                                                                             tGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGt                                 >  2:470927/1‑150 (MQ=255)
                                                                                                             tGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGAATATTAGCGGTTGTTATTGCAACTTTGAAACGTAAAAAATGGTGTGTATTCGATATTTCTCTCGATTTTGt                                 >  1:1137523/1‑150 (MQ=255)
                                                                                                              ggTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTc                                <  1:1637758/150‑1 (MQ=255)
                                                                                                              ggTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTc                                <  2:1420582/150‑1 (MQ=255)
                                                                                                              ggTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTc                                <  2:444977/150‑1 (MQ=255)
                                                                                                              ggTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTc                                <  2:775786/150‑1 (MQ=255)
                                                                                                               gTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCa                               <  2:681148/150‑1 (MQ=255)
                                                                                                                   tcaACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTc                           >  2:769348/1‑150 (MQ=255)
                                                                                                                   acaacaaGTAGAACGTTGATCAAACGAACATCACCATAAATACCGATAACGCTACTGATAGTTAGACATTCGTATATTAGCGGTTTTTATTGCAACTTAGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTc                           <  2:1063649/149‑1 (MQ=255)
                                                                                                                   acaacaaGTACAACGTTGATCAAACGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGACATTCGTATATTAGCGCTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTc                           <  2:1063650/149‑1 (MQ=255)
                                                                                                                     aacaaGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGa                         >  1:418302/1‑150 (MQ=255)
                                                                                                                      acaaGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGAt                        >  1:899020/1‑150 (MQ=255)
                                                                                                                       caaGTAGAATGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGAtt                       <  2:1276291/150‑1 (MQ=255)
                                                                                                                       caaGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGAtt                       <  1:1696266/150‑1 (MQ=255)
                                                                                                                        aaGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTg                      >  1:404994/1‑150 (MQ=255)
                                                                                                                          gTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTGTc                    <  1:548592/150‑1 (MQ=255)
                                                                                                                           tAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTGTcc                   >  2:659326/1‑150 (MQ=255)
                                                                                                                           tAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTGTcc                   >  2:1366229/1‑150 (MQ=255)
                                                                                                                                   tGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTGTCCACAGTTAt           <  1:289069/150‑1 (MQ=255)
                                                                                                                                    gATCAAAGGATTATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTGTCCACAGTTATg          <  1:696817/150‑1 (MQ=255)
                                                                                                                                    gATCAAAGGATTATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTGTCCACAGTTATg          <  1:624483/150‑1 (MQ=255)
                                                                                                                                    gATCAAAGGATTATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTGTCCACAGTTATg          <  1:696795/150‑1 (MQ=255)
                                                                                                                                     aTCAAACGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTCGTCATTCGGATATGAGCGGTTTTTAGTGCAACATTGAAACGTTAAAAATCGTGGGTATTCGATTTTTCTCTCGATTTTGTCAATTCGATTGTCCACAGAGATgg         >  2:1334066/1‑150 (MQ=255)
                                                                                                                                        aaaGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTGTCCACAGTTATGGTCa      >  1:1390343/1‑150 (MQ=255)
                                                                                                                                            gAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTGTCCACAGTTATGGTCACGcc  >  1:611141/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
ATTTGTCCGGCACCACTTAACGGCGTCTTCACCGCACGATGCCTCACCGACGACTTTTATACCCTTTATATCTTCCAGAATGCGTCGTATCCCTGCGCGCACCAGTTCGTGGTCATCAACAAGTAGAACGTTGATCAAAGGAATATCTCCAGAAATAGGGATAACGCTACTGATAGTTAGTCATTCGTATATTAGCGGTTTTTATTGCAACTTTGAAACGTTAAAAATGGTGGGTATTCGATTTTTCTCTCGTTTTTGTCAATTCGATTGTCCACAGTTATGGTCACGCC  >  NC_000913/1995222‑1995511

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