Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,560,085 (T)6→7 intergenic (+19/‑172) intZ → / → yffL CPZ‑55 prophage; putative phage integrase IntZ/CPZ‑55 prophage; uncharacterized protein YffL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,0791.T98.4% 213.7 / ‑3.6 63intergenic (+13/‑178)intZ/yffLCPZ‑55 prophage; putative phage integrase IntZ/CPZ‑55 prophage; uncharacterized protein YffL
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (28/34);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTA  >  NC_000913/2559940‑2560210
                                                                                                                                            |                                                                                                                                    
gAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGc                                                                                                                                       <  2:646474/140‑1 (MQ=255)
       aTAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCttttttt                                                                                                                              >  2:1059572/1‑139 (MQ=255)
       aTAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCttttttt                                                                                                                              >  1:425301/1‑139 (MQ=255)
        tAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCtttttttg                                                                                                                             >  1:1607926/1‑138 (MQ=255)
              gCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAgc                                                                                                                        >  1:898460/1‑140 (MQ=255)
               cGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAgcg                                                                                                                       <  2:267375/140‑1 (MQ=255)
               cGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAgcg                                                                                                                       <  2:350588/140‑1 (MQ=255)
                gTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAgcgc                                                                                                                      >  1:90017/1‑140 (MQ=255)
                      aaCCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAAtttt                                                                                                                <  1:460107/140‑1 (MQ=255)
                         cAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCa                                                                                                             <  2:1503686/140‑1 (MQ=255)
                         cAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCa                                                                                                             <  2:538116/140‑1 (MQ=255)
                         cAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCa                                                                                                             <  1:25744/140‑1 (MQ=255)
                         cAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACAGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCCATTTTGCa                                                                                                             <  2:1419469/140‑1 (MQ=255)
                         cAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGAAGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCa                                                                                                             <  2:49567/140‑1 (MQ=255)
                          aGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCAt                                                                                                            <  2:915185/140‑1 (MQ=255)
                           ggCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCAtt                                                                                                           <  1:1604553/140‑1 (MQ=255)
                           ggCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCAtt                                                                                                           >  2:501919/1‑140 (MQ=255)
                           ggCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCAtt                                                                                                           <  1:1253929/140‑1 (MQ=255)
                            gCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCAttt                                                                                                          >  2:666178/1‑140 (MQ=255)
                            gCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCAttt                                                                                                          >  2:513503/1‑140 (MQ=255)
                            gCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCAttt                                                                                                          >  1:277583/1‑140 (MQ=255)
                            gCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAAAGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCAttt                                                                                                          >  2:1548382/1‑140 (MQ=255)
                            gCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAAACGGGATGCAAAATTCAGCGGCTTTTTTTGCAAAGCGCAATTTTGCAttt                                                                                                          >  2:689416/1‑140 (MQ=255)
                              tGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTg                                                                                                        >  2:653100/1‑140 (MQ=255)
                                ccTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGta                                                                                                      >  1:1246111/1‑140 (MQ=255)
                                  tACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGtata                                                                                                    <  1:1143722/140‑1 (MQ=255)
                                          gCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTAtttg                                                                                            >  1:788580/1‑140 (MQ=255)
                                          gCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGAAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTAtttg                                                                                            >  1:915098/1‑140 (MQ=255)
                                              cGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTAtttgtttt                                                                                        >  1:1196336/1‑140 (MQ=255)
                                                 cGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTtgtg                                                                                     <  2:768210/140‑1 (MQ=255)
                                                 cGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTtgtg                                                                                     <  2:982569/140‑1 (MQ=255)
                                                 cGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTtgtg                                                                                     <  2:702248/140‑1 (MQ=255)
                                                 cGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTtgtg                                                                                     <  1:682633/140‑1 (MQ=255)
                                                 cGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTtgtg                                                                                     <  1:782106/140‑1 (MQ=255)
                                                  gAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTtgtgt                                                                                    <  2:842767/140‑1 (MQ=255)
                                                   aGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTtgtgtg                                                                                   >  2:1611365/1‑140 (MQ=255)
                                                      aaTGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGtt                                                                                <  1:510563/140‑1 (MQ=255)
                                                      aaTGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGtt                                                                                <  1:653100/140‑1 (MQ=255)
                                                        tgatgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTa                                                                              <  1:133643/140‑1 (MQ=255)
                                                         gatgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTaa                                                                             <  1:1059572/140‑1 (MQ=255)
                                                         gatgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTaa                                                                             <  2:1364015/140‑1 (MQ=255)
                                                         gatgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTaa                                                                             <  2:1607926/140‑1 (MQ=255)
                                                          atgCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAAc                                                                            <  1:1549592/140‑1 (MQ=255)
                                                                  ttACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCt                                                                    >  2:924827/1‑140 (MQ=255)
                                                                  ttACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCt                                                                    <  1:708026/140‑1 (MQ=255)
                                                                    aCGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGa                                                                  >  1:730848/1‑140 (MQ=255)
                                                                          aCGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGaa                                                            >  2:1430186/1‑140 (MQ=255)
                                                                          aCGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGaa                                                            >  1:331016/1‑140 (MQ=255)
                                                                                             ggCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCg                                         >  2:1335131/1‑140 (MQ=255)
                                                                                                   tAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCg                                   >  2:1132495/1‑140 (MQ=255)
                                                                                                   tAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCg                                   >  1:78900/1‑140 (MQ=255)
                                                                                                   tAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGGATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCg                                   >  2:380185/1‑140 (MQ=255)
                                                                                                            tCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTTTATTCAATTTGCCGGCTTCGTCCGGcttt                          >  1:991978/1‑140 (MQ=255)
                                                                                                            tCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGcttt                          >  1:737404/1‑140 (MQ=255)
                                                                                                                     aCCGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCg                 <  2:157301/140‑1 (MQ=255)
                                                                                                                     aCCGGAATGAATAATTAAGCGGCTTTTTGTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCg                 <  1:689416/140‑1 (MQ=255)
                                                                                                                      ccGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGc                >  1:1042011/1‑140 (MQ=255)
                                                                                                                       cGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCt               <  1:1611365/140‑1 (MQ=255)
                                                                                                                       cGGAATGAATAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCt               <  2:1154078/140‑1 (MQ=255)
                                                                                                                         gaatgaatAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTaa             >  2:1306622/1‑140 (MQ=255)
                                                                                                                         gaatgaatAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTACATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTaa             >  1:1341899/1‑140 (MQ=255)
                                                                                                                            tgaatAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAatta          <  2:317426/140‑1 (MQ=255)
                                                                                                                            tgaatAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAatta          <  1:134561/140‑1 (MQ=255)
                                                                                                                             gaatAATTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAattat         <  2:115028/140‑1 (MQ=255)
                                                                                                                              aataatTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAattatt        >  2:341212/1‑140 (MQ=255)
                                                                                                                                 aatTAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAattattcca     <  2:1196336/140‑1 (MQ=255)
                                                                                                                                    tAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAattattccatta  <  1:501919/140‑1 (MQ=255)
                                                                                                                                    tAAGCGGCTTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAattattccatta  >  1:945443/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                    
GAAGGGAATAAGGGGCGTTTATAACCAGGCTGCCTACCTGAAGCAACGGCGAGCAATGATGCAGGATTACGCCGACGCTATTGATTCAATATTGGCAGGTAACGGCAATCCTCTAGAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTA  >  NC_000913/2559940‑2560210

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