Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,197,555 A→G F228S (TTC→TCC)  ymfE ← e14 prophage; uncharacterized protein YmfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,5550AG98.8% 306.8 / ‑4.0 83F228S (TTC→TCC) ymfEe14 prophage; uncharacterized protein YmfE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (47/35);  minor base T (0/1);  total (47/36)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTC  >  NC_000913/1197419‑1197691
                                                                                                                                        |                                                                                                                                        
ctcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGgaaa                                                                                                                                       <  1:1057077/140‑1 (MQ=255)
 tcAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGgaaaa                                                                                                                                      <  1:460044/140‑1 (MQ=255)
    aTGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTc                                                                                                                                   >  1:982703/1‑140 (MQ=255)
      gCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCtt                                                                                                                                 <  1:1286249/140‑1 (MQ=255)
                 aCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCtt                                                                                                                                 >  2:155151/1‑129 (MQ=255)
                 aCTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCtt                                                                                                                                 <  1:155151/129‑1 (MQ=255)
                  cTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGt                                                                                                                     >  1:411755/1‑140 (MQ=255)
                   tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTc                                                                                                                    >  2:173793/1‑140 (MQ=255)
                   tCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTc                                                                                                                    <  2:1106057/140‑1 (MQ=255)
                    cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTcc                                                                                                                   <  2:232362/140‑1 (MQ=255)
                    cGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTcc                                                                                                                   >  1:671597/1‑140 (MQ=255)
                       aaaGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCggg                                                                                                                >  2:102064/1‑140 (MQ=255)
                       aaaGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCggg                                                                                                                >  1:744167/1‑140 (MQ=255)
                       aaaGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCggg                                                                                                                >  1:590243/1‑140 (MQ=255)
                                   agGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgact                                                                                                    <  2:128108/140‑1 (MQ=255)
                                    gggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgactt                                                                                                   >  2:1034918/1‑140 (MQ=255)
                                     ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttg                                                                                                  >  1:405507/1‑140 (MQ=255)
                                     ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttg                                                                                                  >  1:411374/1‑140 (MQ=255)
                                     ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttg                                                                                                  >  1:443575/1‑140 (MQ=255)
                                     ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttg                                                                                                  >  1:1285608/1‑140 (MQ=255)
                                     ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttg                                                                                                  <  1:1034918/140‑1 (MQ=255)
                                     ggCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttg                                                                                                  >  1:976566/1‑140 (MQ=255)
                                       cATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgac                                                                                                >  2:648008/1‑140 (MQ=255)
                                       cATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgacttgac                                                                                                >  2:276313/1‑140 (MQ=255)
                                          cGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACttt                                                                                             <  1:173793/140‑1 (MQ=255)
                                               gAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGtt                                                                                        <  2:67675/140‑1 (MQ=255)
                                                 ggTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCt                                                                                      <  2:411374/140‑1 (MQ=255)
                                                     ttGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgact                                                                                                    <  1:1254015/122‑1 (MQ=255)
                                                     ttGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCgact                                                                                                    >  2:1254015/1‑122 (MQ=255)
                                                      tgctgcATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCt                                                                                 >  1:822204/1‑140 (MQ=255)
                                                      tgctgcATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCt                                                                                 >  1:237781/1‑140 (MQ=255)
                                                        ctgcATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTa                                                                               <  2:691585/140‑1 (MQ=255)
                                                        ctgcATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTa                                                                               <  2:671597/140‑1 (MQ=255)
                                                             ttCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGt                                                                          <  2:405507/140‑1 (MQ=255)
                                                                   gAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagc                                                                    <  1:762035/140‑1 (MQ=255)
                                                                        aaTATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGa                                                               >  2:386123/1‑140 (MQ=255)
                                                                        aaTATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGa                                                               <  2:822204/140‑1 (MQ=255)
                                                                        aaTATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGa                                                               <  2:590243/140‑1 (MQ=255)
                                                                         atatTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGat                                                              >  1:516132/1‑140 (MQ=255)
                                                                          tatTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagc                                                                    >  1:156997/1‑133 (MQ=255)
                                                                          tatTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagc                                                                    <  2:156997/133‑1 (MQ=255)
                                                                              aaCTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATAc                                                         >  2:730100/1‑140 (MQ=255)
                                                                                cTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAATGATAATACca                                                       <  1:648008/140‑1 (MQ=255)
                                                                                cTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACca                                                       <  2:237781/140‑1 (MQ=255)
                                                                                   cGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactc                                                    <  1:276313/140‑1 (MQ=255)
                                                                                   cGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactc                                                    >  2:888212/1‑140 (MQ=255)
                                                                                    gTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactca                                                   >  1:97778/1‑140 (MQ=255)
                                                                                     tGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagca                                                                   <  2:245102/123‑1 (MQ=255)
                                                                                     tGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagca                                                                   >  1:245102/1‑123 (MQ=255)
                                                                                     tGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagc                                                                    >  2:674949/1‑122 (MQ=255)
                                                                                     tGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTcaaagc                                                                    <  1:674949/122‑1 (MQ=255)
                                                                                     tGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcac                                                  >  2:1024574/1‑140 (MQ=255)
                                                                                     tGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcac                                                  >  1:1243378/1‑140 (MQ=255)
                                                                                     tGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcac                                                  >  1:475560/1‑140 (MQ=255)
                                                                                     tGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcac                                                  >  2:786570/1‑140 (MQ=255)
                                                                                     tGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcac                                                  >  2:1051753/1‑140 (MQ=255)
                                                                                     tGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcac                                                  >  2:389857/1‑140 (MQ=255)
                                                                                      gCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACcactcact                                                 <  2:411755/140‑1 (MQ=255)
                                                                                         tattatAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCTCTCACTgag                                              <  2:820351/140‑1 (MQ=255)
                                                                                           ttatAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGaa                                            >  1:74922/1‑140 (MQ=255)
                                                                                                  tctaaGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTTAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATaaa                                     >  2:510654/4‑140 (MQ=255)
                                                                                                    aaaGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAAt                                   >  2:532707/1‑140 (MQ=255)
                                                                                                      aGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTTAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTa                                 <  1:510654/140‑1 (MQ=255)
                                                                                                      aGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTa                                 >  2:927660/1‑140 (MQ=255)
                                                                                                      aGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTa                                 >  2:1067094/1‑140 (MQ=255)
                                                                                                      aGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTa                                 >  2:205827/1‑140 (MQ=255)
                                                                                                      aGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTa                                 >  2:951036/1‑140 (MQ=255)
                                                                                                      aGCAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTa                                 >  1:795840/1‑140 (MQ=255)
                                                                                                        cAAGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACa                               <  2:443575/140‑1 (MQ=255)
                                                                                                          aGTACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaa                             <  1:1067094/140‑1 (MQ=255)
                                                                                                            tACCATTCAATTAATGTTATTTTTAGTGTAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaaaa                           <  1:370577/140‑1 (MQ=255)
                                                                                                            tACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaaaa                           <  2:744167/140‑1 (MQ=255)
                                                                                                            tACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaaaa                           <  2:74922/140‑1 (MQ=255)
                                                                                                            tACCATTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaaaa                           <  2:976566/140‑1 (MQ=255)
                                                                                                                aTTCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGaa                       >  2:103370/1‑140 (MQ=255)
                                                                                                                  tCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCGAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATa                     >  1:500120/1‑140 (MQ=255)
                                                                                                                  tCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaa                             <  2:981838/132‑1 (MQ=255)
                                                                                                                  tCAATTAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACaaa                             >  1:981838/1‑132 (MQ=255)
                                                                                                                      ttAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATaaa                   <  1:27978/138‑1 (MQ=255)
                                                                                                                      ttAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATaaa                   >  2:27978/1‑138 (MQ=255)
                                                                                                                      ttAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGa                 >  2:351123/1‑140 (MQ=255)
                                                                                                                      ttAATGTTATTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGa                 <  2:1243378/140‑1 (MQ=255)
                                                                                                                              aTTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAg         >  1:402248/1‑140 (MQ=255)
                                                                                                                              aTTTTTAGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAg         <  2:1285608/140‑1 (MQ=255)
                                                                                                                                    aGTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCtt   <  1:888212/140‑1 (MQ=255)
                                                                                                                                     gTGGAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTc  <  1:386123/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
CTCAATGCTTCTAATATACTCGCAAAGATATCTAGAGGGCATCGTATGAGGTTTTGCTGCATTCTCTGAGCGAATATTAACTCCGTGCATATTATAGTTCAAAGCAAGTACCATTCAATTAATGTTATTTTTAGTGAAAAATTCTTTTATTTTATCGTCCGGGATTTTTCCGACTTGACTTTCAGTTCTGGGCTTAACTGTCAAAGCAAAGATAATACCACTCACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTC  >  NC_000913/1197419‑1197691

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