Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,3210CT21.8% 365.5 / 100.2 189A22V (GCT→GTT) yffLCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (68/79);  new base T (16/25);  total (84/105)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 203 total aligned reads displayed.

GGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTC  >  NC_000913/2560223‑2560419
                                                                                                  |                                                                                                  
ggAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgc                                                                                                    <  2:1430658/99‑1 (MQ=255)
 gAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGACAATAGGATCAACACAAATAAAAGCCCCCGGAAtgct                                                                                                   >  1:4238616/1‑99 (MQ=255)
 gAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgct                                                                                                   >  1:2269726/1‑99 (MQ=255)
  aaTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctg                                                                                                  <  1:3940039/99‑1 (MQ=255)
   aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                                 >  2:3317936/1‑99 (MQ=255)
   aTAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                                                                                                 <  1:733108/99‑1 (MQ=255)
    tAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCa                                                                                                <  2:3053873/99‑1 (MQ=255)
     aaaTTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                                                                                               <  2:1919218/99‑1 (MQ=255)
        ttattCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACAt                                                                                            <  1:823868/99‑1 (MQ=255)
           ttCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTcc                                                                                         <  2:2713141/99‑1 (MQ=255)
            tCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCg                                                                                        >  2:2485382/1‑99 (MQ=38)
              aaTATCAAGAACTTAAATCCCTGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCACCCGGACTGCTGCAACATGCCcgg                                                                                      >  2:4169647/1‑99 (MQ=255)
              aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                      <  2:2148715/99‑1 (MQ=38)
              aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                      <  1:4010739/99‑1 (MQ=38)
              aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                      <  2:2435518/99‑1 (MQ=38)
              aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                      <  1:3358036/99‑1 (MQ=38)
              aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                      <  1:3450135/99‑1 (MQ=38)
              aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggg                                                                                      <  2:3070353/99‑1 (MQ=255)
              aaTATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTGCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCggg                                                                                      <  1:1026315/99‑1 (MQ=38)
               atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                     <  2:521281/99‑1 (MQ=38)
               atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                     <  1:4781179/99‑1 (MQ=38)
               atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                     <  2:4717542/99‑1 (MQ=38)
               atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                     <  1:2067395/99‑1 (MQ=38)
               atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGCATGTTGCAGCATTCCgggg                                                                                     <  1:197935/99‑1 (MQ=38)
               atatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                                                                                     >  2:4513456/1‑99 (MQ=255)
                tatCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCggggg                                                                                    >  2:3268997/1‑99 (MQ=255)
                  tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                                                                                  <  1:39882/99‑1 (MQ=255)
                  tCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                                                                                  <  1:327801/99‑1 (MQ=255)
                   cAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCgg                                                                                 <  1:1223003/99‑1 (MQ=255)
                    aaGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGt                                                                                >  2:1181654/1‑99 (MQ=255)
                     aGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa                                                                               >  1:1813793/1‑99 (MQ=255)
                       aaCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATg                                                                             <  1:1358143/99‑1 (MQ=255)
                        aCTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATga                                                                            <  1:4588112/99‑1 (MQ=255)
                            aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaaa                                                                        <  2:1779285/99‑1 (MQ=255)
                            aaaTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaaa                                                                        <  1:4079886/99‑1 (MQ=255)
                              aTCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGAAAGCTGCAACATGCCGGGGGGGGGATGAGaaaaaa                                                                     >  1:3543210/1‑99 (MQ=255)
                              aTCCATGTTTGCAACGAAAGATCCAGAATTCGAAAATAGGATCAACACAAATAAAAGACCCCGGAATGCGGCAACATGCCGGGGGCGGAATGAGAcaca                                                                      >  1:169499/1‑99 (MQ=255)
                               tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                     >  1:1358634/1‑99 (MQ=255)
                               tCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACaa                                                                     >  2:1163910/1‑99 (MQ=255)
                                   tGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTa                                                                 <  2:3568689/99‑1 (MQ=255)
                                           aCGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGgagac                                                         >  1:4394568/1‑98 (MQ=255)
                                            cGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGTAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAt                                                        <  1:4672852/99‑1 (MQ=255)
                                             gAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAtt                                                       >  2:4278667/1‑99 (MQ=255)
                                              aaaGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAAttt                                                      <  1:1911633/99‑1 (MQ=255)
                                                aGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtct                                                    <  1:4190435/99‑1 (MQ=255)
                                                 gATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTtctc                                                   >  1:3064211/1‑99 (MQ=255)
                                                  aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                  <  2:4007356/99‑1 (MQ=255)
                                                  aTCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCa                                                  >  2:855466/1‑99 (MQ=255)
                                                      agagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                              >  1:4334410/1‑99 (MQ=255)
                                                      agagTGCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTc                                              >  1:918378/1‑99 (MQ=255)
                                                       gagTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCt                                             >  1:900996/1‑99 (MQ=255)
                                                         gTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGa                                           <  1:2616669/99‑1 (MQ=255)
                                                          ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                          <  2:3420282/99‑1 (MQ=255)
                                                          ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                          >  2:3346310/1‑99 (MQ=255)
                                                          ttCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGat                                          >  1:4272127/1‑99 (MQ=255)
                                                            cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGatat                                        <  1:1503657/99‑1 (MQ=255)
                                                            cGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGatat                                        >  2:1527446/1‑99 (MQ=255)
                                                             gAAAATAGGATGAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATg                                       <  1:2786993/99‑1 (MQ=255)
                                                              aaaaTAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGc                                      >  2:2047096/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:2772222/99‑1 (MQ=255)
                                                                      aTCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTa                              <  1:1452415/99‑1 (MQ=255)
                                                                       tCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAg                             <  2:262635/99‑1 (MQ=255)
                                                                         aaCACAAATAAAAGCCCCCGGAATGGTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGaa                           <  1:2688918/99‑1 (MQ=255)
                                                                          acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                          >  1:5146025/1‑99 (MQ=255)
                                                                          acacAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAAc                          >  1:3195/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)
                                                                               aaTAAAAGCTCCCGGAATGCTGCAACATATAGCGGGAGGTATGACAAACAAGCTAAACAAAAATTTCCGAAGTATGGAATGATAGCAGTAGAACAAGaa                     >  2:2238739/1‑99 (MQ=255)
                                                                                aTAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaa                    <  1:540917/99‑1 (MQ=255)
                                                                                 tAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaac                   <  2:1556880/99‑1 (MQ=255)
                                                                                  taccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaaca                  >  1:3788786/5‑99 (MQ=37)
                                                                                  taccgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaaca                  >  2:2335426/5‑99 (MQ=37)
                                                                                    aaGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaa                <  1:5040640/99‑1 (MQ=255)
                                                                                     cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac               >  1:3585619/2‑99 (MQ=37)
                                                                                     cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac               >  1:908946/2‑99 (MQ=37)
                                                                                     cgCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGaaacaaac               >  1:646742/2‑99 (MQ=37)
                                                                                      gCCCCCGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca              >  1:236476/1‑99 (MQ=38)
                                                                                      gCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAaacaaaca              >  2:738349/1‑99 (MQ=255)
                                                                                       cccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa             >  1:3633672/1‑99 (MQ=38)
                                                                                       cccccGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAacaaacaa             <  2:1748568/99‑1 (MQ=255)
                                                                                        ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGCAACAAGGTAAAGGAGAATTTCTCATGACTGATATGCCAGCAGCAGAACAAGAAACAAACAat            >  2:4609916/1‑99 (MQ=25)
                                                                                        ccccGGCATGTTGCAGCATTCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAat            >  1:275981/1‑99 (MQ=38)
                                                                                        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)
                                                                                              aaTGCTGCAACATGCCGGGGGCGGTATGAGAAATAAGCGAAATGAGAATTTCTCATGTCTGATATGCTAGCAATAGAACAAGAAACAAACACTGATGtt      >  2:1861131/1‑99 (MQ=255)
                                                                                               aTGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTc     <  1:3703578/99‑1 (MQ=255)
                                                                                                  ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGCATCTCTCCTGCCTGATATGCTAGCAGTAGAACAAGCACCAAACCATGATGTTCGTc  >  2:1708482/1‑99 (MQ=255)
                                                                                                  ctgcAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTc  >  1:1319398/1‑99 (MQ=255)
                                                                                                  |                                                                                                  
GGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAGCTAAAGGAGAATTTCTCATGTCTGATATGCTAGCAGTAGAACAAGAAACAAACAATGATGTTCGTC  >  NC_000913/2560223‑2560419

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: