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% 139.4 / 8.5 69intergenic (+87/‑383)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (24/32);  new base A (7/6);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAA  >  NC_000913/2560841‑2561131
                                                                                                                                                |                                                                                                                                                  
ggCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCaaa                                                                                                                                                 <  2:12182/148‑1 (MQ=255)
   aGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGCATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGGCAAAAGtt                                                                                                                                             <  5:172412/149‑1 (MQ=255)
   aGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGtt                                                                                                                                             <  1:158954/149‑1 (MQ=255)
    gAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTg                                                                                                                                            >  2:291450/1‑149 (MQ=255)
      aaGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGcc                                                                                                                                          <  3:27094/149‑1 (MQ=255)
      aaGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGc                                                                                                                                           <  5:449743/148‑1 (MQ=255)
         tgAGTGCAACGAATACCAGTTATTATTCAAGCCGTACTACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTTGGTCAAAAGTTGCCGtt                                                                                                                                       <  8:191136/148‑1 (MQ=255)
         ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGt                                                                                                                                        >  1:184206/1‑148 (MQ=255)
          gAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTCGGACAAATGTAGCCCtt                                                                                                                                       >  5:105967/1‑148 (MQ=255)
           aGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCaaaa                                                                                                                                                <  4:292251/138‑1 (MQ=255)
           aGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCaaaa                                                                                                                                                >  3:292251/1‑138 (MQ=255)
           aGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                     <  3:436810/149‑1 (MQ=255)
           aGTGCAACGAATACCAGGTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTGGCCGGTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                     <  6:442644/149‑1 (MQ=255)
               cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                  <  7:443965/148‑1 (MQ=255)
               cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTg                                                                                                                                 <  1:291450/149‑1 (MQ=255)
                aaCTAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                                <  8:132321/149‑1 (MQ=255)
                aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                                <  3:461910/149‑1 (MQ=255)
                 aCGAATACCACTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGAGTGCTGGGACAAAAGTTGCCGTGAAGGTggg                                                                                                                               >  6:349977/1‑149 (MQ=255)
                  cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCTGGTGGGTCAAAAGTTGCGGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                              <  1:208987/149‑1 (MQ=255)
                  cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCATTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                              >  5:377002/1‑149 (MQ=255)
                       aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGAGTGTGaaaa                                                                                                                         >  1:47395/1‑149 (MQ=255)
                        ccAGTTATTATTCTAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                         <  5:349977/148‑1 (MQ=255)
                        ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaatg                                                                                                                        >  1:189760/1‑147 (MQ=255)
                         cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGATGCCGTAAATTTGGTTCAAAAGt                                                                                                                       >  6:457711/1‑149 (MQ=255)
                          aGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGTTGTGTCAAAATTTGCCTTTAAGGATGAGCAAAAGtt                                                                                                                      >  4:205385/1‑149 (MQ=255)
                          aGTTATTATTCAAGCCGTACAACAAGAAATAACGGAATCAGTCGGTCAATTTTCCAGCCTTTTCTTTGCCTTAATTCCATGAATTAAGGTGGGTCATAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                                      <  6:377002/149‑1 (MQ=255)
                           gTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaattttg                                                                                                                     >  5:387878/1‑144 (MQ=255)
                           gTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                                                                                     >  7:435687/1‑149 (MQ=255)
                            cttccgttCTAGCCTTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTGTTTTCCTTAATTGCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                                                                                                    <  5:13281/143‑1 (MQ=255)
                              attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTAAAAATTTGACGTTAAGGAGGGTCAAAATTTGACg                                                                                                                  >  1:376661/1‑149 (MQ=255)
                                  ttCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCTTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATGGGTGGGTAAAAAGTTGCCGTAATGTTGGGTCAAAAGTTCCCTTTat                                                                                                              >  7:478539/1‑148 (MQ=255)
                                  ttCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATATAGTTGGGTAAAAAGTTGCAGTTTAGGAGGAGCAAAAGTTGGCTTTaa                                                                                                              >  1:170655/1‑149 (MQ=255)
                                   tCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGTGTAAAAAGTTGCCGTGTAGGAGTGTCAAATGTTGACGTTAAg                                                                                                             >  4:144575/1‑149 (MQ=255)
                                    cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAg                                                                                                             <  7:141327/148‑1 (MQ=255)
                                    cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTACATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGATGCCGTTAAGGTGGGTGAAAAATTTCCGTTTAgg                                                                                                            >  3:27098/1‑149 (MQ=255)
                                     aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGAGGGTCAAAATTTTCAGATAAGGt                                                                                                           >  5:490781/1‑149 (MQ=255)
                                     aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAgg                                                                                                            <  1:329453/148‑1 (MQ=255)
                                     aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGACAAAAGTTACCGTTAAGGt                                                                                                           >  1:11734/1‑149 (MQ=255)
                                       gCCGTACAAGATGAAATAACGGAATAATGCGTGCAATTTTCCAGCCTTTTCTTTTCCTTAATTCCATGAATTCAGGTGGGTAAAAAGTTGCCGTTAGTTGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                         <  3:100723/149‑1 (MQ=37)
                                       gCCGTACAACATGAATTAACGGAATCAGGCTTGCATTTTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAATGTTGCCGTTAGTGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                         <  2:189760/149‑1 (MQ=38)
                                       gCCGTACAACATGAAATATCGGATTCAGGCGGGCAATGTTCACGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGGGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                           <  5:457711/147‑1 (MQ=39)
                                         cGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCTTTAAGGAGGGTCAAAAGTTGACGTTAAggggggg                                                                                                       >  8:165610/1‑148 (MQ=255)
                                                caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACGTTAATGTGTGTCaaaatt                                                                                                >  7:237181/1‑147 (MQ=255)
                                                 aTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGGCAAACTTGGCCGTTATGGTGTCTCAAAAAGGTCGGTTGTTGTGGGTCAAAAggt                                                                                               >  5:59544/1‑147 (MQ=255)
                                                 aTGAAATAACGGAATCAGGCGGACAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAATTTTGGCTTTAAGGTGGGGGTAAAGtt                                                                                               >  7:411241/1‑149 (MQ=38)
                                                         aCGGACTCCGTCGGGCAATGTTCCAGCCTTTGCTTTTCCTTAAATACAGGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGCGGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGt                                                                                        <  3:418113/148‑1 (MQ=37)
                                                         aCGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGAGATTGGTGGGTAAAATTTTGCAGTTAAGGAGTGGAAAAATTTTAAGATTAGGTGTGTGAAAAGTTGACGAGt                                                                                        >  8:232456/1‑148 (MQ=255)
                                                           ggAATCAGTCGGGCAAGGTTCCGTGATTTTTTTTGGCTTACTTCCATGCATTCAGGTGGGTCAATAGGTGGCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                       <  5:319151/131‑1 (MQ=21)
                                                                cAGTCGGCCAATCTTACCGCCTTCTCTTTGCCTTAATTACCTGAATTCTTGGGGGTCCAAATTAGGCCTTTCTGGTGGGTAAATTCTTGCCGTTAAGGTGGGTCAAAAGTTTCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTc                                                                                <  2:11734/149‑1 (MQ=16)
                                                                    agggCATTGTTCCCGACTTTTCTTTTCCTTAATTCAATGAATTCAGGTGTGTAAAAAGTTGCCTTTATTGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctac                                                                             <  7:206197/147‑1 (MQ=25)
                                                                                      tttCTTTGCCTTAATTCCATGAATTAAGGTGGGTCAAAAGTTGCCGTTAATGTTGTGTAAAAATTTGCCTTTAAGGAGGGTCAAAATTTGAAGTGTAGGTGGCTGTATATTTGAAGAGAAAGGATTCGGCTATTATTATTTCTAGGc                                                            >  4:282515/1‑147 (MQ=12)
                                                                                                           caTTCAGGTGGGTCTAAATTCGGCGTTCCGTGTGGGTAACTCCCTGCCGTTTCGTTGGGTCAAAAGTTGCCGTTCAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTcc                                       <  5:442818/146‑1 (MQ=21)
                                                                                                           aaTTCAGGGGGCTCAAAAGATGCCGTTACTGTTGGGTCAAAAGCTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcg                                     <  3:205385/149‑1 (MQ=37)
                                                                                                                aggggtgTCAACACTTGCCTTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcgg                                <  3:282515/143‑1 (MQ=37)
                                                                                                                  gTGGGTCAAAAGAAGCCATTACTGGTGGGTAAAAAGTAGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcg                                 <  7:232456/146‑1 (MQ=37)
                                                                                                                       tCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCAGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTACCGGTCGGTCATag                          >  3:331298/1‑148 (MQ=255)
                                                                                                                             gTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTATGGTGGGTCAAAAGTTGACTTGTAAGAATTCTACTATTTTGAATTATAGGCGAGGGGAGTTTTTTa                                             >  6:246059/1‑123 (MQ=25)
                                                                                                                                 cctttATTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa               <  8:435687/146‑1 (MQ=255)
                                                                                                                                     tcagggggggTAAAATGTAGCCGTTACGGTGGGTCAAAAGTTGCATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga             <  5:176432/140‑1 (MQ=38)
                                                                                                                                      tgtggtgTGTCAATAGCTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGt           <  7:185769/147‑1 (MQ=255)
                                                                                                                                      aggggtggCTAAAAAGCTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGt           <  7:165610/145‑1 (MQ=38)
                                                                                                                                       gtggtgGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTACCGTTAAGGTGGGTCAAAAGTTTCGGAGAAAGAATTCTACTATTAAGAATTCTAGGCGACGGGACTTTTTTAGACTCCCGGGCGGTGGTAGAGAGTCGGCTAGAGt           >  2:83091/1‑147 (MQ=25)
                                                                                                                                          gtggGTCAAAAGCTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaaa       <  1:83091/148‑1 (MQ=255)
                                                                                                                                              ggcAAACGCTGCCGTTCAGGTGGGTCAATAGATGCCGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaa  <  1:501449/147‑1 (MQ=255)
                                                                                                                                              gTCAACATTTGCCGTTACGGTGGGTCAAACGTTGCCGTTAAGGTGTGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaa  <  3:48973/149‑1 (MQ=255)
                                                                                                                                               tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGcc                 <  7:363160/133‑1 (MQ=255)
                                                                                                                                               tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGcc                 >  8:363160/1‑133 (MQ=255)
                                                                                                                                               tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaaa       >  6:249115/1‑143 (MQ=255)
                                                                                                                                               tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaaa       <  5:249115/143‑1 (MQ=255)
                                                                                                                                               tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag    >  2:358732/1‑146 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAA  >  NC_000913/2560841‑2561131

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: