Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,996 G→T 14.9% intergenic (+98/‑372) yffL → / → yffM CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,9960GT14.9% 124.7 / 2.5 59intergenic (+98/‑372)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (31/18);  new base T (5/5);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCC  >  NC_000913/2560851‑2561144
                                                                                                                                                 |                                                                                                                                                    
gAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGtt                                                                                                                                                    >  3:137229/1‑148 (MQ=255)
   tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTTGGTAAAATTTTGCCGTTAAgg                                                                                                                                                >  3:35029/1‑149 (MQ=255)
     cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                               <  8:277616/148‑1 (MQ=255)
     cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGTGTGGTTGGTCATATTTTGCGGtttttgtg                                                                                                                                              >  5:306026/1‑143 (MQ=255)
     cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATAAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTg                                                                                                                                              >  7:119474/1‑149 (MQ=255)
      aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACTGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                                             <  3:314896/149‑1 (MQ=255)
         gAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGt                                                                                                                                           >  7:46488/1‑148 (MQ=255)
            tACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGTGGGGTAAAAAGTTGCCGTTAAGGTGTGTaaa                                                                                                                                        >  5:246290/1‑148 (MQ=255)
              ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTTGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                                      <  5:159866/148‑1 (MQ=255)
              ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                                      <  7:70480/148‑1 (MQ=255)
               cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGt                                                                                                                                    >  4:263266/1‑149 (MQ=255)
                  ttattaTTTAAGCCGTACCACATGAATTAACTGAATCCGGCGGGCAATGTTCCGGCCTTTTATTTGCATTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTTGTGGTTAAATGTTGCATTTAAGGTGGGTCAAAAGTTGcc                                                                                                                                <  4:256519/149‑1 (MQ=25)
                    attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTCAAAATTTGCAGTTAAGTAGGGTCAAAAGtttacg                                                                                                                               >  1:412538/1‑145 (MQ=255)
                    attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                                                                                                                <  3:263266/148‑1 (MQ=255)
                    attaatCAAGCAGTAGAACATGAAATAACGGAATCAGTCTGGCAATGTTCCAGCATTTTCTTTTCCTTAATTCCATGAATTCAGGTGGGTAAAAAGTTGCCCTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                                                                                                                <  2:380332/148‑1 (MQ=37)
                     ttattCAACCCGTACAACATGAAATTAAGGAATCAGTCTTGCAATTTTCCCGCATTTTGTTTTCCTTAATTCCATGAATTCAGGTGGGTAAAAAGTTGCCGTTAGTGTTGGGTCAAAAGTTTCAGTTAAGGTGGGTCAAAAGTTGCCg                                                                                                                               <  4:137229/148‑1 (MQ=25)
                          cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGAAGTTAAgg                                                                                                                         >  7:315037/1‑149 (MQ=255)
                           atgcggTACTACATGAAATTACGGAATCAGGGGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAg                                                                                                                          <  3:99467/142‑1 (MQ=255)
                           aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCAGTTAGTGGTCTGTCAAAAGTTGCCGTTAAGGACGGGCAAAAGATGCAGTTAAgg                                                                                                                         >  4:99467/1‑148 (MQ=255)
                                 tACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGAGTGGTGGGGCAAAAGTTGCGGGTTAGGAGGGGCAAAAGTTGCAGTTAAGGGGGGTc                                                                                                                   >  1:293396/1‑148 (MQ=255)
                                      caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGTGTCAAAAGTTGCCTTTAAGGAGGGTCAAAATTTTCCGTGTAGGTGGGTCAAAAGt                                                                                                             >  3:221509/1‑149 (MQ=255)
                                      caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                                             <  8:280395/117‑1 (MQ=255)
                                      caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                                             >  7:280395/1‑117 (MQ=255)
                                           aaTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                                      <  1:385054/119‑1 (MQ=255)
                                           aaTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                                      >  2:385054/1‑119 (MQ=255)
                                               aCGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGTGTGTTGGGTCAAATGTTGACTTGAAGGAGCGGCTATTGTTGTCGTGAAGGTGGGGGTAATGTTGCGGAgg                                                                                                     >  3:159853/1‑147 (MQ=255)
                                                             gCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTGAAGGAGGGTCAAAAGTTGACGTGTAGGTGGCTGTAAAGTTGAGGAGAAAGCATTctacaag                                                                                       >  8:239699/1‑145 (MQ=37)
                                                                       cGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAATTAGACTTTAATGTGTCTCAAAATTTTAAGTTAAGGTGTGTGAAAATTAGCATAGTAAAAATTATAAAATTATGAATACTa                                                                            >  1:240852/1‑149 (MQ=255)
                                                                                       tAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCAGTTTACGGGGGTCAAAAGTTTCCGAGTAAGAATTCGACTATTAAGAATTCTATGCGTCGGGAGtgttg                                                            >  5:294915/1‑148 (MQ=25)
                                                                                             cATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGACAAAAGTTGCCTTTAAGGAGGGTCAAAATTTGCAGTATAGGTGGCTCAAAAGTTGACGAGGAAGAATTCTAATATTATCAATGCTAGGCGACGGGAGTTTTTTATCCt                                                      >  8:68227/1‑149 (MQ=25)
                                                                                                   ttCAGCTGGGCAAAAAGCAGCCCTTACGGGTGGCTAAAAAGAGTCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggt                                                <  5:250320/149‑1 (MQ=25)
                                                                                                        ggtggTCAAAAGTTGCCGTTAGTGGTGGGTCAAAACTTGCCGTTCAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCggtcggtc                                           <  8:157552/146‑1 (MQ=255)
                                                                                                                  agctggcGTTATTGGGGGGTCATTAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTc                                 <  8:198529/143‑1 (MQ=38)
                                                                                                                    ttGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGAGGGTCAAAAGTTGCCGTTAAGGTCGGTCAAAAGTTTCCGAGTAAGAATACTACTATTAATAATTCAAGGCTACTGGAGTTTTTTAGCATACAGGACGGAAATAGAGAGtag                                >  8:201499/1‑146 (MQ=21)
                                                                                                                         gTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa                            >  8:136440/1‑147 (MQ=255)
                                                                                                                             tttgtggGTAATAAGTTGCCGTTAAGGTGGGTAAAATGTTTCCTTTATGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGtt                       <  8:315037/145‑1 (MQ=255)
                                                                                                                             gtggtgGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGtt                       >  6:325550/1‑148 (MQ=255)
                                                                                                                             gtggtgGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAAGTCGGTCAAAAGTTGCCTAGGAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTAAGCCTCCCGCTCGGTCATAGAGAGTCGCCTAGAGtg                       >  4:340717/1‑147 (MQ=37)
                                                                                                                                gtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCTCCTAGAGTTaaa                    >  6:33690/1‑148 (MQ=255)
                                                                                                                                gtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGACTTTTTTAGCCTCCAGGTCGGTCATAGAGAGACGCCTAGAGTTaaa                    >  6:322907/1‑148 (MQ=255)
                                                                                                                                 tggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaaa                    >  1:218640/1‑147 (MQ=255)
                                                                                                                                   ggTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag                 <  5:322907/148‑1 (MQ=255)
                                                                                                                                    gtaaaaaGTTGCCTTTAAGGTGGGTAAAAAGTAGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaa               <  5:315159/146‑1 (MQ=255)
                                                                                                                                    gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaa               <  7:68227/149‑1 (MQ=255)
                                                                                                                                     tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAATAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag                 >  7:100053/1‑146 (MQ=255)
                                                                                                                                     tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag                 >  4:334397/1‑146 (MQ=255)
                                                                                                                                     tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag                 >  6:765/1‑146 (MQ=255)
                                                                                                                                      cAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaag              >  2:138416/1‑148 (MQ=255)
                                                                                                                                       aaaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGtttttt                                                           <  1:16350/102‑1 (MQ=255)
                                                                                                                                       aaaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGtttttt                                                           >  2:16350/1‑102 (MQ=255)
                                                                                                                                        aaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                                          <  1:411474/102‑1 (MQ=255)
                                                                                                                                        aaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                                          >  2:411474/1‑102 (MQ=255)
                                                                                                                                        aaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAat           >  6:298651/1‑149 (MQ=255)
                                                                                                                                         aaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatt          >  6:224508/1‑149 (MQ=255)
                                                                                                                                          aggtGCCTTTTAGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatta         <  8:119474/146‑1 (MQ=255)
                                                                                                                                           caggCCGTTAAGGTGTGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattat        <  5:362890/146‑1 (MQ=255)
                                                                                                                                              gCCGTTAATGAGCTTCAAAAGTTGAAGTAAAGGTGGGACAAAAGTTTCCGAGTAAGAATACTACTATTAACAATTCTAGGCGACGCGAGTTTTTTAGCCTCCCGGTCGGTCAAAGAGATTAGCCTATAGTTAAACATAAGAaatattac     >  6:91297/1‑147 (MQ=255)
                                                                                                                                               ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                                          <  3:137877/95‑1 (MQ=255)
                                                                                                                                               ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                                          >  4:137877/1‑95 (MQ=255)
                                                                                                                                                 gTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgcc  >  1:373483/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCC  >  NC_000913/2560851‑2561144

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: