Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,977 C→A 67.5% intergenic (+40/+36) 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,9770CA67.5% 225.8 / 254.3 216intergenic (+40/+36)matP/ompATer macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base C (32/37);  new base A (74/64);  total (109/107)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-03

Only 100 of 221 total aligned reads displayed.

CTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAG  >  NC_000913/1018831‑1019124
                                                                                                                                                  |                                                                                                                                                   
cTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTAATAAAGGCAAAAAAAAacc                                                                                                                                                   >  3:493717/1‑149 (MQ=255)
  gAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCg                                                                                                                                                 >  4:705826/1‑149 (MQ=255)
  gAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCg                                                                                                                                                 >  1:299661/1‑149 (MQ=255)
  gAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGGTCAATTTACTAAAGGCAAAAAAAACCCCg                                                                                                                                                 >  3:458978/1‑149 (MQ=255)
   aaaCGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCgc                                                                                                                                                >  4:183365/1‑149 (MQ=255)
        aTTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGcagcag                                                                                                                                           >  3:1035881/1‑149 (MQ=255)
            ttCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCggg                                                                                                                                       >  2:236701/1‑149 (MQ=255)
            ttCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCggg                                                                                                                                       <  1:1129809/149‑1 (MQ=255)
             tCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCgggg                                                                                                                                      >  3:1221705/1‑149 (MQ=255)
                gCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGttt                                                                                                                                   <  2:357449/149‑1 (MQ=255)
                gCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGttt                                                                                                                                   >  2:114054/1‑149 (MQ=255)
                   gATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTc                                                                                                                                >  1:1336149/1‑149 (MQ=255)
                   gATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCGCAGCAGCGGGGTTTTTCt                                                                                                                               >  4:492447/1‑149 (MQ=255)
                    aTTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCt                                                                                                                               <  1:482946/149‑1 (MQ=255)
                      tGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTAc                                                                                                                             <  3:1165672/149‑1 (MQ=255)
                      tGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTAc                                                                                                                             <  2:539057/149‑1 (MQ=255)
                             gCGGAAAACAAAGATAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACga                                                                                                                      <  3:719860/149‑1 (MQ=255)
                             gCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACga                                                                                                                      <  2:967170/149‑1 (MQ=255)
                               ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAGAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACgaga                                                                                                                    >  3:1410735/1‑149 (MQ=255)
                               ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGttt                                                                                                                                    >  3:147525/1‑133 (MQ=255)
                               ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCTGCTGCGGGttt                                                                                                                                    <  4:147525/133‑1 (MQ=255)
                               ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACgaga                                                                                                                    >  2:562280/1‑149 (MQ=255)
                               ggAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACgaga                                                                                                                    >  1:375309/1‑149 (MQ=255)
                                gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCTGCTGCGGGtttttt                                                                                                                                 >  4:1074077/1‑135 (MQ=255)
                                gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                   >  3:1147370/1‑149 (MQ=255)
                                gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTGCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGaa                                                                                                                   >  3:1405564/1‑149 (MQ=255)
                                    aCAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                              >  4:1225236/1‑149 (MQ=255)
                                     cAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTACAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                              <  3:1276946/149‑1 (MQ=255)
                                      aaaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAg                                                                                                             >  1:431916/1‑149 (MQ=255)
                                        agagAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAACGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGcc                                                                                                           >  4:907763/1‑149 (MQ=255)
                                                aCGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTg                                                                                                   <  3:940943/149‑1 (MQ=255)
                                                  gcgAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAg                                                                                                 >  3:222075/1‑149 (MQ=255)
                                                    gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTGAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCggg                                                                                                                                       <  2:1022646/109‑1 (MQ=255)
                                                    gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGtt                                                                                               >  2:1390459/1‑149 (MQ=255)
                                                      aTAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAGAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTAc                                                                                             <  4:1410735/149‑1 (MQ=255)
                                                        gaaaTGTCTTCTCTCAAGCAGGATCTGCAGGCACTGCTGGGTAAGGAATCACTGATGCAAGTCAGTTCAGTTTACTAAAGGCAAAAGAAGACCCGCAGCAGCGGGGTTTTTCTACCAGTCGAGAACTTAATCCTGCGGCTGAGTTACaa                                                                                           <  4:1109045/148‑1 (MQ=255)
                                                        aaaaTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAATTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACaa                                                                                           <  1:520711/149‑1 (MQ=255)
                                                              tcttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                    <  2:410282/149‑1 (MQ=255)
                                                               cttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                    <  3:264395/149‑1 (MQ=255)
                                                               cttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                    <  2:1262026/149‑1 (MQ=255)
                                                                ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGACAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTAGAACGTCTTTg                                                                                   >  2:758910/1‑149 (MQ=255)
                                                                ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTg                                                                                   >  4:1443456/1‑149 (MQ=255)
                                                                ttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                    >  1:1413479/1‑148 (MQ=255)
                                                                   tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATa                                                                                <  4:1082249/149‑1 (MQ=255)
                                                                   tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATa                                                                                >  4:32462/1‑149 (MQ=255)
                                                                   tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCgg                                                                                                                                        >  4:688847/1‑93 (MQ=255)
                                                                   tctGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATAc                                                                               >  4:704303/1‑149 (MQ=255)
                                                                        aGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTATTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCttt                                                                           >  3:798573/1‑149 (MQ=255)
                                                                        aGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                  >  1:647576/1‑142 (MQ=255)
                                                                        aGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGa                                                                                  <  2:647576/142‑1 (MQ=255)
                                                                        aGCAGGATCTGCAGGCATTGCTGGGTAAAGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCttt                                                                           <  1:560728/149‑1 (MQ=255)
                                                                             gATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACtt                                                                      >  2:468690/1‑149 (MQ=255)
                                                                                 tGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAt                                                                  <  2:212698/149‑1 (MQ=255)
                                                                                  gCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtc                                                                 >  1:860893/1‑149 (MQ=255)
                                                                                    aGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtctc                                                               >  3:182089/1‑149 (MQ=255)
                                                                                       cATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAc                                                            <  2:519079/149‑1 (MQ=255)
                                                                                         ttGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTAcgc                                                          <  2:373824/149‑1 (MQ=255)
                                                                                            cTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcgac                                                       >  3:538140/1‑149 (MQ=255)
                                                                                             tgagTAAGGAATAACGGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGa                                                                                                                         <  3:999438/79‑1 (MQ=255)
                                                                                             tGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcgacg                                                      >  3:138250/1‑149 (MQ=255)
                                                                                                 tAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATcc                                                  >  4:936930/1‑149 (MQ=255)
                                                                                                    ggAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGa                                               >  4:218921/1‑149 (MQ=255)
                                                                                                     gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAg                                              >  1:1493452/1‑149 (MQ=255)
                                                                                                     gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAg                                              >  1:1267081/1‑149 (MQ=255)
                                                                                                     gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAg                                              >  4:1360704/1‑149 (MQ=255)
                                                                                                       aTAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGcc                                            <  2:446784/149‑1 (MQ=255)
                                                                                                       aTAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGcc                                            >  1:135095/1‑149 (MQ=255)
                                                                                                         aaCTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAg                                          <  1:1183447/149‑1 (MQ=255)
                                                                                                         aaCTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAg                                          <  4:6601/149‑1 (MQ=255)
                                                                                                          aCTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAgg                                         <  2:1054907/149‑1 (MQ=255)
                                                                                                            tGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCt                                                                                                          >  3:1070183/1‑82 (MQ=255)
                                                                                                                gAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACTTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCg                                   >  4:758102/1‑149 (MQ=255)
                                                                                                                gAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCg                                   >  3:534787/1‑149 (MQ=255)
                                                                                                                  aaGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGAt                                 <  3:1424026/149‑1 (MQ=255)
                                                                                                                   aGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATc                                <  1:964756/149‑1 (MQ=255)
                                                                                                                    gTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCa                               <  4:1147370/149‑1 (MQ=255)
                                                                                                                     tCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtc                                                                 >  2:332727/1‑114 (MQ=255)
                                                                                                                     tCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGAtc                                                                 <  1:332727/114‑1 (MQ=255)
                                                                                                                     tCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAg                              >  4:1046600/1‑149 (MQ=255)
                                                                                                                      cAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCGAGGCAGTCGATCAGt                             <  4:866026/149‑1 (MQ=255)
                                                                                                                       aGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTg                            <  4:376292/149‑1 (MQ=255)
                                                                                                                        gTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCtt                                                                            >  1:223966/1‑100 (MQ=255)
                                                                                                                         ttCAATTTACTAAAGGCAAAAAAATACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                          <  1:1003650/149‑1 (MQ=255)
                                                                                                                         ttCAATTTACTAAAGGCAAAAAAAAGCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                          <  3:228942/149‑1 (MQ=255)
                                                                                                                         ttCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAg                                      >  4:1395148/1‑137 (MQ=255)
                                                                                                                         ttCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                          >  2:1329599/1‑149 (MQ=255)
                                                                                                                         ttCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgca                          >  2:726571/1‑149 (MQ=255)
                                                                                                                          tCAATTTACTAAAGGCAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcagc                        <  1:1045285/149‑1 (MQ=255)
                                                                                                                           cAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcagc                        <  1:271524/149‑1 (MQ=255)
                                                                                                                            aaTTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcagca                       >  2:1116941/1‑149 (MQ=255)
                                                                                                                                  cTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGt                 <  4:256814/149‑1 (MQ=255)
                                                                                                                                     aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGa                                               <  3:763302/116‑1 (MQ=255)
                                                                                                                                     aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGa                                               <  3:351385/116‑1 (MQ=255)
                                                                                                                                     aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt               <  4:618342/148‑1 (MQ=255)
                                                                                                                                     aaGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt               >  3:618342/1‑148 (MQ=255)
                                                                                                                                          aaaaaaaaACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGtt         <  1:1460463/149‑1 (MQ=255)
                                                                                                                                           aaaaaaaaCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTg        >  2:910888/1‑149 (MQ=255)
                                                                                                                                           aaaaaaaaCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTg        >  2:917788/1‑149 (MQ=255)
                                                                                                                                               aaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTcac    >  4:1365344/1‑149 (MQ=255)
                                                                                                                                                 aCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAg  <  4:514829/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAG  >  NC_000913/1018831‑1019124

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