Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,771,183 A→C 6.2% A12A (GCA→GCC)  yfjT → CP4‑57 prophage; putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,771,1830AC 6.2% 327.6 / 11.9 97A12A (GCA→GCC) yfjTCP4‑57 prophage; putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (41/50);  new base C (2/4);  total (43/54)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

TACGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTA  >  NC_000913/2771037‑2771331
                                                                                                                                                  |                                                                                                                                                    
taCGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCGGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGCTTGTACTGGCaac                                                                                                                                                    <  4:827195/149‑2 (MQ=255)
   gaaagaaaAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTa                                                                                                                                                  <  2:581948/148‑1 (MQ=255)
       gaaaAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaa                                                                                                                                              <  2:479770/148‑2 (MQ=255)
       gaaaAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAAt                                                                                                                                              >  2:306869/1‑148 (MQ=255)
          aaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCTTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGtt                                                                                                                                           <  2:758847/148‑1 (MQ=255)
          aaaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGtt                                                                                                                                           >  1:739767/1‑148 (MQ=255)
           aaTACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaaccggc                                                                                                                                         <  3:150148/149‑7 (MQ=255)
               cTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCa                                                                                                                                      >  1:106682/1‑148 (MQ=255)
                ttCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGc                                                                                                                                    >  4:562691/1‑149 (MQ=255)
                 tCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCATCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCt                                                                                                                                   <  3:188026/149‑1 (MQ=255)
                   gAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaac                                                                                                                                             <  4:863382/137‑3 (MQ=255)
                   gAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTTGCCAGTACaaac                                                                                                                                             >  3:863382/1‑135 (MQ=255)
                            cATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCctc                                                                                                                         >  2:508338/1‑148 (MQ=255)
                               gACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCg                                                                                                                     <  4:574416/149‑1 (MQ=255)
                                   gCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTcc                                                                                                                 >  3:620262/1‑149 (MQ=255)
                                     aaaCTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTccc                                                                                                                <  1:20419/148‑1 (MQ=255)
                                         tGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGa                                                                                                            <  2:106682/148‑1 (MQ=255)
                                          gATGGGGCGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGAc                                                                                                           >  2:517903/1‑148 (MQ=255)
                                            tGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGActc                                                                                                         <  1:296322/148‑1 (MQ=255)
                                             gggACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGActct                                                                                                        <  2:278685/148‑1 (MQ=255)
                                              ggACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTg                                                                                                       <  1:714477/148‑1 (MQ=255)
                                               gACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTgg                                                                                                      <  1:755239/148‑1 (MQ=255)
                                                    ctgctgATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAg                                                                                                <  3:322372/149‑1 (MQ=255)
                                                    ctgctgATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAg                                                                                                <  3:848000/149‑1 (MQ=255)
                                                      gctgATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATTTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAggg                                                                                              <  4:465076/149‑1 (MQ=255)
                                                      gctgATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAggg                                                                                              <  4:872889/149‑1 (MQ=255)
                                                        ttaTACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGt                                                                                            >  4:459592/3‑149 (MQ=255)
                                                               gTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCa                                                                                     >  4:797317/1‑149 (MQ=255)
                                                                tCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAg                                                                                    <  4:751061/149‑1 (MQ=255)
                                                                tCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAg                                                                                    >  3:948297/1‑149 (MQ=255)
                                                                 cAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAgg                                                                                   <  3:375298/149‑1 (MQ=255)
                                                                     cccTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCaa                                                                               >  4:580054/1‑149 (MQ=255)
                                                                     cccTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCa                                                                                <  2:341209/148‑1 (MQ=255)
                                                                      ccTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCaa                                                                               >  2:875367/1‑148 (MQ=255)
                                                                        tGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGc                                                                            >  3:811088/1‑149 (MQ=255)
                                                                         gATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCg                                                                           <  3:413091/149‑1 (MQ=255)
                                                                          aTAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCg                                                                           >  2:262207/1‑148 (MQ=255)
                                                                           tAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGACAATATGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGa                                                                         >  3:862620/1‑149 (MQ=255)
                                                                              ccGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAAt                                                                       <  1:875367/148‑1 (MQ=255)
                                                                              ccGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAAt                                                                       <  2:831950/148‑1 (MQ=255)
                                                                                gggAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGActct                                                                                                        <  3:934193/113‑1 (MQ=255)
                                                                                gggAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGActct                                                                                                        >  4:934193/1‑113 (MQ=255)
                                                                                 ggAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACa                                                                    <  1:84294/148‑1 (MQ=255)
                                                                                 ggAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAg                                                                   <  4:626672/149‑1 (MQ=255)
                                                                                  gAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAgt                                                                  <  4:811088/149‑1 (MQ=255)
                                                                                      aaaaaCAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGATTACAGtgtg                                                               >  2:964690/1‑148 (MQ=255)
                                                                                      aaaaaCAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGa                                                              <  3:747900/149‑1 (MQ=255)
                                                                                         aacaacGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACgg                                                           >  3:816208/1‑149 (MQ=255)
                                                                                         aacaacGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACg                                                            >  2:114796/1‑148 (MQ=255)
                                                                                          acaacGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGa                                                          <  3:718905/149‑1 (MQ=255)
                                                                                            aacGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATg                                                        >  4:690747/1‑149 (MQ=255)
                                                                                                 gaAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGc                                                   >  3:799376/1‑149 (MQ=255)
                                                                                                 gaAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGc                                                   >  4:548824/1‑149 (MQ=255)
                                                                                                   aaGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCg                                                  >  1:274773/1‑148 (MQ=255)
                                                                                                      gagaCATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGaa                                               <  1:381670/148‑1 (MQ=255)
                                                                                                         aCATAAATTAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGCCTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGAt                                           <  4:112457/149‑1 (MQ=255)
                                                                                                         aCATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGAt                                           <  3:757944/149‑1 (MQ=255)
                                                                                                         aCATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGCGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGAt                                           <  3:719199/149‑1 (MQ=255)
                                                                                                           aTAAATGAAAATACTTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATg                                          <  1:94567/148‑1 (MQ=255)
                                                                                                           aTAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGt                                         >  3:626847/1‑149 (MQ=255)
                                                                                                            tAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGTGTACGGTCAGGGCAATGCGGAATACATTGTGACGGATGCCAGCGGAAAgttgtt                                        >  4:737672/1‑149 (MQ=255)
                                                                                                             aaaTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGttt                                       <  4:734593/149‑1 (MQ=255)
                                                                                                               aTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTAc                                     >  3:556388/1‑149 (MQ=255)
                                                                                                               aTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTAc                                     <  4:799376/149‑1 (MQ=255)
                                                                                                                 gAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCa                                   <  3:548824/149‑1 (MQ=255)
                                                                                                                    aaTACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGaaa                                              >  3:164918/1‑135 (MQ=255)
                                                                                                                    aaTACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGaaa                                              <  4:164918/135‑1 (MQ=255)
                                                                                                                    aaTACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATg                                          >  2:232586/1‑139 (MQ=255)
                                                                                                                    aaTACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATg                                          <  1:232586/139‑1 (MQ=255)
                                                                                                                       gccatCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGTACTCTGGCAATAGGGGTACGGTCAGGGCGATGCGGAATACAGTGTGACGGATGCGAGCGGAAAGATGTTTACCGTCAAc                              <  2:841108/144‑1 (MQ=255)
                                                                                                                       aCGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGAAAAGATGTTTACCATCAAc                              <  1:747394/148‑1 (MQ=255)
                                                                                                                         gTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAgt                                                                  <  3:34026/110‑1 (MQ=255)
                                                                                                                         gTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAgt                                                                  >  4:34026/1‑110 (MQ=255)
                                                                                                                         gTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGc                           >  3:178344/1‑149 (MQ=255)
                                                                                                                           tCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGa                                                          <  1:148207/116‑1 (MQ=255)
                                                                                                                           tCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGa                                                          >  2:148207/1‑116 (MQ=255)
                                                                                                                             aCTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACg                        <  1:705074/148‑1 (MQ=255)
                                                                                                                              cTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACgg                       >  1:645708/1‑148 (MQ=255)
                                                                                                                                aaGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTa                                      <  1:696517/131‑1 (MQ=255)
                                                                                                                                aaGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTa                                      >  2:696517/1‑131 (MQ=255)
                                                                                                                                 aGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGaa                    <  1:262207/148‑1 (MQ=255)
                                                                                                                                  gCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAAt                  <  3:690747/149‑1 (MQ=255)
                                                                                                                                   ccGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATc                 <  3:49021/149‑1 (MQ=255)
                                                                                                                                     ggTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATcc                <  2:934838/148‑1 (MQ=255)
                                                                                                                                      gTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTTACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCgg              >  3:739505/1‑149 (MQ=255)
                                                                                                                                      gTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCgg              >  4:942468/1‑149 (MQ=255)
                                                                                                                                       tttGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGa             <  3:942468/149‑1 (MQ=255)
                                                                                                                                       tttGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGa             <  4:178344/149‑1 (MQ=255)
                                                                                                                                       tttGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGGCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCgg              <  2:184044/148‑1 (MQ=255)
                                                                                                                                       attGTACTTGCCAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGa                                                              >  4:605894/2‑100 (MQ=255)
                                                                                                                                       attGTACTTGCCAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGCGGTACGGTCAGGGCAATGCGGAATACAGTGTGa                                                              <  3:605894/99‑1 (MQ=255)
                                                                                                                                            aCTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGa        >  4:799703/1‑149 (MQ=255)
                                                                                                                                              tGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCgg                                                 >  2:925669/1‑106 (MQ=255)
                                                                                                                                              tGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCgg                                                 <  1:925669/106‑1 (MQ=255)
                                                                                                                                              tGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAAc      >  4:660586/1‑149 (MQ=255)
                                                                                                                                              tGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAAc      >  4:5612/1‑149 (MQ=255)
                                                                                                                                               ggCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACg     >  3:12188/1‑149 (MQ=255)
                                                                                                                                                 cAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACgg    >  1:595716/1‑148 (MQ=255)
                                                                                                                                                  aaGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGt   <  2:295705/148‑1 (MQ=255)
                                                                                                                                                  aaGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTa  <  4:556388/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TACGAAAGAAAAATACTTCGAATATTCCCATGACGGCAAACTGATGGGACGGCTGCTGATACCGTCAAACCCTGATAACCGGGAATAAAAACAACGAGAAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTA  >  NC_000913/2771037‑2771331

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