Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,983 A→T 5.8% intergenic (+46/+30) matP → / ← ompA Ter macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,9830AT 5.8% 1087.9 / 39.0 306intergenic (+46/+30)matP/ompATer macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base A (131/155);  new base T (12/6);  total (145/161)
Fisher's exact test for biased strand distribution p-value = 9.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 321 total aligned reads displayed.

AAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGC  >  NC_000913/1018834‑1019131
                                                                                                                                                     |                                                                                                                                                    
        gTTCAGCTTATTGATGATGCTTAAAACAAAGATAAATACGCTAATAAATTTTCTTCTCTGAAGCAGGATCTGCAGGCTTTGGTGGTTAAGTAATAACTTACTAATTTCATTTAAATTTAAAAACGGAAAAAAAAACCCCGCAGCAGCggg                                                                                                                                              <  1:1846519/150‑1 (MQ=255)
        gTTCAGCTGATTGACGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGTGTATGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCggg                                                                                                                                              <  1:2893039/150‑1 (MQ=255)
        gTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCggg                                                                                                                                              <  2:1564086/150‑1 (MQ=255)
        gTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCggg                                                                                                                                              <  2:2768051/150‑1 (MQ=255)
           cAGCTGATTGAAGATGCGCAAAACAAAGATAAATACGCGAATAAACTGTCTTCTCTGAAGCAGGATCTGCCGGCTTTGCTGGTTAAGGAATAACTGACTAAAGACAGTACAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                           <  1:3206628/150‑1 (MQ=255)
              cTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGttttt                                                                                                                                        >  1:1547768/1‑150 (MQ=255)
                gATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCt                                                                                                                                      <  2:346627/150‑1 (MQ=255)
                gATGGAAGCTGCGGAAAACAAAGAGAAAGACGCGAATAAAATGTCTTCTCTGACGCCGGAGGTGCAGGCTTTGCTGGGTAAGGAATACCTGACGAAAGTCTGTTCAATTTACTAACGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCt                                                                                                                                      <  1:985243/150‑1 (MQ=255)
                 aTTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCTGGCATTGCTGGGTAAGGAATAACTGACGTAAGTCAGTTCAATTTACTAAAGGCAATATAATCCCCGCAGCTGCGGGGTTTTTGTa                                                                                                                                     >  2:4173036/1‑150 (MQ=255)
                      aGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGa                                                                                                                                <  2:3996388/150‑1 (MQ=255)
                      aGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGa                                                                                                                                <  2:1328538/150‑1 (MQ=255)
                         tGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACga                                                                                                                             <  2:1228083/150‑1 (MQ=255)
                            ggAAAACAAAGATAAATACGCGAATAAAATGTCTTCTCTGAAGCAGTATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                          <  1:673119/150‑1 (MQ=255)
                            ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCCGCAGCGGGGTGTTTCTACCCGACgagac                                                                                                                          >  1:3441319/1‑149 (MQ=255)
                            ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                          >  1:3650495/1‑150 (MQ=255)
                            ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAAGTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                          >  2:3078361/1‑150 (MQ=255)
                            ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGATGGGTAAGTAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                          >  2:1071204/1‑150 (MQ=255)
                             gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAAc                                                                                                                         >  2:1915724/1‑150 (MQ=255)
                                 aCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                     <  2:2503267/150‑1 (MQ=255)
                                    aaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGcc                                                                                                                  >  2:4142834/1‑150 (MQ=255)
                                     agagAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCt                                                                                                                 <  2:173034/150‑1 (MQ=255)
                                       agaAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGc                                                                                                               >  1:783781/1‑150 (MQ=255)
                                       agaAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGc                                                                                                               >  2:70360/1‑150 (MQ=255)
                                               gcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGt                                                                                                       >  1:1303624/1‑150 (MQ=255)
                                                  aaTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATACCTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTAc                                                                                                    <  2:178794/150‑1 (MQ=255)
                                                   aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACa                                                                                                   >  2:54918/1‑150 (MQ=255)
                                                       aaTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGt                                                                                               >  2:1597821/1‑150 (MQ=255)
                                                        aTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGa                                                                                                         <  1:2450662/139‑1 (MQ=255)
                                                        aTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTc                                                                                              <  1:3923016/150‑1 (MQ=255)
                                                          gTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCtt                                                                                            <  1:240257/150‑1 (MQ=255)
                                                           tcttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                           >  1:1567338/1‑150 (MQ=255)
                                                             ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGACCTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                         >  1:3818003/1‑150 (MQ=255)
                                                             ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                         >  1:1452588/1‑150 (MQ=255)
                                                             ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                         >  2:26803/1‑150 (MQ=255)
                                                             ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                         >  2:2364138/1‑150 (MQ=255)
                                                             ttctCTGAAGCAGGATCTGCAGGCAGTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                         <  2:2901225/150‑1 (MQ=255)
                                                              tctcTGAAGCCGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGAt                                                                                        <  1:3998511/150‑1 (MQ=255)
                                                              tctcTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGAt                                                                                        <  2:725176/150‑1 (MQ=255)
                                                                tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAc                                                                                      >  1:3686310/1‑150 (MQ=255)
                                                                  tGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAATTCATTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCt                                                                                    <  1:2654852/150‑1 (MQ=255)
                                                                      gCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTaa                                                                                <  1:1832227/150‑1 (MQ=255)
                                                                      gCAGGATCTGCAGGCAGTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTaa                                                                                <  1:1145256/150‑1 (MQ=255)
                                                                          gATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAACAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTCAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTc                                                                            <  2:603672/150‑1 (MQ=255)
                                                                          gATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTc                                                                            >  1:3753718/1‑150 (MQ=255)
                                                                          gATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTc                                                                            >  1:3552472/1‑150 (MQ=255)
                                                                                 aGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATctct                                                                     >  1:1372764/1‑150 (MQ=255)
                                                                                  ggCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAc                                                                   <  1:705563/150‑1 (MQ=255)
                                                                                       tGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcga                                                               <  1:3863553/150‑1 (MQ=255)
                                                                                        gCTGGGTAAGTAATAACTGACTAAAGTCAGTTAAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcgac                                                              <  1:3340342/150‑1 (MQ=255)
                                                                                                aGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGa                                                      <  1:4342670/150‑1 (MQ=255)
                                                                                                  gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc                                                    >  2:1444598/1‑150 (MQ=255)
                                                                                                     taccTGACGAACGTCAGTTCAATTTACTAACGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAg                                                 <  2:2314105/147‑1 (MQ=255)
                                                                                                     tAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAg                                                 <  1:1597821/150‑1 (MQ=255)
                                                                                                     tAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAg                                                 <  2:1141397/150‑1 (MQ=255)
                                                                                                      aaCTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAgg                                                <  2:470337/150‑1 (MQ=255)
                                                                                                        cTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCa                                              >  2:3824019/1‑150 (MQ=255)
                                                                                                        cTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCa                                              <  2:3896503/150‑1 (MQ=255)
                                                                                                        cTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCa                                              <  2:2465134/150‑1 (MQ=255)
                                                                                                        cTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCa                                              <  2:1985937/150‑1 (MQ=255)
                                                                                                         tGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAg                                             >  1:3259230/1‑150 (MQ=255)
                                                                                                         tGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAg                                             >  1:3017092/1‑150 (MQ=255)
                                                                                                          gACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGt                                            >  2:2838069/1‑150 (MQ=255)
                                                                                                           aCGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTc                                           <  1:1922501/150‑1 (MQ=255)
                                                                                                           aCGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTc                                           <  2:3636378/150‑1 (MQ=255)
                                                                                                              aaaGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAt                                        >  1:2783427/1‑150 (MQ=255)
                                                                                                                aGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGTGGTTTTTCTACCAGACGAGAACTTTAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCa                                      <  1:814574/150‑1 (MQ=255)
                                                                                                                aGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCa                                      <  1:929096/150‑1 (MQ=255)
                                                                                                                  tCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGt                                    >  1:3979840/1‑150 (MQ=255)
                                                                                                                   catttCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTg                                   <  1:2195639/147‑1 (MQ=255)
                                                                                                                   cAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTg                                   >  1:3490606/1‑150 (MQ=255)
                                                                                                                       tCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcagc                               <  1:1958692/150‑1 (MQ=255)
                                                                                                                           tttACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGc                           >  1:3713177/1‑150 (MQ=255)
                                                                                                                             tACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTg                         <  2:2546351/150‑1 (MQ=255)
                                                                                                                               cTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGtt                       <  1:2907481/150‑1 (MQ=255)
                                                                                                                                tAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt                      >  1:393347/1‑150 (MQ=255)
                                                                                                                                tAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTGTGTCTACCAGACGAGAACTTAAGCCGGCGGCTGAGTTACAACGGCTTTGAGACCTGTAACTGCGATCTCTCCGCGACGATCCGGAGCCAGGCCGTCGAGCCGTGCCGCACGCTgggg                      >  2:2132764/1‑147 (MQ=255)
                                                                                                                                  aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCa                    >  2:310933/1‑150 (MQ=255)
                                                                                                                                  aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCa                    >  2:2178750/1‑150 (MQ=255)
                                                                                                                                  aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCa                    >  1:708484/1‑150 (MQ=255)
                                                                                                                                   aGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCAc                   <  1:2355512/150‑1 (MQ=255)
                                                                                                                                   aGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCAc                   <  2:3235516/150‑1 (MQ=255)
                                                                                                                                    ggCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACg                  <  1:4232498/150‑1 (MQ=255)
                                                                                                                                      cAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGtt                <  1:1685724/150‑1 (MQ=255)
                                                                                                                                        aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt              >  1:3137658/1‑150 (MQ=255)
                                                                                                                                        aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt              >  1:4150277/1‑150 (MQ=255)
                                                                                                                                            aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcact          >  2:2578715/1‑149 (MQ=255)
                                                                                                                                            aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAGCAGTGCAGCACGCTGTTTCACGTTGTcaca          >  2:1257532/1‑150 (MQ=255)
                                                                                                                                            aaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAATCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca          >  2:3776188/1‑150 (MQ=255)
                                                                                                                                            aaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca          <  2:262901/150‑1 (MQ=255)
                                                                                                                                            aaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca          <  1:1890596/150‑1 (MQ=255)
                                                                                                                                            aaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca          <  1:583292/150‑1 (MQ=255)
                                                                                                                                            aaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca          <  1:2531444/150‑1 (MQ=255)
                                                                                                                                            aaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcaca          <  2:3048178/150‑1 (MQ=255)
                                                                                                                                             aaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAg         <  2:4223955/150‑1 (MQ=255)
                                                                                                                                               ccccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgt       <  1:1641729/150‑1 (MQ=255)
                                                                                                                                                cccGCCGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCGGCGGCTGAGTTACACCGGCGTTGATACCGGTAACTTCGATCTCTACGCGACGCTCCGGCGCCAGGCAGTCGCTCAGTGCAGCACGCTGTTTCACGGTGTCACCgggg      >  1:2090824/1‑148 (MQ=255)
                                                                                                                                                 ccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgtgt     >  2:143320/1‑150 (MQ=255)
                                                                                                                                                   gcagcaGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTg   <  2:2271327/150‑1 (MQ=255)
                                                                                                                                                   gcagcaGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTg   <  1:3611252/150‑1 (MQ=255)
                                                                                                                                                    cagcagCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGc  <  2:3346045/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGC  >  NC_000913/1018834‑1019131

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