Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,985 C→A 17.2% intergenic (+87/‑383) 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,9850CA17.2% 131.5 / 3.8 61intergenic (+87/‑383)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (22/27);  new base A (7/5);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAG  >  NC_000913/2560842‑2561129
                                                                                                                                               |                                                                                                                                                
gCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGGTGTGGGTCaaaac                                                                                                                                             >  2:220301/1‑148 (MQ=255)
  aGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGtt                                                                                                                                           <  8:107710/149‑1 (MQ=255)
  aGAAAGGGAGGGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAACAGTTGCCGTTAGTGGTGGGTCAAAAGtt                                                                                                                                           <  6:308361/149‑1 (MQ=255)
   gAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATTTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCCTTTCTGGTGTGAAAAAAGTTg                                                                                                                                          >  5:198904/1‑149 (MQ=255)
    aaaGGGAGTGCAACGATTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGc                                                                                                                                         <  8:205033/149‑1 (MQ=255)
          aGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGTCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                   <  5:24396/149‑1 (MQ=255)
          aGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                   <  8:119991/149‑1 (MQ=255)
           gTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTGGCCGTTAAg                                                                                                                                  >  3:280187/1‑149 (MQ=255)
            tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAgg                                                                                                                                 >  1:75885/1‑149 (MQ=255)
              cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGTCAATCTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGACGTTATAGATGGGTAAAATTAAGACGTTAAGGt                                                                                                                                >  7:81723/1‑148 (MQ=255)
              cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                <  7:238346/148‑1 (MQ=255)
                tcGAATACCTGTTATTATTCAAGCCGTACAACATGAAATAAAGGAAGCAGGCGGGCATTGTTCCCGCCTTTTGTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCTTTAGTGGGGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                                             <  3:38730/148‑1 (MQ=255)
                aCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                                             <  1:238841/149‑1 (MQ=255)
                 cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCaaa                                                                                                                                               <  5:233283/130‑1 (MQ=255)
                 cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCaaa                                                                                                                                               >  6:233283/1‑130 (MQ=255)
                 cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                            <  2:598/149‑1 (MQ=255)
                 cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                            >  3:107898/1‑149 (MQ=255)
                       gcatttattaTTCAAGCCGTACAACATGAAATAACGGAATGAGGCGGGCAATGTTCCCGCCTTTTCTTTTCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                       <  4:347209/144‑1 (MQ=255)
                       ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTGAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                       <  4:107898/148‑1 (MQ=255)
                       ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                   <  6:3259/136‑1 (MQ=255)
                        cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTTAAGGTGGGGCAATAGt                                                                                                                     >  5:328700/1‑149 (MQ=255)
                        cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAATTTTGACGTTAAGGTGTGTGaaaatt                                                                                                                     >  1:59901/1‑147 (MQ=255)
                         aGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGACGTTAAGGTGGGTCAAATGtt                                                                                                                    >  7:85171/1‑149 (MQ=255)
                           ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAAGTTGCCTTTAAGGAGTGTCTAAAGTTGc                                                                                                                  >  1:323669/1‑149 (MQ=255)
                              ttattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAATGTTGCAGt                                                                                                               >  2:290874/1‑149 (MQ=255)
                                 ttCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGCATTGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGGTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTaa                                                                                                            <  1:220301/149‑1 (MQ=255)
                                 ttCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGACggt                                                                                                              >  2:76295/1‑145 (MQ=255)
                                    aaccagTACAACATTAAATAACGGAATAAGTCGTGCAATTTTAAAGCATTTTCTTTGCCTTTATTTCATGTATTCAGGTGGGTCAAAAGTTGGCGTTATTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAgg                                                                                                          <  2:75885/143‑1 (MQ=21)
                                    aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGGCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                         <  1:76295/149‑1 (MQ=255)
                                      gCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGAGGAGGGTAAAAAGTTGCAGTTAAGGAGGTTAAAAAGTTTAAGATAAGGTgg                                                                                                       >  1:178457/1‑149 (MQ=255)
                                             aacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTGGCGGTTAAGGTGGGTCAAATGGTGACGGTAAGGTGGGTCaa                                                                                                  >  3:268093/1‑147 (MQ=255)
                                             aacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTAAAAAGTTGCCTTGTAGGTGGGTAAAAATTTGAAGTTAACGGGTGTGaaaa                                                                                                >  2:271912/1‑149 (MQ=255)
                                               caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCaaaa                                                                                                                                              <  5:374927/101‑1 (MQ=255)
                                               caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCaaaa                                                                                                                                              >  6:374927/1‑101 (MQ=255)
                                                aTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGTGTCAAAATTTGCCTTTAAGGAGGGGCAAAATTTTCAGAGAAGGTGGCTGaaatttt                                                                                             >  8:308356/1‑145 (MQ=255)
                                                       aaCGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTTAAGGAGGGTCAAAAGTTGCAGTTTAGGGGGCGGAAGAGTTGCCGAg                                                                                       >  1:120744/1‑148 (MQ=38)
                                                          ggAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAAAGTTGCCGTTATGGGGGGTCAATAGtttactagt                                                                                      >  1:420320/1‑139 (MQ=39)
                                                            aaTCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAg                                                                                                                                             >  2:398366/1‑89 (MQ=255)
                                                            aaTCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAg                                                                                                                                             <  1:398366/89‑1 (MQ=255)
                                                                 ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCCGTTAAGGTGGGTCTAAATTTGACGTATACGAATTcg                                                                             >  3:290874/1‑147 (MQ=37)
                                                                                agCCGTTTCTATGCCTTAATTCCCTGAATTCAGGGGGGTCAAAGGGAGCCTTTAGTGGTGGGTCAATAGTTGCCGTTAAGGTGGGTCAAAAGTGGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa                                                             <  1:424716/148‑1 (MQ=25)
                                                                                                 atttacctGACTTCAGGTGCGTCAACAGCTGCCGTTAGCGGGGGGTCAAAAGTTGCCGTTAATGTGGGTCAAAAGTTTCATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGtttt                                              <  3:291690/140‑1 (MQ=25)
                                                                                                 aaTTCCATTCATTCAGGCGGGTCAAAAGTACGCGTTAGTGGTGGGTAAAAACCTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGtttttt                                            <  7:373045/149‑1 (MQ=37)
                                                                                                  aTTGCATGAATTCAGGGTCGCAAAAAGAAGCCGTTACTGGTGGCTCAAAACCTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                           <  1:271912/149‑1 (MQ=25)
                                                                                                  aTTCCATGAATTCAGGTGCGTAAAAATTTGCCGTTTGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGtttttt                                            <  8:418609/148‑1 (MQ=255)
                                                                                                                     ttCAAAAGTTGCCGTCAGTGCTGGGCCAAAAGTTGCCGTCACGGTGGGTCAAAAGTGGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGa                                                   <  1:152463/121‑1 (MQ=37)
                                                                                                                      taaaaaGATGGCCTTACTGGTGGGTAAATACTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATaga                       <  8:85171/147‑1 (MQ=37)
                                                                                                                      tCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCAGTTAAGGTGGGTCAAAATTTGCAGATAAGGTGGGGCAAAATTTGACTAGAAAGAATTGTAAAATTAAGAATTCTAGGAGACGTGAGTTTTTTAGACTAACGGTCGGTCATa                         >  6:30078/1‑147 (MQ=255)
                                                                                                                      tCAAAAGTTGCCGTTAGGGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg                                                           <  1:317036/113‑1 (MQ=255)
                                                                                                                        ataaGTTGCCTTTTTTGTTGGGTAAATAGTTGCCTTTTAGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagaga                     <  4:290874/147‑1 (MQ=38)
                                                                                                                         aaatgtGCCGTTAGTCGTTGTTCAAAAGTTGCGGTTAAGGTTTGTAAAAAGTGGCCGTTAAGGTGGTTAAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagaga                     <  8:411157/143‑1 (MQ=255)
                                                                                                                         aaaGTTGCCGTTAGTGGTGGGTAAAAATTTGACTTTAAGGAGTGTCAAAATTTTACGTTTAGGTGTGTCAAATTTTGACTAGTAAGAATTGTACTATTAAGAATTCTAGGCGAAGGGAGATTTTTAGAATCCCGTAAGCTCATaga                       >  7:411157/1‑146 (MQ=255)
                                                                                                                                ccGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGcc               >  8:174854/1‑147 (MQ=255)
                                                                                                                                  gTTACTGGCGGGTCAAAAGTGGCCTTTAAGGTGGGTCAAAAGTTCCCGTTCAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga           <  1:132761/149‑1 (MQ=255)
                                                                                                                                  cttACTGGTGGGTCAAAAGTAGCCGTTACGGTGGGTCAAAAGTTGCCGTTCAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga           <  1:114708/148‑1 (MQ=255)
                                                                                                                                          tggGTCAAAACTTGCCGTTAAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACa   <  7:219574/149‑1 (MQ=255)
                                                                                                                                          tcggTAAAAAGTTGCCGTTAATGTGGGTCAATCGGTGCCGGTTCGGTGGGTCAAAAGTTGCCGAGTAACAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCa                           <  1:66130/123‑1 (MQ=255)
                                                                                                                                            ggTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGCCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag  >  8:125579/1‑148 (MQ=255)
                                                                                                                                             gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag  >  3:199939/1‑147 (MQ=255)
                                                                                                                                              tCAATAGGTGGCGTTAAGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAATGTTGCCGAGTAAGATTTCTACTTTGAAGAATTCTAGACGACGGGAGTGTTTTAGCCTCACGGTCGGTCTGAGAGAGTCTCAGAGATTTAAagag  >  6:369023/1‑146 (MQ=255)
                                                                                                                                              tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag  >  8:307055/1‑146 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAG  >  NC_000913/2560842‑2561129

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: