Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3260AG23.1% 284.4 / 88.3 148T24A (ACA→GCA) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (45/68);  new base G (20/14);  total (66/82)
Fisher's exact test for biased strand distribution p-value = 7.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

Only 100 of 153 total aligned reads displayed.

AAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTT  >  NC_000913/2560228‑2560425
                                                                                                  |                                                                                                    
aaaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAAc                                                                                                     <  1:101980/100‑1 (MQ=255)
 aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                                    <  2:2378843/100‑1 (MQ=255)
 aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                                                                                                    >  2:2561687/1‑100 (MQ=255)
  attattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACAt                                                                                                   <  1:1383540/100‑1 (MQ=255)
       tCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgg                                                                                             <  1:12822/100‑1 (MQ=38)
       tCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgg                                                                                             <  1:2898111/100‑1 (MQ=38)
        cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                            <  2:2473428/100‑1 (MQ=38)
        cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCACCATGcc                                                                                                >  2:636395/1‑97 (MQ=255)
        cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                            >  2:1580723/1‑100 (MQ=255)
         aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                           <  1:1066831/100‑1 (MQ=38)
          atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                          <  1:3819/100‑1 (MQ=38)
          atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                          <  1:2446359/100‑1 (MQ=38)
          atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCAAGAGTTCGAAACTAGGATCAACACACATCAAAGCCCCAGGCCTGTTGCCGCATTCCggggg                                                                                          <  2:446551/100‑1 (MQ=255)
           tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCGAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                         >  2:1637818/1‑100 (MQ=255)
           tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGc                                                                                         <  2:7437/100‑1 (MQ=38)
           tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                                                                                         <  1:1264433/100‑1 (MQ=255)
            atCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                                                                                        <  1:984968/100‑1 (MQ=255)
             tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                       <  1:2599976/100‑1 (MQ=255)
             tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                       <  1:2306916/100‑1 (MQ=255)
              cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGt                                                                                      <  1:982861/100‑1 (MQ=255)
              cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGt                                                                                      >  1:154583/1‑100 (MQ=255)
                aGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTAt                                                                                    >  2:728494/1‑100 (MQ=255)
                  aaCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                                  <  1:909262/100‑1 (MQ=255)
                   aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                                 <  1:308232/100‑1 (MQ=255)
                   aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgag                                                                                 >  1:2443024/1‑100 (MQ=255)
                    cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGcttttatttg                                                                                <  1:459154/100‑10 (MQ=21)
                    cTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATgaga                                                                                >  1:524184/1‑100 (MQ=255)
                     ttAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa                                                                               <  2:1545081/100‑1 (MQ=255)
                     ttAAATCCATGTTTGCAACGAAAAATCCAGAGTTGGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCGGGGGcttttatt                                                                                  <  2:570940/97‑8 (MQ=21)
                      tAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaaa                                                                              >  2:2863172/1‑100 (MQ=255)
                      tAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaaa                                                                              >  1:2321110/1‑100 (MQ=255)
                       aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAc                                                                             <  1:37494/100‑1 (MQ=255)
                       aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAc                                                                             <  1:688556/100‑1 (MQ=255)
                        aaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACa                                                                            <  1:1427368/100‑1 (MQ=255)
                        aaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACa                                                                            >  2:1142663/1‑100 (MQ=255)
                        aaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACa                                                                            >  2:525650/1‑100 (MQ=255)
                         aTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                           <  1:799409/100‑1 (MQ=255)
                           ccATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGc                                                                         <  1:17424/100‑1 (MQ=255)
                            cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCt                                                                        >  2:1288748/1‑100 (MQ=255)
                            cATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCt                                                                        <  2:1676575/100‑1 (MQ=255)
                                  tGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGga                                                                  <  1:613566/100‑1 (MQ=255)
                                    cAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgaga                                                                >  1:2740418/1‑100 (MQ=255)
                                      aCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAt                                                              <  1:1030381/100‑1 (MQ=255)
                                          aaGATCCAGAGTGCGAAAATAGGATCAACAGAAATAAAAGCCCCCGGAATGCTGCAACATGTCGGGGGCGGTATGAGTAACAAGCTAAAGGAGAAGTtct                                                          <  2:1334453/100‑1 (MQ=255)
                                             aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                                       >  1:663210/1‑100 (MQ=255)
                                             aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAt                                                       <  2:1767470/100‑1 (MQ=255)
                                                  gagTTCGAAAATAGGATCAACACAAATAAAAGCCCCGGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTg                                                  <  2:2026890/100‑1 (MQ=255)
                                                  gagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTg                                                  <  1:1974944/100‑1 (MQ=255)
                                                   agTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGa                                                 <  1:1043454/100‑1 (MQ=255)
                                                    gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                                >  2:317133/1‑100 (MQ=255)
                                                    gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                                <  2:353117/100‑1 (MQ=255)
                                                    gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                                <  2:101755/100‑1 (MQ=255)
                                                     ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAAGATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGata                                               >  2:1699925/1‑100 (MQ=255)
                                                     ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGata                                               >  2:3055871/1‑100 (MQ=255)
                                                     ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGata                                               >  1:2499233/1‑100 (MQ=255)
                                                      tCGAAAATAGTATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGATACGAGCTAAAGGAGAATTTCTCATGTCTGatat                                              <  2:2000896/100‑1 (MQ=255)
                                                       cGAGAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGata                                               <  2:2576818/98‑1 (MQ=255)
                                                         aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                           >  2:226284/1‑100 (MQ=255)
                                                         aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                           <  1:1288760/100‑1 (MQ=255)
                                                         aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCt                                           <  1:906351/100‑1 (MQ=255)
                                                            aTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGc                                        <  1:1588043/100‑1 (MQ=255)
                                                            aTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGc                                        <  2:814999/100‑1 (MQ=255)
                                                              aGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAg                                      <  1:2653873/100‑1 (MQ=255)
                                                               ggATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGt                                     >  1:1389421/1‑100 (MQ=255)
                                                                gATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                                    <  1:2317591/100‑1 (MQ=255)
                                                                   cAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGaa                                 >  1:2097124/1‑100 (MQ=255)
                                                                   cAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGaa                                 <  2:799325/100‑1 (MQ=255)
                                                                   cAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGaa                                 >  1:2128862/1‑100 (MQ=255)
                                                                    aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                                >  1:1331704/1‑100 (MQ=255)
                                                                    aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                                >  2:2717019/1‑100 (MQ=255)
                                                                    aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                                >  1:2218512/1‑100 (MQ=255)
                                                                    aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                                >  1:3052624/1‑100 (MQ=255)
                                                                     acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                               <  1:1590696/100‑1 (MQ=255)
                                                                     acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACa                               <  1:2408899/100‑1 (MQ=255)
                                                                       ttctcataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAg                             >  1:72300/11‑100 (MQ=21)
                                                                        caAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGa                            <  1:780564/100‑1 (MQ=255)
                                                                         aaaTAAAAGCCCCGGAATGCTGCAACATGCCGGGGGGCGGTATGAGAAACAAGCTAAAGGATAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaa                           <  1:2151522/100‑1 (MQ=255)
                                                                           aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                         <  1:2820358/100‑1 (MQ=255)
                                                                            ataccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACCAGaaaca                        >  1:1804325/6‑100 (MQ=25)
                                                                             aaaaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaa                       <  1:2717056/100‑1 (MQ=255)
                                                                               ccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac                     >  2:65720/3‑100 (MQ=37)
                                                                               ccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac                     >  2:2754615/3‑100 (MQ=37)
                                                                                cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca                    >  1:2414832/2‑100 (MQ=38)
                                                                                cgCCCCCGGCATGTTGCAGCATGCCGGGGGCGGTATGAGAAACAAGCGAAAGGAGAATTTATCATGTCTGATATGATAGCAGTAGAACAAGaaacaaaaa                    >  2:2132015/2‑98 (MQ=37)
                                                                                aGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca                    <  2:2897668/100‑1 (MQ=255)
                                                                                aGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca                    >  2:2218282/1‑100 (MQ=255)
                                                                                 gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                   >  1:1189134/1‑100 (MQ=38)
                                                                                 gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa                   >  1:1085659/1‑100 (MQ=38)
                                                                                  cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat                  >  2:2842370/1‑100 (MQ=38)
                                                                                  cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTATGATATGCTAGCAGTAGAACAAGAAacaaacaa                   >  1:2280778/1‑99 (MQ=38)
                                                                                  cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCAGGTCTGGTATGCTAGGAGTAGAACAAGAAACAAACAat                  >  2:2053283/1‑100 (MQ=37)
                                                                                      cGGAATGCTGCAACATGCCGGGGGCGGTATGAGATACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgatg              <  2:1614269/100‑1 (MQ=255)
                                                                                          aTGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCg          <  1:2037685/100‑1 (MQ=255)
                                                                                             ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTCAAGGAGAATTTCTCATGTCTGATATGCTAGCAGAAGAACAAGAAACAAACAATGATGTTCGTCa       >  2:890406/1‑100 (MQ=255)
                                                                                             ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCa       >  1:544327/1‑100 (MQ=255)
                                                                                              tgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCaa      >  1:1879806/1‑100 (MQ=255)
                                                                                              tgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCaa      <  2:2442992/100‑1 (MQ=255)
                                                                                                  aCATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAtttt  >  2:2390593/1‑100 (MQ=255)
                                                                                                  aCATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAtttt  >  2:2176977/1‑100 (MQ=255)
                                                                                                  aCATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAAtt    >  2:2161170/1‑98 (MQ=255)
                                                                                                  |                                                                                                    
AAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTCAATTTT  >  NC_000913/2560228‑2560425

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