Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,561,012 A→T 14.3% intergenic (+114/‑356) yffL → / → yffM CPZ‑55 prophage; uncharacterized protein YffL/CPZ‑55 prophage; uncharacterized protein YffM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,561,0120AT14.3% 119.5 / 2.4 56intergenic (+114/‑356)yffL/yffMCPZ‑55 prophage; uncharacterized protein YffL/CPZ‑55 prophage; uncharacterized protein YffM
Reads supporting (aligned to +/- strand):  ref base A (28/19);  new base T (5/3);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTG  >  NC_000913/2560864‑2561161
                                                                                                                                                    |                                                                                                                                                     
aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAg                                                                                                                                                      >  1:478235/1‑150 (MQ=255)
 ccAGTTATTATTCAAGCCGTACAACATGAAATAAAGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGACGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTAAAAAGt                                                                                                                                                     >  1:710177/1‑150 (MQ=255)
   aGTTATTATTCAAGCCGTACAACATGAAATAACGTAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                                                                                                                   <  1:351406/150‑1 (MQ=255)
   aGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                                                                                                                   <  1:346053/150‑1 (MQ=255)
     ttattaTTGAAGCCGTATAACATGAAATAACGGAATCAGGCTTGCAATTTTTCCGCCTTTTGTTTTCCTTTATTCCATGAATTCAGGTGGGTAAAATTTTGCCGTTTGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGcc                                                                                                                                                 <  2:478235/150‑1 (MQ=25)
       attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTTGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGt                                                                                                                                               >  1:693410/1‑150 (MQ=255)
         tattCAAGCCGTACAACATGAAATAACGGAATCAGGCTGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTa                                                                                                                                             <  1:1119191/150‑1 (MQ=255)
          attCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAATTTTGCCGTGTAGGAGGGGCAATTGTTGGCgggaa                                                                                                                                            >  1:1313835/1‑146 (MQ=255)
             cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                         >  1:1327556/1‑150 (MQ=255)
             cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCAATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTAGCCGTTAAGGTGGGTCAAATGTTGACGTTAAGGt                                                                                                                                         >  1:470236/1‑150 (MQ=255)
              aaGCCGTACATCATGAAATTACGGAATAAGGCGGGCAATGTTCCCTCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAATTTGCCGTTAGTGGTGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTg                                                                                                                                        <  2:1143931/150‑1 (MQ=255)
                tccttACAAGATTAAATAACGGAATCAGGCGGGCAATGTTCCCGCATTTTCTTTGCCTTATTTCCATGAATTCAGGTGGGTTAAAATTTGGCGTTAGTGTTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                                                      <  2:1327556/146‑1 (MQ=38)
                gCCGTACAACATGAATTATCGGATTCAGGCGTGCATTGTTCCCTAATTTTCTTTGACTTTATTCCATGAATTCAGTTGGGGAAAAAGTTGCCGTTAGTGGTTGGTTAATAGTTTCCGTTTAGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                                                      <  2:662837/150‑1 (MQ=255)
                  cGTACAACATCAAATAACGGAATCCGTCGGGCACTTTTCCCGCCTTTTGTTTGCCATAAATCCATGAATTCTGGTGGGTAAAAATTTGCCTTTATTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                                                                                                    <  2:1172300/150‑1 (MQ=25)
                         caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTGAAGGTGGGTCAAAAGTTGCCGTGAAGGTGGGGGAAATGtt                                                                                                                             >  1:915710/1‑150 (MQ=255)
                         caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGACGTTAAGGTGGGTCAAATGTTGGAGTGTAGGGGGGTCGAAAGtt                                                                                                                             >  1:421856/1‑150 (MQ=255)
                         caTGAAATAACGGAATCAGGCCGGCAATGTTCCCGCCTTTTCTTTTCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                                             <  1:277123/150‑1 (MQ=255)
                           tGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTACTGGTGGGTCAAAAGTAGCCGTTAAGGAGGGTCAAATGTTGAAGTTAAGGTGGGGCAAAAGTAGc                                                                                                                           >  1:582687/1‑150 (MQ=255)
                             aaaTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGGGGGTCAAAATTTGCAGTTTAGGTGGGTGAAATGTTgagg                                                                                                                         >  1:580197/1‑147 (MQ=255)
                                     gAATCCGGCGGGCAATCTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGGGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGaa                                                                                                                 <  1:974407/150‑1 (MQ=255)
                                           ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTGAAAAGTTGAATAGTAAGAATTctaa                                                                                                           >  1:107099/1‑149 (MQ=255)
                                                 cAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCATTTAAGGTGGGTCAAAAGTTGCAGTAAACGTGGGTCAAAATTTGAAGAGTAAGAATTGTAAAATTaa                                                                                                     >  2:22656/1‑150 (MQ=25)
                                                       tCCAGCCTTTTGTTTTCCTTTAATTAATGAATTCAGGTTGGTCAAAAGCTGGCTTTTTTTTTGGGTAATAACTTGCCTTTAAGGTGGGTAAAATGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTc                                                                                               <  2:1185419/150‑1 (MQ=255)
                                                        cccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGTTGGGTCAAAAGTTGCATTTAAGGAGGCTCAAATTTTGGCGTATAGGTGGGGGAAATGTTGACTAGTAAGAATTGTACTATTATTAAttgt                                                                                              >  1:8906/1‑148 (MQ=255)
                                                         ccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACGAGTAAGAATTCTACTATTAAGAATTCTa                                                                                             >  1:547398/1‑150 (MQ=255)
                                                           gCCTTTTGTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGCTGCCGTTAGGGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg                                                                                           <  2:951515/150‑1 (MQ=255)
                                                            ccTTTTCTTTTGCTTAATTCCATGAATTCAGGTTGGTAAAAAGTAGCCTTTACTGGTGGGTCAAAAGTTGCCTTTAAGGTGTGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGc                                                                                          <  1:917971/150‑1 (MQ=37)
                                                               tttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTGAAAATTTGCAGAGTAAGAATTGTTCTTTTAAGATTTCTAGGCGag                                                                                       >  1:846745/1‑149 (MQ=37)
                                                                 tCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACAATTAAGAATTCTAGGCGACgg                                                                                     >  2:775929/1‑150 (MQ=255)
                                                                  cTTTCCCTTAATTCCATGAATTCAGTGGGGTCAAAAGAACCCCTTATGGGTGGGTCAAATCCTGCCTTTACGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACggg                                                                                    <  1:1248752/150‑1 (MQ=25)
                                                                         ttAATTACATGAATTCTGTTGGGTTATTTTTTGGCCTTTTTGTTGGCTAATATGTTGCCTTTTAGGTGGGTTAAAAGTTGCCTTTAAGGTGTGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGttttt                                                                             <  2:337768/150‑1 (MQ=255)
                                                                             ttCCATGAATTCAGGTGGGTCAACAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAGGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGc                                                                         <  1:1150410/150‑1 (MQ=255)
                                                                             ttCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCAGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGCATTCTACGATTATGAATTCTAGGCGACGGGAGTTTTTTAga                                                                         >  1:1347753/1‑149 (MQ=38)
                                                                              tCCATGAATTCCGGCGCGTCCAAATTTCGCCTTCCGGGGGGGTCACAAGCGGCCGTTAAGGTGGGTCAAAAGTTGCCGTTCAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGcc                                                                        <  1:125523/150‑1 (MQ=25)
                                                                                cATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTTAGGGGGGGGAAATTTTGCCGCGTAAGAATTGTACTATTAAGATTTCTTGGCGGCGTGAGTTTTTTAgggtt                                                                      >  1:448256/1‑146 (MQ=21)
                                                                                 aTGAATTCAGGTGCGTCAAAATATGCCGTTCCGGGGGGGTCAAAAGCTGCCTTTTAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTcc                                                                     <  1:824589/150‑1 (MQ=37)
                                                                                    aaTTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGGGGGTGAAAAGTTGGCGAGTAAGAATTCTAGTTTTAAGAATTGTACGGGACGGGAGTTAATTAGCATGCcgg                                                                  >  1:990166/1‑150 (MQ=25)
                                                                                        cAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCCGTTAAGGTGGGGCAAAAGTTGACGAGTAAGAATTCTACAATTAAGAATTCTACGCGACGGGAGTTTTTTAGCATCCcggtcgg                                                              >  1:47729/1‑150 (MQ=37)
                                                                                               gTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATaga                                                       >  1:18545/1‑150 (MQ=255)
                                                                                               gTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATaga                                                       >  1:626594/1‑150 (MQ=255)
                                                                                                 cAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTGAAAAGTTGACGAGTAAGAATTCTACTATTAAGAATTCTAGGCGAAGGGAGTTTTTTAGCCTCCCGGTCGGTCATagaga                                                     >  2:614607/1‑150 (MQ=255)
                                                                                                     aGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCg                                                 >  1:409927/1‑150 (MQ=255)
                                                                                                                gtggtgGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaa                                      >  1:894963/1‑150 (MQ=255)
                                                                                                                  ggtggGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc                                    <  2:278052/150‑1 (MQ=255)
                                                                                                                  ggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc                                    >  1:1050637/1‑150 (MQ=255)
                                                                                                                  ggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc                                    >  2:1330956/1‑150 (MQ=255)
                                                                                                                   gtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACa                                   <  2:626594/150‑1 (MQ=255)
                                                                                                                                tGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTc                      >  2:1053486/1‑150 (MQ=255)
                                                                                                                                    gTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgccg                  >  1:328861/1‑150 (MQ=255)
                                                                                                                                       tatgTGGGTCAAATGTTGGCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCt               <  2:8906/147‑1 (MQ=255)
                                                                                                                                           tGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATa           >  2:1053818/1‑150 (MQ=255)
                                                                                                                                            gggTCAAATGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTATACTTAATGTCTGTGATAGATAGTAGATTATATTTAAACATAAGAATTTTTTCTCAGGATGATAg          <  2:470236/150‑1 (MQ=255)
                                                                                                                                              gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGACTAGAGTTAAACAGAAGAATTATTTCGCCGACTGATAGGt        >  1:1372224/1‑150 (MQ=255)
                                                                                                                                              gTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGCCGACTAAGAATTCTACTATTAAGAATTCTAGCCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGACTACAGTTAAACAGAAGAATTAATTCGCCGCCTGATacgg        >  1:170511/1‑147 (MQ=255)
                                                                                                                                               tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTc       <  1:427/150‑1 (MQ=255)
                                                                                                                                                cAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCg      >  2:387321/1‑150 (MQ=255)
                                                                                                                                                 aaaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATCCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCgg     >  2:751701/1‑150 (MQ=255)
                                                                                                                                                  aaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGc    <  2:1050637/150‑1 (MQ=255)
                                                                                                                                                  aaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAG‑‑TTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTg  <  1:1182145/150‑1 (MQ=255)
                                                                                                                                                   aaGTTGCCGTTCAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGGCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCt   <  2:519358/150‑1 (MQ=255)
                                                                                                                                                    aGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATa           >  2:398244/1‑141 (MQ=255)
                                                                                                                                                    aGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATa           <  1:398244/141‑1 (MQ=255)
                                                                                                                                                    aGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTg  <  2:399164/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTG  >  NC_000913/2560864‑2561161

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: