Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,415,148 G→T 20.5% P747T (CCG→ACG)  recE ← Rac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,415,1480GT20.5% 230.4 / 53.0 93P747T (CCG→ACG) recERac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base T (11/8);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-02

TGCACTCCAAACTGTGCTTCATAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCA  >  NC_000913/1415000‑1415293
                                                                                                                                                    |                                                                                                                                                 
tGCACTCCAAACTGTGCTTCATAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCgg                                                                                                                                                   >  1:280869/1‑149 (MQ=255)
     tCCAAACTGTGCTTCATAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGc                                                                                                                                              >  3:681514/1‑149 (MQ=255)
     tCCAAACTGTGCTTCATAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGc                                                                                                                                              >  3:311236/1‑149 (MQ=255)
     tCCAAACTGTGCTTCATAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGc                                                                                                                                              >  3:46971/1‑149 (MQ=255)
            tgtgCTTCATAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGaca                                                                                                                                       <  2:281661/149‑1 (MQ=255)
               gCTTCATAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACaa                                                                                                                                    <  3:823621/149‑1 (MQ=255)
                 ttCATAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACaaaa                                                                                                                                  >  3:280956/1‑149 (MQ=255)
                    aTAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTc                                                                                                                               <  1:881971/149‑1 (MQ=255)
                    aTAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTc                                                                                                                               >  2:628517/1‑149 (MQ=255)
                    aTAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTc                                                                                                                               <  2:280869/149‑1 (MQ=255)
                    aTAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTc                                                                                                                               >  1:422299/1‑149 (MQ=255)
                     tAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTcc                                                                                                                              >  2:80919/1‑149 (MQ=255)
                     tAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTcc                                                                                                                              >  4:599867/1‑149 (MQ=255)
                     tAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTcc                                                                                                                              <  1:638004/149‑1 (MQ=255)
                      aaCCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCt                                                                                                                             <  2:214932/149‑1 (MQ=255)
                               tGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGa                                                                                                                    <  3:599867/149‑1 (MQ=255)
                               tGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGa                                                                                                                    >  2:749826/1‑149 (MQ=255)
                               tGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGa                                                                                                                    >  1:6214/1‑149 (MQ=255)
                                gTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCTTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGAt                                                                                                                   <  1:661735/149‑1 (MQ=255)
                                 tAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATc                                                                                                                  >  3:283260/1‑149 (MQ=255)
                                 tAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATc                                                                                                                  <  4:505992/149‑1 (MQ=255)
                                 tAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATc                                                                                                                  >  3:179439/1‑149 (MQ=255)
                                    aaTGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTc                                                                                                               <  4:311236/149‑1 (MQ=255)
                                     aTGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTcc                                                                                                              <  2:6214/149‑1 (MQ=255)
                                     aTGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTcc                                                                                                              <  1:303207/149‑1 (MQ=255)
                                     aTGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTcc                                                                                                              <  3:317992/149‑1 (MQ=255)
                                      tGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTccc                                                                                                             <  2:337456/149‑1 (MQ=255)
                                       gCATCCTGAACGTGGTAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATcatggac                                                                                                            <  2:178179/149‑7 (MQ=255)
                                         aTCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAgg                                                                                                                     >  2:851768/1‑138 (MQ=255)
                                         aTCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAgg                                                                                                                     <  1:851768/138‑1 (MQ=255)
                                         aTCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGt                                                                                                          >  4:777250/1‑149 (MQ=255)
                                          tCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTa                                                                                                         >  4:521805/1‑149 (MQ=255)
                                             tGAACGGGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAAt                                                                                                      >  1:810882/1‑149 (MQ=255)
                                              gAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAAtt                                                                                                     <  4:46971/149‑1 (MQ=255)
                                               aaCGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATtg                                                                                                    >  4:396001/1‑149 (MQ=255)
                                                 cGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATtgat                                                                                                  <  4:845718/149‑1 (MQ=255)
                                                 cGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATtgat                                                                                                  <  3:396001/149‑1 (MQ=255)
                                                 cGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATtgat                                                                                                  <  4:266392/149‑1 (MQ=255)
                                                   tGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATtgatga                                                                                                <  3:518487/149‑1 (MQ=255)
                                                      tAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTc                                                                                             <  2:267516/149‑1 (MQ=255)
                                                       aGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCa                                                                                            >  4:150772/1‑149 (MQ=255)
                                                         cGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGc                                                                                          <  4:323450/149‑1 (MQ=255)
                                                               tCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTccc                                                                                                             >  2:446437/1‑124 (MQ=255)
                                                               tCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTccc                                                                                                             <  1:446437/124‑1 (MQ=255)
                                                                 gTAATAAGCGGTTTTGAAGCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCATCCCAGTAAATTGATGATTCAGCGTGTCCgg                                                                                  >  1:273054/1‑149 (MQ=255)
                                                                   aataaGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGgcg                                                                                >  4:290437/1‑149 (MQ=255)
                                                                   aataaGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACAGAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGgcg                                                                                >  3:798019/1‑149 (MQ=255)
                                                                      aaGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAAt                                                                                                      <  1:145917/124‑1 (MQ=255)
                                                                      aaGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAAt                                                                                                      >  2:145917/1‑124 (MQ=255)
                                                                      aaGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCtt                                                                             <  3:862200/149‑1 (MQ=255)
                                                                      aaGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCtt                                                                             >  1:365554/1‑149 (MQ=255)
                                                                       aGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCGGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCttt                                                                            >  4:586429/1‑149 (MQ=255)
                                                                       aGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCttt                                                                            <  4:657061/149‑1 (MQ=255)
                                                                           gTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGTTCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTcaac                                                                        >  4:699547/1‑149 (MQ=255)
                                                                             tttGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTcaacaa                                                                      <  4:321981/149‑1 (MQ=255)
                                                                                 aaTCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCa                                                                  >  2:435842/1‑149 (MQ=255)
                                                                                     gTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATtgat                                                                                                  <  3:732287/113‑1 (MQ=255)
                                                                                     gTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATtgat                                                                                                  >  4:732287/1‑113 (MQ=255)
                                                                                     gTTCAATACCCGCCGTAGTTTTCACGTCCATGATCCAGTGAACTTCAGGGATAATTAGGTTCGGACGGCTCCGACTCAAGATTCCTGTTTCAGGATCTTCCCAGTAAATTGACGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGc                                                              <  3:17147/149‑1 (MQ=255)
                                                                                      ttGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATtga                                                                                                   <  4:730821/111‑1 (MQ=255)
                                                                                      ttGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGcc                                                             >  3:654828/1‑149 (MQ=255)
                                                                                      ttGAATATCCGCCGTAGTTCTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATtga                                                                                                   >  3:730821/1‑111 (MQ=255)
                                                                                            atCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCg                                                       >  4:89178/1‑149 (MQ=255)
                                                                                                gccgTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCaa                                                   >  1:147126/1‑149 (MQ=255)
                                                                                                 ccgTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCaaa                                                  >  4:204697/1‑149 (MQ=255)
                                                                                                 ccgTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCaaa                                                  <  4:843943/149‑1 (MQ=255)
                                                                                                 ccgTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCaaa                                                  <  2:545070/149‑1 (MQ=255)
                                                                                                    tAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGcc                                               >  1:644505/1‑149 (MQ=255)
                                                                                                       ttttCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATa                                            >  1:157235/1‑149 (MQ=255)
                                                                                                       ttttCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATa                                            >  2:916203/1‑149 (MQ=255)
                                                                                                       ttttCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATa                                            >  2:751922/1‑149 (MQ=255)
                                                                                                          tCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACg                                         >  1:821773/1‑149 (MQ=255)
                                                                                                            aCGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCt                                       >  3:93018/1‑149 (MQ=255)
                                                                                                            aCGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCt                                       >  1:26784/1‑149 (MQ=255)
                                                                                                             cGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCtt                                      <  1:749826/149‑1 (MQ=255)
                                                                                                                ccATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGa                                   <  2:753534/149‑1 (MQ=255)
                                                                                                                 cATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTcc                                                                                                                              <  2:305314/57‑1 (MQ=255)
                                                                                                                 cATGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTcc                                                                                                                              >  1:305314/1‑57 (MQ=255)
                                                                                                                  aTGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATa                                 <  1:435842/149‑1 (MQ=255)
                                                                                                                   tGATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGacac                                >  3:613709/1‑149 (MQ=255)
                                                                                                                    gATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATa                                 >  2:607939/1‑147 (MQ=255)
                                                                                                                    gATCCAGTGAAATTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATa                                 <  1:607939/147‑1 (MQ=255)
                                                                                                                    gATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACa                               <  2:147126/149‑1 (MQ=255)
                                                                                                                     aTCCAGTGAAATTCACGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGTCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACAt                              <  2:273054/149‑1 (MQ=255)
                                                                                                                      tCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATg                             >  3:147174/1‑149 (MQ=255)
                                                                                                                               aTTCAGGGATAATTTTGTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAAtt                    >  1:469785/1‑149 (MQ=255)
                                                                                                                                  cAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTc                 >  1:504455/1‑149 (MQ=255)
                                                                                                                                      gATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGc                                                     >  3:555479/1‑109 (MQ=255)
                                                                                                                                      gATAATTTTGTCCGGACGGCACCGACACAAAATTCCCGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGc                                                     <  4:555479/109‑1 (MQ=255)
                                                                                                                                      gATAATTTTGTCCGGACGGCACCGACACAAAATTCCAGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGc             >  4:875160/1‑149 (MQ=255)
                                                                                                                                            tttGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCcttctt       <  3:777250/149‑1 (MQ=255)
                                                                                                                                               gTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCaaa                                                  >  2:285861/1‑103 (MQ=255)
                                                                                                                                               gTCCGTACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCaaa                                                  <  1:285861/103‑1 (MQ=255)
                                                                                                                                                 ccGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCa  <  2:26784/149‑1 (MQ=255)
                                                                                                                                                    gACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACa                               <  4:226504/117‑1 (MQ=255)
                                                                                                                                                    gACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACa                               >  3:226504/1‑117 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGCACTCCAAACTGTGCTTCATAACCGTCACTGTAGAATGCATCCTGAACGTGATAGCGGTAGTCGTAATAAGCGGTTTTGAATCGTTGAATATCCGCCGTAGTTTTCACGTCCATGATCCAGTGAAATTCAGGGATAATTTTGTCCGGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCA  >  NC_000913/1415000‑1415293

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: