Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,574,0510GT44.0% 161.2 / 38.9 85V384V (GTC→GTA) pqqLputative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base G (5/42);  new base T (37/0);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.83e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGC  >  NC_000913/1573907‑1574201
                                                                                                                                                |                                                                                                                                                          
tttCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTAcatc                                                                                                                                                       >  2:172757/1‑150 (MQ=255)
  tCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGAcatcat                                                                                                                                                     <  2:1257843/150‑1 (MQ=255)
  tCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGAcatcat                                                                                                                                                     <  1:858432/150‑1 (MQ=255)
     aaaGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTAcatca                                                                                                                                                      >  1:869194/1‑146 (MQ=255)
     aaaGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGAATTTAcaacattgtg                                                                                                                                                 >  2:1067481/1‑147 (MQ=255)
       aGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑Gtt                                                                                                                                             <  2:822131/151‑1 (MQ=255)
       aGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAAATTTTAAATTATCG‑A‑Gtt                                                                                                                                             >  1:211120/1‑151 (MQ=255)
              atGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTcagc                                                                                                                                      <  1:206674/151‑1 (MQ=255)
               tGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTcagc                                                                                                                                      <  1:418307/150‑1 (MQ=255)
                gTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTcagcag                                                                                                                                    <  2:28217/151‑1 (MQ=255)
                gTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTcagca                                                                                                                                     <  1:554136/150‑1 (MQ=255)
                   tctTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAAATTTTACATCATCG‑A‑GGTCTTCAGCAGaa                                                                                                                                  >  2:989134/1‑150 (MQ=255)
                     ttCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGTGAGATTTTACATCATCG‑C‑GAGCTTTCGCAGaaa                                                                                                                                 >  2:764730/1‑149 (MQ=255)
                        cGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCa                                                                                                                            <  1:501552/151‑1 (MQ=255)
                         ggCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCAt                                                                                                                           <  2:13844/151‑1 (MQ=255)
                           cGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCG‑A‑GTTCTTTTGCAGAAAAAACa                                                                                                                            >  2:1312004/1‑148 (MQ=255)
                            gACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTTTTTAGCAGaaaaaac                                                                                                                             >  1:58639/1‑144 (MQ=255)
                                  aaCGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTAATCAACGAA‑TTTTTTCAACA‑AAAAAACATTCTGAt                                                                                                                     >  2:1382309/1‑148 (MQ=255)
                                   aCGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTATATTTTACATCATTC‑A‑GTTTTTCAGCAGGAAAAACaagttggt                                                                                                                     >  1:99857/1‑140 (MQ=255)
                                   aCGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAt                                                                                                                <  2:285947/151‑1 (MQ=255)
                                   aCGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAt                                                                                                                <  1:764602/151‑1 (MQ=255)
                                    cGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑g                                                                                                              <  1:752699/151‑1 (MQ=255)
                                    cGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTGCATTATCC‑A‑GTTTTTCAGAAGAAAAAACAAGCTGGTC‑AAt                                                                                                                >  2:152459/1‑150 (MQ=255)
                                        gTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTg                                                                                                          <  1:691330/151‑1 (MQ=255)
                                        gTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTg                                                                                                          <  2:506914/151‑1 (MQ=255)
                                         tattatTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTAAATCATCG‑A‑GTTTTTTTGCAGAAAAAACATGCTGGTC‑AAT‑gggt                                                                                                           >  2:230793/1‑147 (MQ=255)
                                               ttAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAAct                                                                                                    <  1:762668/150‑1 (MQ=255)
                                                tAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTCATCATCC‑A‑GGTTTTTTGCAGAAAAAACAAACTGATC‑AATGGGTTGCCa                                                                                                       >  2:81696/1‑147 (MQ=255)
                                                 aaTGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCG‑A‑GTTTTTCAGCCGAAAAAACATGCTGATC‑AAT‑GGGTGccca                                                                                                      >  2:688207/1‑144 (MQ=255)
                                                 aaTGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAGTTTTCATCATCG‑A‑GTTCTTCAGCAGAAAAAACATGCTGATC‑AAT‑TGGTGCCCAatct                                                                                                  >  2:380477/1‑150 (MQ=255)
                                                       cTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCAACA‑C‑GTTTTTCATCTGAAAAAACATGCTGATC‑AAT‑CGTTGGCCACACTGCCAtt                                                                                            >  1:1417421/1‑150 (MQ=255)
                                                               cAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCG‑A‑GGTTTTTAGCAGAAAAACCAAGCTGGTT‑AAA‑TGGTGCCAACTCTGCCCTTAATGCAtt                                                                                    >  1:531213/1‑150 (MQ=255)
                                                               cAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCG‑A‑GGTTTTCAGCAGAAAAAACATGCTGCTC‑AAT‑GGTTTCCCACTCTGCCCTTAAAGCAtt                                                                                    >  1:1220546/1‑150 (MQ=255)
                                                               cAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTc                                                                                   <  1:1003171/151‑1 (MQ=255)
                                                                aGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCg                                                                                  <  1:984023/151‑1 (MQ=255)
                                                                   cGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTCATCATCG‑A‑GTTTTTTAGCAGAAAAAACAAGCTGATC‑AAA‑GGTTTCCCAATCTTCCCTTTATGCATTCgc                                                                                 >  1:902740/1‑149 (MQ=255)
                                                                   cGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCgcag                                                                               <  2:303670/151‑1 (MQ=255)
                                                                   cGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCgcag                                                                               <  1:423445/151‑1 (MQ=255)
                                                                     aCTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCG‑A‑GGTCTTTAGCAGAAAAAACATGCTTGTC‑AAA‑GGTTGCCCACCCTGCCCTTAAAGCATTCccagc                                                                              >  1:1334524/1‑150 (MQ=255)
                                                                      cTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCG‑A‑GGTCTTTAGCAGAAAAAACATGCTGGTT‑AAT‑GGGTGCCCAATCTTCCATTAATGCAtt                                                                                    >  1:1169123/1‑143 (MQ=255)
                                                                       tGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGCTAAATTTCATATAATCA‑A‑ATTCTTTTCTAGAAAAAACAAGCTGATC‑AAA‑TGTTGCCAACCCTCCCATTAATGcgtattcccc                                                                              >  1:203813/1‑139 (MQ=255)
                                                                         gTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGCAGGTTTTGAATCATCC‑C‑GCTTTTCTGCCGAAAAACCACGCTGATC‑AAA‑TGGTGCTAAATCTGCCATTTATTCATTTGTAGCATcc                                                                          >  2:306924/1‑150 (MQ=255)
                                                                           cAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTg                                                                        <  1:15899/150‑1 (MQ=255)
                                                                           cAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTAATCATCG‑A‑GGTTTTTCTCAGAAAAAACAATCTGGTT‑AAT‑TGGTGGCAAATCCTCCATTTATGCATTCGCAGCATcc                                                                          >  2:472004/1‑148 (MQ=255)
                                                                             gCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTGAATCATCG‑G‑GGTTTTTTTCAGAAAAAACATGCTGATC‑AAA‑TGTTGCCAAATCTTCCCTTAATTCAtt                                                                                    >  2:1002434/1‑136 (MQ=255)
                                                                                tACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTGCATCATCG‑A‑GTTCTTTAACAGAAAAAACATGCTGATC‑AAT‑GGTTGCCAACTCTTCCATTAATGCATTCGCCGCATCCTTCAtatt                                                                  >  1:627947/1‑151 (MQ=255)
                                                                                  cGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCAtatt                                                                  <  2:345134/149‑1 (MQ=255)
                                                                                      aaTCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCG‑A‑GTTCTTTTTCAGAAAAAACCAGCTGGTC‑AAT‑GGTTGCCCACTCTGCCCTTAATGCATTTCCAGCATCCTGCCTATTATTGt                                                             >  2:1091214/1‑150 (MQ=255)
                                                                                         cACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGtctc                                                          <  2:1333843/150‑1 (MQ=255)
                                                                                           cGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCATCG‑A‑GTTCTTTTGCAGAAAAAACAATCCGATC‑AAA‑GGTTGCCAACTCTTCCCTTAATTCATTTGCCGCATCCTCCATATTATCTtctc                                                          >  2:317821/1‑148 (MQ=255)
                                                                                               cAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCAt                                                    <  2:26190/150‑1 (MQ=255)
                                                                                                aGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCAttt                                                  <  1:570105/151‑1 (MQ=255)
                                                                                                aGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCATCC‑AGGTTTTTTCGCCGAAAAA‑CACGCTGTTT‑AAT‑TGTTTCCCACCCCTCCACTTAAGCATTCTCCGCATTCTCCCTATTTTCTTCTCCGGCAttt                                                  >  2:1383402/1‑151 (MQ=255)
                                                                                                 gCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑G‑GGTTTTTTGCCGAAAAAACATGCTGATC‑AAT‑GGGTGCCCACTCTTCCATTTATGCATTTTCAGCATCCTGCATATTATCTTCTCGTGCAttt                                                  >  2:22154/1‑150 (MQ=255)
                                                                                                  cTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTGCATCAACG‑A‑GGTTTTTTACAGAAAAAACAAGCTGGTCAAAT‑GGTT‑CCCACTCCTCCCTTTATTAATTCGCAGCAACCTGCATTATATTGTCTCggggaatt                                                  >  1:590134/1‑142 (MQ=255)
                                                                                                    tGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTaca                                               <  1:571960/150‑1 (MQ=255)
                                                                                                     gCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACg                                             <  1:1037106/151‑1 (MQ=255)
                                                                                                     gCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACg                                             <  1:501677/151‑1 (MQ=255)
                                                                                                     gCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTGCATCATCG‑A‑GTTTTTCAGCAGAAAAAACAAGCTGATC‑AATGGGTTGCC‑ACTCTTCCCTTAATGCATTTGCAGCATCCTGCATATTTTCGTCTCGTGCAttt                                                  >  2:596328/1‑146 (MQ=255)
                                                                                                      cTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAAATTTTACATTAACC‑A‑GTTTTTCTCTAGAAAAACCATTCTGATC‑AAT‑GGGTGTCAACCCTGCCCTTTATGCATTTGCAGGATCCTTAATATTATTGTCTCCTGCATTTTCACg                                             >  1:878640/1‑150 (MQ=255)
                                                                                                         aTCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTCATCATCG‑A‑GGTTTTTAGCAGAAAAAACAAGCTGATC‑AAA‑GGTTGGCAACTCTTCCCTTTATACAATTGCAGCAACCCTCCTATTATTGTCTCCTGCATTTTCAAggga                                          >  1:701649/1‑148 (MQ=255)
                                                                                                          tCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCG‑A‑GTTCTTTAGCAGAAAAAACATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGaaa                                         >  1:314623/1‑150 (MQ=255)
                                                                                                           cAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGaaaa                                        <  1:414318/150‑1 (MQ=255)
                                                                                                           cAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTTAATCATCG‑G‑GTTCTTCAGCAGAAAAAACATGCTGATC‑AAA‑GGGTGGCCAATCTGCCATTTATGCATTTGCAGCATCCTGCATAATTTCTTTTCGTGCATTTTAACGGaaaa                                        >  1:1162290/1‑150 (MQ=255)
                                                                                                             aCCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCG‑A‑GGTTTTTAGCAGAAAAAACATGCTGATC‑AAA‑GGGTGCCAACTCCGCCCTTTATTCATTCGCCGCATCCCGCCTATTATTGTCTCCTGCATTTTCACCGaaaaa                                       >  2:347974/1‑149 (MQ=255)
                                                                                                                gCATTTTTCAGCCAGGTGAGGCGGTTAAATTTTACAAAATCG‑A‑GGTCTTCAACAGAAAAAACCTGCTGATC‑AAT‑AGTTGCCAACTCTTCCCTTAATGCATTCGCAGCATTCTGCCTATTTTCGTCTCCTGCATTTTCACGGAAAAACa                                     >  1:863129/1‑148 (MQ=255)
                                                                                                                 cATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGa                                 <  1:726733/151‑1 (MQ=255)
                                                                                                                 cATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCG‑A‑GTTTTTTTGCCGAAAAAACATGCTGATC‑AAT‑GGGTGCCCAATCTGCCATTAATGCATTCGCCGCATCCTGCCTATTATTGTTTCGTGCCTTTTAACGGgaaaa                                       >  2:767932/1‑145 (MQ=255)
                                                                                                                    ttttCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAA‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACt                               <  1:706886/150‑1 (MQ=255)
                                                                                                                     tttCAGCCAGGTGAGGCGGGGAGATTTTACATCATCC‑A‑GGTTTTTAGCCGAAAAAACCAGCTGATC‑AAT‑GGTTGCCAACTCTTCCCTTAATGCATTTGCAGCATCCCTCATATTATTCTCTCCTGCATTTTTACGGgaaaaaa                                     >  1:426776/1‑141 (MQ=255)
                                                                                                                      ttCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGa                             <  1:508366/150‑1 (MQ=255)
                                                                                                                        caggcaggTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATaa                          <  1:448036/151‑1 (MQ=255)
                                                                                                                         aGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAAt                         <  1:1297933/151‑1 (MQ=255)
                                                                                                                         aGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAAt                         <  1:172764/151‑1 (MQ=255)
                                                                                                                            cAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCgg                      <  1:8438/151‑1 (MQ=255)
                                                                                                                              ggTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGt                    <  1:750836/151‑1 (MQ=255)
                                                                                                                                  aGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCaa                <  1:891823/151‑1 (MQ=255)
                                                                                                                                   ggCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAt               <  1:252897/151‑1 (MQ=255)
                                                                                                                                     cGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAttt             <  2:58639/151‑1 (MQ=255)
                                                                                                                                      gggTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAAtttt            <  2:898788/151‑1 (MQ=255)
                                                                                                                                         tAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAAc         >  1:624958/1‑151 (MQ=255)
                                                                                                                                          aGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACg        <  2:371280/151‑1 (MQ=255)
                                                                                                                                          aGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACg        >  1:687430/1‑151 (MQ=255)
                                                                                                                                            aTTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCt      <  2:1304604/151‑1 (MQ=255)
                                                                                                                                               tGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcg   >  1:509516/1‑151 (MQ=255)
                                                                                                                                                gACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcgc  <  2:90996/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                          
TTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCG‑A‑GTTCTTCAGCAGAAAAACCATGCTGATC‑AAT‑GGTTGCCAACTCTGCCATTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGC  >  NC_000913/1573907‑1574201

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