Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3170AC20.7% 398.0 / 118.4 198N21H (AAT→CAT) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (74/83);  new base C (16/25);  total (90/108)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 202 total aligned reads displayed.

AGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTT  >  NC_000913/2560220‑2560415
                                                                                                 |                                                                                                  
  tGGAATAAATTACTCAATATCAAGCACTTAAATCCATGTGTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCCACACAACTAAAACCCCCCGGactc                                                                                                 >  2:2171317/1‑96 (MQ=255)
   ggAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgc                                                                                                <  2:1430658/99‑1 (MQ=255)
    gAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGACAATAGGATCAACACAAATAAAAGCCCCCGGAAtgct                                                                                               >  1:4238616/1‑99 (MQ=255)
    gAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgct                                                                                               >  1:903164/1‑99 (MQ=255)
    gAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgct                                                                                               >  1:2269726/1‑99 (MQ=255)
      aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                             >  2:3317936/1‑99 (MQ=255)
       tAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                            <  2:3053873/99‑1 (MQ=255)
         aaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAAc                                                                                          <  2:2743779/99‑1 (MQ=255)
                cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgg                                                                                   <  2:1759687/99‑1 (MQ=38)
                cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgg                                                                                   <  1:2528805/99‑1 (MQ=38)
                cAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgg                                                                                   <  1:4006011/99‑1 (MQ=255)
                 aaTATCAAGAACTTAAATCCCTGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCACCCGGACTGCTGCAACATGCCcgg                                                                                  >  2:4169647/1‑99 (MQ=255)
                 aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                  <  2:2444139/99‑1 (MQ=38)
                 aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                  <  1:3450135/99‑1 (MQ=38)
                 aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                  <  2:2148715/99‑1 (MQ=38)
                 aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                  <  1:3358036/99‑1 (MQ=38)
                 aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                  <  2:3070353/99‑1 (MQ=255)
                 aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                  >  1:2708713/1‑99 (MQ=255)
                 aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTGCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                  <  1:1026315/99‑1 (MQ=38)
                  atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                 <  2:4717542/99‑1 (MQ=38)
                  atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                 <  1:4781179/99‑1 (MQ=38)
                  atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                 <  1:197935/99‑1 (MQ=38)
                  atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                 <  1:2067395/99‑1 (MQ=38)
                  atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                 <  2:2612473/99‑1 (MQ=38)
                  atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                 <  2:521281/99‑1 (MQ=38)
                  atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                                                                                 >  2:4513456/1‑99 (MQ=255)
                   tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggggg                                                                                <  1:4053732/99‑1 (MQ=38)
                   tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggggg                                                                                >  2:3268997/1‑99 (MQ=255)
                     tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                                                                              <  1:327801/99‑1 (MQ=255)
                      cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                             <  1:1223003/99‑1 (MQ=255)
                       aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGt                                                                            <  2:270614/99‑1 (MQ=255)
                       aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGt                                                                            >  2:1181654/1‑99 (MQ=255)
                        aGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                           >  1:1813793/1‑99 (MQ=255)
                         gAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTAt                                                                          <  1:2630632/99‑1 (MQ=255)
                           aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                        <  1:4588112/99‑1 (MQ=255)
                                 aTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACa                                                                  >  1:1328520/1‑99 (MQ=255)
                                 aTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCC‑GAAAGCTGCAACATGCCGGGGGGGGGATGAGaaaaaa                                                                 >  1:3543210/1‑99 (MQ=255)
                                 aTCCATGTTTGCAACGAAAGATCCAGAATTCGAAAATAGGATCAACACAAATAAAAGACCCCGGAATGCGGCAACATGCCGGGGGCGGAATGAGAcaca                                                                  >  1:169499/1‑99 (MQ=255)
                                  tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                 >  2:1163910/1‑99 (MQ=255)
                                      tGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTa                                                             <  2:3568689/99‑1 (MQ=255)
                                         ttGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAg                                                          >  1:4178144/1‑99 (MQ=255)
                                             aaCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgaga                                                      >  1:936671/1‑99 (MQ=255)
                                              aCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgagac                                                     >  1:4394568/1‑98 (MQ=255)
                                              aCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGaa                                                     >  1:913190/1‑99 (MQ=255)
                                               cGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGTAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAt                                                    <  1:4672852/99‑1 (MQ=255)
                                                gAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAtt                                                   >  2:4278667/1‑99 (MQ=255)
                                                 aaaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAttt                                                  <  1:1911633/99‑1 (MQ=255)
                                                  aaGAGCCAGAGGTAGAAAATAGGATCAACACAAATAAAAGCCCCAGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtc                                                 <  2:1782933/99‑1 (MQ=255)
                                                   aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                <  1:4190435/99‑1 (MQ=255)
                                                    gATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                               >  1:3064211/1‑99 (MQ=255)
                                                     aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                              >  2:855466/1‑99 (MQ=255)
                                                     aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                              <  2:4007356/99‑1 (MQ=255)
                                                         agagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                          >  1:4334410/1‑99 (MQ=255)
                                                             ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                      >  1:4272127/1‑99 (MQ=255)
                                                             ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                      >  1:3972489/1‑99 (MQ=255)
                                                             ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                      >  2:3346310/1‑99 (MQ=255)
                                                             ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                      <  2:3420282/99‑1 (MQ=255)
                                                               cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGatat                                    >  2:1527446/1‑99 (MQ=255)
                                                               cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGatat                                    <  1:1503657/99‑1 (MQ=255)
                                                                gAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATg                                   >  1:3708959/1‑99 (MQ=255)
                                                                 aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                  >  2:2047096/1‑99 (MQ=255)
                                                                   aaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCC‑GGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                               <  2:2258962/99‑1 (MQ=255)
                                                                   aaTAGGAACAACACAAATAAAAGCCCCCGGAATGCTGAAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTa                                <  2:1871045/99‑1 (MQ=255)
                                                                    aTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                               >  1:2784342/1‑99 (MQ=255)
                                                                    aTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAg                               <  2:3131396/99‑1 (MQ=255)
                                                                    aTAGGATCAACACAAACAAAAGCCCACGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCCAg                               <  2:467836/99‑1 (MQ=255)
                                                                      aGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCa                             <  2:4874982/99‑1 (MQ=255)
                                                                       ggATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAg                            >  2:174085/1‑99 (MQ=255)
                                                                         aTCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                          <  1:1452415/99‑1 (MQ=255)
                                                                            aaCACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGaa                       >  1:1358611/1‑99 (MQ=255)
                                                                             acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                      >  1:3195/1‑99 (MQ=255)
                                                                             acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                      >  1:5146025/1‑99 (MQ=255)
                                                                             acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                      >  1:1523976/1‑99 (MQ=255)
                                                                             acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                      >  1:4999273/1‑99 (MQ=255)
                                                                               acaAATAAAAGCCCCCGTAATGCTGCAACATGCCGCGGGCGGTATGAGAAACAAGCTAAAGGAAAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                    >  2:1597439/1‑99 (MQ=255)
                                                                               acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCCTGTCTGATATGCTAGCAGTAGAACaa                    >  1:482972/1‑99 (MQ=255)
                                                                               acaAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACaa                    >  2:1042774/1‑99 (MQ=255)
                                                                                   aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                <  1:766711/99‑1 (MQ=255)
                                                                                   aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                <  1:540917/99‑1 (MQ=255)
                                                                                    tAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac               <  2:1556880/99‑1 (MQ=255)
                                                                                       aaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaa            <  1:5040640/99‑1 (MQ=255)
                                                                                        cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac           >  1:646742/2‑99 (MQ=37)
                                                                                        cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac           >  1:3585619/2‑99 (MQ=37)
                                                                                         gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca          >  1:236476/1‑99 (MQ=38)
                                                                                         gCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca          >  2:738349/1‑99 (MQ=255)
                                                                                          cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa         >  2:2401023/1‑99 (MQ=38)
                                                                                          cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa         >  1:3633672/1‑99 (MQ=38)
                                                                                          cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa         >  1:3817270/1‑99 (MQ=38)
                                                                                           ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGCAACAAGGTAAAGGAGAATTTCTCATGACTGATATGCCAGCAGCAGAACAAGAAACAAACAat        >  2:4609916/1‑99 (MQ=25)
                                                                                           ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat        >  2:2002041/1‑99 (MQ=38)
                                                                                           ccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat        >  1:2802411/1‑99 (MQ=255)
                                                                                           ccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat        <  2:4343490/99‑1 (MQ=255)
                                                                                           ccccGACATGTTGCAGCCTTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat        >  1:4317106/1‑99 (MQ=255)
                                                                                            cccGGAATGCTGCAACATGCCGGGGGCGGTATGTGAAACAAGCTAAAGGAGAGTTTCTCATCTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg       >  1:4822847/1‑99 (MQ=255)
                                                                                            cccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatg       >  2:3852314/1‑99 (MQ=255)
                                                                                              cGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAATCAatgat     <  2:2998815/99‑1 (MQ=255)
                                                                                              cGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgat     <  2:3783637/99‑1 (MQ=255)
                                                                                               ggAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgatg    >  2:3466945/1‑99 (MQ=255)
                                                                                               ggAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAatgatg    >  1:113659/1‑99 (MQ=255)
                                                                                                 aaTGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGtt  <  1:84017/99‑1 (MQ=255)
                                                                                                 |                                                                                                  
AGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTT  >  NC_000913/2560220‑2560415

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