Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,326 A→G 15.4% T24A (ACA→GCA)  yffL → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3260AG15.4% 443.1 / 64.7 182T24A (ACA→GCA) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (71/83);  new base G (16/12);  total (87/95)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 184 total aligned reads displayed.

TAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAT  >  NC_000913/2560227‑2560422
                                                                                                   |                                                                                                
tAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                                                                                                  >  1:3254847/1‑100 (MQ=255)
  aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                                <  1:2821546/100‑1 (MQ=255)
     tattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGc                                                                                             <  1:1148066/100‑1 (MQ=255)
         cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                         <  2:1032608/100‑1 (MQ=38)
         cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                         <  1:2438152/100‑1 (MQ=255)
         cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTAGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                         <  2:616555/100‑1 (MQ=38)
          aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                        <  2:1751561/100‑1 (MQ=38)
          aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                                                                                        >  2:1154605/1‑100 (MQ=255)
             atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGct                                                                                     <  2:2582786/100‑2 (MQ=38)
             atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                                                                                     >  2:1661559/1‑100 (MQ=255)
              tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                    >  1:3048698/1‑100 (MQ=255)
              tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                    <  2:302222/100‑1 (MQ=255)
              tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                    >  1:1889903/1‑100 (MQ=255)
               caacaaCTTCAATCCCTCTTTGCATCTTATGATCCTACGTCCACAAACAGGATCACCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                    <  2:4080503/99‑1 (MQ=255)
                    aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                              <  2:3458884/100‑1 (MQ=255)
                     tttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                             <  1:3839942/99‑1 (MQ=255)
                     cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                             >  1:804197/1‑100 (MQ=255)
                     cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                             <  1:1614272/100‑1 (MQ=255)
                     cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTGCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                             <  2:1879526/100‑1 (MQ=255)
                      ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                            <  2:366410/100‑1 (MQ=255)
                       tAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaaa                                                                           >  2:722641/1‑100 (MQ=255)
                       naaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaaa                                                                           >  2:721706/2‑100 (MQ=255)
                        aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAc                                                                          <  1:1742438/100‑1 (MQ=255)
                         aaTCCATGTTTGCAACGAGAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACa                                                                         <  2:2079152/100‑1 (MQ=255)
                           tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAg                                                                       >  2:2293038/1‑100 (MQ=255)
                           tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAg                                                                       >  1:30616/1‑100 (MQ=255)
                           tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAg                                                                       >  1:3484669/1‑100 (MQ=255)
                            ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                      <  1:155624/100‑1 (MQ=255)
                            ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                      <  2:3396432/100‑1 (MQ=255)
                            ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                      <  2:3071805/100‑1 (MQ=255)
                             cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACGAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                      >  1:2380008/1‑99 (MQ=255)
                             cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCGACATGCCGGGGGCGGTATGAGAAACAAGCt                                                                     >  1:210624/1‑100 (MQ=255)
                               tGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaa                                                                   <  1:3964255/100‑1 (MQ=255)
                                gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaca                                                                  >  2:2514571/1‑98 (MQ=255)
                                gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                                  >  2:2394913/1‑100 (MQ=255)
                                gTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTaaa                                                                  >  1:15971/1‑100 (MQ=255)
                                gTTTGCAACGAAAGATCCAGAGCTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCAGGGGCGGTATGAGAAACAAGCTaa                                                                   >  2:849528/1‑99 (MQ=255)
                                 tttGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                                 >  1:708555/1‑100 (MQ=255)
                                  ttGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                                 <  2:90818/99‑1 (MQ=255)
                                   tGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAACAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGga                                                               <  2:3137638/100‑1 (MQ=255)
                                   tGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGga                                                               >  2:2370279/1‑100 (MQ=255)
                                      aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                            >  1:3650712/1‑100 (MQ=255)
                                      aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                            <  1:2419027/100‑1 (MQ=255)
                                       aCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAt                                                           <  1:4031216/100‑1 (MQ=255)
                                       aCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGCATGAGAAACAAGCTAAAGGAGAAt                                                           >  2:1435061/1‑100 (MQ=255)
                                           aaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                       >  2:499818/1‑100 (MQ=255)
                                           aaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                       <  2:1933166/100‑1 (MQ=255)
                                            cttccgatctttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                      <  1:1129350/90‑1 (MQ=255)
                                            aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                      <  2:2443025/100‑1 (MQ=255)
                                            aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                      <  1:1504573/100‑1 (MQ=255)
                                             tccgatctgTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                     <  1:2565936/92‑1 (MQ=255)
                                              aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                     >  2:3087978/1‑99 (MQ=255)
                                                 cAGAGTTCGAAAATAGGATCAACACAGATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                                 >  1:2512211/1‑100 (MQ=255)
                                                 cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                                 <  2:752120/100‑1 (MQ=255)
                                                 cAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                                 >  1:1380201/1‑100 (MQ=255)
                                                  agagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCt                                                >  1:2416671/1‑100 (MQ=255)
                                                     gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                             >  2:3271259/1‑100 (MQ=255)
                                                     gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                             <  2:2151890/100‑1 (MQ=255)
                                                     gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                             <  2:1577294/100‑1 (MQ=255)
                                                     gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGan                                             <  2:115950/100‑2 (MQ=255)
                                                      ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGata                                            <  1:2487140/100‑1 (MQ=255)
                                                         gAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                         >  1:2683777/1‑100 (MQ=255)
                                                            aaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                                      <  1:311435/100‑1 (MQ=255)
                                                             aTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGc                                     <  1:2947745/100‑1 (MQ=255)
                                                             aTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGc                                     <  2:2803588/100‑1 (MQ=255)
                                                                 gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                                 <  2:1543125/100‑1 (MQ=255)
                                                                     aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                             >  1:2449059/1‑100 (MQ=255)
                                                                     aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                             <  1:1952052/100‑1 (MQ=255)
                                                                     aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                             >  2:2932807/1‑100 (MQ=255)
                                                                      acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                            >  2:1593206/1‑100 (MQ=255)
                                                                       cacaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                           >  2:329518/1‑100 (MQ=255)
                                                                        acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                          >  1:2117365/1‑100 (MQ=255)
                                                                          aaaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaa                        <  2:3015600/100‑1 (MQ=255)
                                                                           aaTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAAAAGCTAAAGGAGAATTTCTCATGTCTGATATGATAGCAGTAGAACAAGaaa                       >  1:2108570/1‑100 (MQ=255)
                                                                            acaaaaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                      <  2:2505383/98‑1 (MQ=255)
                                                                            aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                      <  2:3769489/100‑1 (MQ=255)
                                                                             ataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaaca                     >  2:3829602/6‑100 (MQ=25)
                                                                                ccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAACCAAGCTAAAGGAGAATTTCTCCTGTCTGATATGCTAGCAGTAGAACTAGaaacaaac                  >  2:2726232/3‑100 (MQ=25)
                                                                                ccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac                  >  1:2495375/3‑100 (MQ=37)
                                                                                caacccccGGAATACTGCAACATGCCGCTGGCGGTATGAGAAGCAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGTCGTAGAACAAGaatcaacc                  <  2:3571888/97‑3 (MQ=255)
                                                                                aaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac                  <  1:1393148/100‑1 (MQ=255)
                                                                                 aGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca                 <  2:804191/100‑1 (MQ=255)
                                                                                 aGCCCCCGGAATGCTGCAAAATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGGATTTCTCATATCTGATATGCTAGCAGTAGAACAATAaacaaaca                 >  2:1370733/1‑100 (MQ=255)
                                                                                  gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                >  1:399479/1‑100 (MQ=38)
                                                                                  gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                >  1:458121/1‑100 (MQ=38)
                                                                                  gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTACCAGTAGAACAAGAaataaaca                 >  2:3982801/1‑99 (MQ=37)
                                                                                  gCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                >  1:3883326/1‑100 (MQ=255)
                                                                                   cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAaaaaacaa                >  1:1464475/1‑99 (MQ=37)
                                                                                   cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat               >  1:1048252/1‑100 (MQ=38)
                                                                                   cccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat               >  1:147861/1‑100 (MQ=255)
                                                                                   cccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat               >  1:1273934/1‑100 (MQ=255)
                                                                                     cccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga             >  2:562040/1‑100 (MQ=38)
                                                                                     cccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga             >  1:2319239/1‑100 (MQ=38)
                                                                                     cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatga             <  1:3198629/100‑1 (MQ=255)
                                                                                        ggAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGt          >  1:1843769/1‑100 (MQ=255)
                                                                                         gAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGtt         >  1:1739783/1‑100 (MQ=255)
                                                                                             gctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTc     <  2:708549/100‑1 (MQ=255)
                                                                                              ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCg       <  2:1594883/97‑1 (MQ=255)
                                                                                              ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCa    <  2:3844955/100‑1 (MQ=255)
                                                                                               tgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCaa   >  1:2355530/1‑100 (MQ=255)
                                                                                                   |                                                                                                
TAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAT  >  NC_000913/2560227‑2560422

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