Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,639,557 A→T 5.6% intergenic (+27/+33) 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,5570AT 5.6% 380.6 / 8.0 133intergenic (+27/+33)creD/arcAinner membrane protein/response regulator in two‑component regulatory system with ArcB or CpxA
Reads supporting (aligned to +/- strand):  ref base A (75/50);  new base T (2/6);  total (77/56)
Fisher's exact test for biased strand distribution p-value = 6.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 141 total aligned reads displayed.

GTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTAC  >  NC_000913/4639408‑4639706
                                                                                                                                                     |                                                                                                                                                      
gtggtggCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTaaaa                                                                                                                                                        >  1:886108/1‑150 (MQ=255)
        gcTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGt                                                                                                                                               >  2:589937/1‑150 (MQ=255)
         cTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGtt                                                                                                                                              >  1:1580835/1‑150 (MQ=255)
           ggCCGGCATGATGTTTGTGACCCGTAATATCGACTGTTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGtttt                                                                                                                                            <  2:1187566/150‑1 (MQ=255)
               ggCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTTTTTAGCGCCGtttttatt                                                                                                                                        <  2:140615/150‑1 (MQ=255)
               ggCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGtttttttt                                                                                                                                        <  2:1098845/150‑1 (MQ=255)
                gCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTTCCGAAAATGAAAGCCAGTAAATAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTg                                                                                                                                       <  2:1299289/150‑1 (MQ=255)
                   tgatgTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTTTTTAGCGCCGtttttatttttc                                                                                                                                    <  2:1284790/150‑5 (MQ=255)
                    gatgTTTGTGACCCGTAATATCGACTGGTATGCGTGTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTGACAACGGCCGAGGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAACCGGCGCTAAAAAGCGCCGTTTTTTTTGACgg                                                                                                                                   >  1:1243457/1‑150 (MQ=255)
                       gTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACggtgg                                                                                                                                >  2:1668068/1‑150 (MQ=255)
                           gtgACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTaaa                                                                                                                            >  2:1772031/1‑150 (MQ=255)
                             gACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGc                                                                                                                          >  2:497456/1‑150 (MQ=255)
                               cccGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCg                                                                                                                        >  1:597264/1‑150 (MQ=255)
                               cccGTAATATCGACTGGTATGCGTTTCCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCg                                                                                                                        >  2:675043/1‑150 (MQ=255)
                                 cGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAgactat                                                                                                                      >  2:1132039/1‑145 (MQ=255)
                                  gTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGAtt                                                                                                                     >  2:1765118/1‑150 (MQ=255)
                                  gTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGAtt                                                                                                                     >  2:752304/1‑150 (MQ=255)
                                    aaTATCGAGTTGTATGCGTTTTCACTGCCGAAATTGAAAGCAAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTTACGGTGGTAAAGCCGATTaa                                                                                                                   <  1:308811/150‑1 (MQ=255)
                                    aaTATCGACTGGTATGCGTATTCACTGCCGAAACTGAAAGCCAGTAAAGAAGTTACACCGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAAAAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTaa                                                                                                                   <  1:1621175/150‑1 (MQ=255)
                                          gACAGGTTTGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTcc                                                                                                             <  2:973432/150‑1 (MQ=255)
                                            cTGGTATGCGTTTTCAGTGAGGAAAATGAAAGCCTGTAAATAATTAAAAAGGGACGATGATTTACTTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAATAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAg                                                                                                           <  1:931819/150‑1 (MQ=255)
                                                   gCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCg                                                                                                    >  1:17403/1‑150 (MQ=255)
                                                    cGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGc                                                                                                   <  2:6001/150‑1 (MQ=255)
                                                       tttCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGa                                                                                                >  1:212311/1‑150 (MQ=255)
                                                       tttCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGa                                                                                                >  1:255404/1‑150 (MQ=255)
                                                       tttCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGa                                                                                                >  1:415213/1‑150 (MQ=255)
                                                         tCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAg                                                                                              >  2:459435/1‑150 (MQ=255)
                                                          cACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGc                                                                                             >  2:501352/1‑150 (MQ=255)
                                                                cGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAAc                                                                                       >  2:690482/1‑150 (MQ=255)
                                                                  aaaaTGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGAt                                                                                                         >  2:16056/1‑130 (MQ=255)
                                                                   aaaTGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCtt                                                                                    >  2:1513449/1‑150 (MQ=255)
                                                                     aTGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCa                                                                                  >  2:1527448/1‑150 (MQ=255)
                                                                     aTGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCa                                                                                  >  1:1304417/1‑150 (MQ=255)
                                                                       gtaaGCCCGTAAAGAATTAACTACGGATGAGTAGTTACGTATCTGGAAATAAGTTGGAAAAAAAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACCGTGGTAAAGCCGATTAATCTTGCAGATCACCGCAGAAGCGATAAGCTTCAcc                                                                                <  2:1118255/148‑1 (MQ=255)
                                                                       gAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCAcc                                                                                >  1:1616556/1‑150 (MQ=255)
                                                                        aaaGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCg                                                                               >  2:1840532/1‑150 (MQ=255)
                                                                        aaaGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCg                                                                               >  2:1808745/1‑150 (MQ=255)
                                                                        aaaGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGGGGTAAAGCCGATTAATCGTCCAGATCACCGCAGAAGCGATTACCTTCACCg                                                                               >  2:378030/1‑150 (MQ=255)
                                                                         aaGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGCCGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGt                                                                              >  1:507982/1‑150 (MQ=255)
                                                                           gCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGa                                                                            >  1:727879/1‑150 (MQ=255)
                                                                            ccAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGaa                                                                           >  2:1547704/1‑150 (MQ=255)
                                                                            ccAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGaa                                                                           >  2:1001042/1‑150 (MQ=255)
                                                                             cAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAAt                                                                          >  1:1784549/1‑150 (MQ=255)
                                                                             cAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAAt                                                                          >  1:582355/1‑150 (MQ=255)
                                                                               gTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAAACTTCACCGTGAAtg                                                                         >  1:913049/1‑150 (MQ=255)
                                                                               gTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAAtgg                                                                        <  1:443489/150‑1 (MQ=255)
                                                                                  aagaagTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGc                                                                    <  1:722869/150‑1 (MQ=255)
                                                                                   cgacgTTACAACGGACGATGCGTTACGTATCTGGAAATAAGGTTGAAAAAAAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGc                                                                    <  2:1530921/146‑1 (MQ=255)
                                                                                   agcagTTACAACGGACGATGAGGTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGc                                                                    <  2:670286/147‑1 (MQ=255)
                                                                                   agaagTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGc                                                                    <  2:1223860/150‑1 (MQ=255)
                                                                                   agaagTTACAACGGACGAAGAGTTACGTATCTGTAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGc                                                                    <  1:1748145/150‑1 (MQ=255)
                                                                                    gaagTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGTAGCGATAACCTTCACCGTGAATGGTGGCg                                                                   <  1:497456/150‑1 (MQ=255)
                                                                                    gaagTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCg                                                                   >  1:1057298/1‑150 (MQ=255)
                                                                                      atttACACCGTACGATGATTTAAGTATCTGTAAATAGGGTGTAAAAAAAAAAACTGCGCTATAACGAGCCGTTTTTTTTGACGGTGGTAAATCCGAGTAATCTTCCAGATCACCTCAGAATCGATAACCTTCACCGTGAATGGTGGCgat                                                                 <  1:562929/148‑1 (MQ=255)
                                                                                       gTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAATCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCgatg                                                                >  1:1446193/1‑150 (MQ=255)
                                                                                       gTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATATTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCgatt                                                                >  1:1632208/1‑149 (MQ=255)
                                                                                        ttACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCgatga                                                               >  2:139512/1‑150 (MQ=255)
                                                                                         tACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCgatgat                                                              <  2:259068/150‑1 (MQ=255)
                                                                                           cAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGAttt                                                            >  2:12068/1‑150 (MQ=255)
                                                                                           cAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGAttt                                                            >  1:1536093/1‑150 (MQ=255)
                                                                                           cAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGGAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGAttt                                                            >  2:1750138/1‑150 (MQ=255)
                                                                                           cAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGATAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCAGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGAttt                                                            >  2:1757171/1‑150 (MQ=255)
                                                                                             aCGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtcc                                                          >  1:374605/1‑150 (MQ=255)
                                                                                                gACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAACACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGCAAAGCCGATTAATCTTCCCGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggc                                                       <  2:459367/150‑1 (MQ=255)
                                                                                                 accaTGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATAtccggcg                                                      <  2:598161/147‑1 (MQ=255)
                                                                                                      taGTTACGTATCTGGAAATAAGTTTTAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatcc                                                 <  1:1411464/149‑1 (MQ=255)
                                                                                                       aGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccg                                                <  1:1808745/150‑1 (MQ=255)
                                                                                                         tttcGTATCTGGAAATAAGGTTGAAAAATAATAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccggc                                              <  2:972759/147‑1 (MQ=255)
                                                                                                            cGTATCTGGAAATATGGTTGATAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTTATCTTCCAGATAACCGCAGAAGCGATAACGTTCACCGTGTATGGTGGCGTTGATTtcgggggtttccggcgtt                                           >  2:1181970/1‑149 (MQ=255)
                                                                                                             gTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAg                                          <  2:1447865/150‑1 (MQ=255)
                                                                                                                tCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACTGCAGAAGCTATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGAtt                                       >  2:170395/1‑150 (MQ=255)
                                                                                                                tCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGAtt                                       >  1:961274/1‑150 (MQ=255)
                                                                                                                  tGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCg                                     <  1:1838064/150‑1 (MQ=255)
                                                                                                                   ggAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGa                                    >  1:1807213/1‑150 (MQ=255)
                                                                                                                   ggAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGa                                    >  1:115880/1‑150 (MQ=255)
                                                                                                                   ggAAATAAGGTAGAAAAATAAAAACGGCGCTAAACGCGCCGTTTTTTTTGACGGAGGTAAAGCCGATTAATCATCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGaa                                   <  1:905345/150‑1 (MQ=255)
                                                                                                                    gAAATAAGGTTGAAAAATAAACACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGCAGATTCGaa                                   <  2:753998/150‑1 (MQ=255)
                                                                                                                    gAAATAAGGTTGAAAAATAAAAACGGCGCTAAAACGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACAGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCAGGCGTAGATTCGaa                                   >  2:119712/1‑150 (MQ=255)
                                                                                                                           ggTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTa                            <  1:745477/150‑1 (MQ=255)
                                                                                                                            gTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAc                           <  2:255404/150‑1 (MQ=255)
                                                                                                                            gTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAc                           >  2:637647/1‑150 (MQ=255)
                                                                                                                            gTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCAAAATGTTTAc                           <  2:845241/150‑1 (MQ=255)
                                                                                                                             ttGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACg                          <  2:841529/150‑1 (MQ=255)
                                                                                                                             ttGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACg                          <  2:585332/150‑1 (MQ=255)
                                                                                                                              tGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGa                         >  1:1602568/1‑150 (MQ=255)
                                                                                                                                     tAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcggc                  <  2:219666/150‑1 (MQ=255)
                                                                                                                                      aaaaaCGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcggcg                 <  1:1550019/150‑1 (MQ=255)
                                                                                                                                      aaaaaCGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcggcg                 <  1:1205161/150‑1 (MQ=255)
                                                                                                                                       aaaaAGGAGCTAAAAAGAGCCGTTTTTTTAGACGGTGGTAAAGCCGATTAATCTTCCAGATCAACGATGAAGCGATAACCTTCAACGTGAATGGTGGCGATGATTTACGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATActgcgg                >  2:481450/1‑150 (MQ=255)
                                                                                                                                           cGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCg            >  1:263756/1‑150 (MQ=255)
                                                                                                                                           cGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCg            <  1:876154/150‑1 (MQ=255)
                                                                                                                                           cGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCg            <  1:1772031/150‑1 (MQ=255)
                                                                                                                                            ggCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGt           <  1:765629/150‑1 (MQ=255)
                                                                                                                                                cTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACg       >  1:1090347/1‑150 (MQ=255)
                                                                                                                                                  aaaaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTc     >  2:1502109/1‑150 (MQ=255)
                                                                                                                                                  aaaaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCAAAGTGAATGGTGGCGATGATTTCCAGAGTATCCTGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGta     >  1:202722/1‑149 (MQ=255)
                                                                                                                                                   aaaaGCGCCGTTTTTTTTGGCGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCt    >  1:1323269/1‑150 (MQ=255)
                                                                                                                                                   aaaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCt    <  2:374605/150‑1 (MQ=255)
                                                                                                                                                   aaaaGCGCCGTTTTTTTTGACGGTGCTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCt    >  1:1323264/1‑150 (MQ=255)
                                                                                                                                                     aaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTAc  <  1:1513449/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTAC  >  NC_000913/4639408‑4639706

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: