Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,713,5390A.37.4% 312.9 / 185.7 147coding (448/564 nt)tag3‑methyl‑adenine DNA glycosylase I, constitutive
Reads supporting (aligned to +/- strand):  ref base A (43/49);  new base . (25/30);  total (68/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 156 total aligned reads displayed.

AAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGGCTGCTGTTGCTATCCGGGAAATAAACCATGATTCGGGAAGCGCAACGTTCA  >  NC_000913/3713402‑3713675
                                                                                                                                         |                                                                                                                                        
aaaTGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGaa                                                                                                                                         >  1:219241/1‑139 (MQ=255)
      aaCAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGaaaaacgt                                                                                                                                   >  2:453592/1‑136 (MQ=255)
           aaCGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGtt                                                                                                                              >  2:1952705/1‑139 (MQ=255)
           aaCGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGtt                                                                                                                              >  2:1687697/1‑139 (MQ=255)
           aaCGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGttt                                                                                                                             >  2:980018/1‑139 (MQ=255)
           aaCGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGttt                                                                                                                             >  2:1805128/1‑139 (MQ=255)
             cGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGtttt                                                                                                                            >  1:68564/1‑139 (MQ=255)
                cGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAg                                                                                                                         >  1:1781661/1‑139 (MQ=255)
                  aaCCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGttt                                                                                                                      <  2:1165048/139‑1 (MQ=255)
                        ttGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCgg                                                                                                                 <  2:955353/139‑1 (MQ=255)
                           tCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGcac                                                                                                              <  2:737781/139‑1 (MQ=255)
                              aCTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACcaca                                                                                                          <  2:1156545/139‑1 (MQ=255)
                               cTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACcaca                                                                                                          <  2:528303/139‑1 (MQ=255)
                                tttGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACaa                                                                                                         <  2:68564/139‑1 (MQ=255)
                                      tGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGtt                                                                                                   <  2:435038/139‑1 (MQ=255)
                                      tGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTa                                                                                                  <  1:1867290/139‑1 (MQ=255)
                                        gTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGaaaaacg                                                                                                                                    >  1:487334/1‑102 (MQ=255)
                                           gTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCt                                                                                             <  1:399365/139‑1 (MQ=255)
                                            tttGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGTGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCtt                                                                                            >  1:1457123/1‑139 (MQ=255)
                                            tttGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCtt                                                                                            >  1:680646/1‑139 (MQ=255)
                                             ttGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAg                                                                                                                         <  2:563740/110‑1 (MQ=255)
                                             ttGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAg                                                                                                                         >  1:563740/1‑110 (MQ=255)
                                              tGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCttt                                                                                           >  1:341446/1‑139 (MQ=255)
                                              tGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTa                                                                                                                           >  2:577121/1‑106 (MQ=255)
                                              tGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTa                                                                                                                           <  1:577121/106‑1 (MQ=255)
                                                tAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATg                                                                                        <  2:476306/139‑1 (MQ=255)
                                                  aaTCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCa                                                                                      <  2:680646/139‑1 (MQ=255)
                                                   atcatcAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCa                                                                                      <  2:341446/139‑1 (MQ=255)
                                                     catcaGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACcaca                                                                                                          <  2:570358/117‑1 (MQ=255)
                                                     catcaGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACcaca                                                                                                          >  1:570358/1‑117 (MQ=255)
                                                      atcaGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGc                                                                                   <  1:961401/139‑1 (MQ=255)
                                                       tcaGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCgg                                                                                                                 >  1:921995/1‑108 (MQ=255)
                                                       tcaGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCgg                                                                                                                 <  2:921995/108‑1 (MQ=255)
                                                       tcaGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGt                                                                                                    >  1:1952021/1‑120 (MQ=255)
                                                       tcaGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCAt                                                                                 >  1:421164/1‑139 (MQ=255)
                                                        caGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTa                                                                                                  <  2:1277627/122‑1 (MQ=255)
                                                        caGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTa                                                                                                  >  1:1277627/1‑122 (MQ=255)
                                                        caGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCAtg                                                                                <  2:1717384/139‑1 (MQ=255)
                                                          gCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCAtgtg                                                                              >  1:1170820/1‑139 (MQ=255)
                                                            cacaGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTggg                                                                            <  2:421164/139‑1 (MQ=255)
                                                             acaGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGc                                                                           <  1:453592/139‑1 (MQ=255)
                                                              caGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGc                                                                           <  1:753597/139‑1 (MQ=255)
                                                               aGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTa                                                                                                  >  1:433252/1‑115 (MQ=255)
                                                                    acacAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGaa                                                                    <  2:720900/139‑1 (MQ=255)
                                                                     cacaAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTa                                                                                                  <  1:1369036/109‑1 (MQ=255)
                                                                      acaAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGcac                                                                                                              <  1:1758637/95‑1 (MQ=255)
                                                                      acaAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACcaca                                                                                                          >  1:1820560/1‑99 (MQ=255)
                                                                      acaAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATg                                                                  <  2:1170820/139‑1 (MQ=255)
                                                                        aaGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACcaca                                                                                                          >  1:989516/1‑97 (MQ=255)
                                                                         aGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACcaca                                                                                                          <  1:1366527/96‑1 (MQ=255)
                                                                         aGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACcaca                                                                                                          >  2:1366527/1‑96 (MQ=255)
                                                                         aGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCt                                                                                             <  2:1438470/109‑1 (MQ=255)
                                                                         aGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATc                                                               <  2:944240/139‑1 (MQ=255)
                                                                           ccACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTa                                                                                                  >  2:1751578/1‑103 (MQ=255)
                                                                           ccACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTa                                                                                                  <  1:1751578/103‑1 (MQ=255)
                                                                           ccACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCAt                                                             >  1:888120/1‑139 (MQ=255)
                                                                              caACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTTAAAAAACGTGGTTTTAAGTTTGTCGGcacca                                                                                                            >  1:440346/1‑90 (MQ=255)
                                                                              caACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGcacca                                                                                                            <  2:440346/90‑1 (MQ=255)
                                                                               aaCGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTa                                                                                                  <  2:238068/99‑1 (MQ=255)
                                                                                    tGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGttt                                                                                                                      <  2:540642/74‑1 (MQ=255)
                                                                                    tGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGttt                                                                                                                      >  1:540642/1‑74 (MQ=255)
                                                                                     gAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGgc                                                    <  1:759151/139‑1 (MQ=255)
                                                                                     gAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGgc                                                    <  1:1345961/139‑1 (MQ=255)
                                                                                     gAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGgct                                                   >  1:1967146/1‑139 (MQ=255)
                                                                                      aGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATc                                                               >  2:1545139/1‑127 (MQ=255)
                                                                                      aGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATc                                                               <  1:1545139/127‑1 (MQ=255)
                                                                                       gCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGgctg                                                  <  2:1227654/139‑1 (MQ=255)
                                                                                       gCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGgctg                                                  <  2:1781661/139‑1 (MQ=255)
                                                                                         gAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCg                                                                                                                  <  1:636463/73‑1 (MQ=255)
                                                                                         gAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCg                                                                                                                  >  2:636463/1‑73 (MQ=255)
                                                                                            aTTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTa                                                                                                  <  2:1903707/85‑1 (MQ=255)
                                                                                            aTTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTa                                                                                                  >  1:1903707/1‑85 (MQ=255)
                                                                                                  acaTCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGGCTGCTGTTGCTATc                                       <  2:1352785/139‑1 (MQ=255)
                                                                                                  acaTCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGGCTGCTGTTGCTATcc                                      <  2:1967146/139‑1 (MQ=255)
                                                                                                     tCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGcacca                                                                                                            >  1:1618236/1‑67 (MQ=255)
                                                                                                     tCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGcacca                                                                                                            <  2:1618236/67‑1 (MQ=255)
                                                                                                        aCGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTa                                                                                                  <  1:383811/74‑1 (MQ=255)
                                                                                                        aCGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGgctg                                                  <  2:2043154/122‑1 (MQ=255)
                                                                                                         cGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGGCTGCTGTTGCTATCCGGGaaa                                >  2:2030876/1‑139 (MQ=255)
                                                                                                         cGTCCGCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCg                                                                                                                  <  1:452612/56‑1 (MQ=255)
                                                                                                              gCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCAtg                                                                                >  1:1009573/1‑85 (MQ=255)
                                                                                                              gCCTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCAtg                                                                                <  2:1009573/85‑1 (MQ=255)
                                                                                                               ccTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCg                                                                                                                  <  1:1709907/50‑1 (MQ=255)
                                                                                                               ccTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCg                                                                                                                  >  2:1170111/1‑50 (MQ=255)
                                                                                                               ccTCCGACGCCCTATCTAAGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCg                                                                                                                  <  1:1170111/50‑1 (MQ=255)
                                                                                                                                 aGGCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTggg                                                                            <  2:853362/70‑1 (MQ=255)
                                                                                                                                  ggCACTG‑AAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGGCTGCTGTTGCTATCCGGGAAATAAACCATGATTCGGGAAGCGCAACg      >  1:1766208/1‑139 (MQ=255)
                                                                                                                                     aCTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGaa                                                                    <  2:364731/75‑1 (MQ=255)
                                                                                                                                       tGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTg                                                                                                     <  2:675678/40‑1 (MQ=255)
                                                                                                                                       tGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTg                                                                                                     >  1:675678/1‑40 (MQ=255)
                                                                                                                                       tGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGc                                                                           >  1:1999464/1‑66 (MQ=255)
                                                                                                                                       tGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGaa                                                                    >  1:286633/1‑73 (MQ=255)
                                                                                                                                       tGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTgg                                                     >  1:166646/1‑88 (MQ=255)
                                                                                                                                       tGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTgg                                                     <  2:166646/88‑1 (MQ=255)
                                                                                                                                       tGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGgctg                                                  >  2:453170/1‑91 (MQ=255)
                                                                                                                                       tGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGGCTGCTGTTGCTATCCgg                                    <  2:413358/105‑1 (MQ=255)
                                                                                                                                       tGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGGCTGCTGTTGCTATCCGGGAAATAAAcc                          <  2:188263/115‑1 (MQ=255)
                                                                                                                                       tGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGGCTGCTGTTGCTATCCGGGAAATAAACCATGATTCGGGAAGCGCAACGTTCa  <  1:1849653/139‑1 (MQ=255)
                                                                                                                                       tGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGGCTGCTGTTGCTATCCGGGAAATAAACCATGATTCGGGAAGCGCAACGTTCa  >  1:1004310/1‑139 (MQ=255)
                                                                                                                                        gAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTggg                                                                            <  2:1841334/64‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
AAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTCCGCCTCCGACGCCCTATCTAAGGCACTGAAAAAACGTGGTTTTAAGTTTGTCGGCACCACAATCTGTTACTCCTTTATGCAGGCATGTGGGCTGGTGAATGATCATGTGGTTGGCTGCTGTTGCTATCCGGGAAATAAACCATGATTCGGGAAGCGCAACGTTCA  >  NC_000913/3713402‑3713675

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑