Read alignment evidence...
  seq id position ref A a var freq score (cons/poly) reads annotation genes product
*NC_0009132,074,0440TTC36.9% 239.5 / 18.7 123S836P (TCC→CCC) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  major Allele T (37/40);  minor allele C (0/45);  total (37/86)
Fisher's exact test for biased strand distribution p-value = 5.42e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

CGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCG  >  NC_000913/2073803‑2074294
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
                                        atcaggggctgtattGGTTTCTTCCGTGTGGGGGGTGACCTGAGGAGACCAGGGGTGCGCGGTTTTTCTTGGCCCGGGGGGGTATATGGTGCGGCGGGCCTTTCTTCCGTGGATGTAAAGGATGAGGCCGCCCCCCGTGCCGGCCCGGGCGGGGATGGGGCGGGCGGCGGGGGCGGATCCCAAATGCGGGACCACCCGCCCCCCGGCATGTGGGCGAACTTTGGGGCGGGGGGAACCCGCCACAGCATGaa                                                                                                                                                                                                           <  2:2105648/237‑1 (MQ=255)
                                                  gCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGtcctc                                                                                                                                                                                                                                                             <  2:3757729/193‑1 (MQ=255)
                                                              gccgtgTGGGGGGTGACGTGAGGAGAGCAGAGGTTGCGGGTATGTCGGTGCCCGGGGGGGTATGTGGGGCTGCTGGCCATTCTCCCGTGGATGTAAAGGAGGAGGCCGCCCCCCGGGCCGGCGGGGGCGGGGGTGAGGCCGGCGGCCGGGGGGAATCCTTGAATCGGGAACACCCGCCCCCCGGGGGGGGGGCTGAAATTGGGGCAAAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGa                                                                                                                                                                                      <  2:2314749/249‑1 (MQ=255)
                                                                  tcgggggggTGACTTGAGGAGAACGGGGGTTGCCGTTTGGTCGGGGCCCGGGGGGGTATTGGGTGCTGCGGCCCTTTCTTCGTTGGATTTAAGGGAGGATGAGGGCCCCGGGGCCGGCAGGGTCGGGGAGGGGGCGGGCAGCCGGGGGGGATCCCTGAATCGGGTACACACGCCCCCCGGCTTGGGGGCTGACATTGGGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTc                                                                                                                                                                                  <  2:3514839/244‑1 (MQ=255)
                                                                  tcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAAcc                                                                                                                                                                                                                         <  2:1211599/211‑1 (MQ=255)
                                                                            gACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGc                                                                                                                                                                                                                                   <  1:665688/191‑1 (MQ=255)
                                                                                 gatgaGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACa                                                                                                                                                                                                                                          <  2:2347151/179‑1 (MQ=255)
                                                                                      gaACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACAt                                                                                                                                                                                                                                         <  2:3213947/175‑1 (MQ=255)
                                                                                      gaACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTACCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGcc                                                                                                                                                                                                                                                       <  2:2640891/162‑1 (MQ=255)
                                                                                             ggTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTgg                                                                                                                                                                                                                                                  <  2:5839547/160‑1 (MQ=255)
                                                                                              gtggcgggTATGTTGTGGCCGGGGGGTGTATAGGGGGCTGCGGCCCTTTCTTCCGTGGTGGGAAGGGAGGATGACGGCGCCCGGGCCGGCACGGTCGGGGAGGAGGCCGGAGGCGTGGGGGGATCCCTGAATCGGGTCCACACGCCCTCGGGCTTGGGGGCTGACTTTGGGCCACGGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCAGGGCTGGGGCTGGCTGGGCt                                                                                                                                                      <  1:5825424/244‑1 (MQ=255)
                                                                                                         ttCTGTGCCCGGGGGGTTATGTGGTGCTGCTGGCCTTTCTCCCGTTGTTGTAAAGGATGATGACGGCCCCCGTGCCGGCCCGGTCGGGGATGACGCCGGCAGCGGGGGCGGATCCCTGAATCGGGTACACCCGCCCCCCGGCCTGGGGGCTGACTTTGTGGCACAGGGAACCCGCCACAGCATGAAAg                                                                                                                                                                                                         <  2:1427019/187‑1 (MQ=255)
                                                                                                           cTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACAt                                                                                                                                                                                                                                         <  1:2284076/154‑1 (MQ=255)
                                                                                                                      ggggTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCTGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCCGGTACACACGTCCTCCGGCCTGTGGGCTGACACTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGCCATCGGACAACAACGACTTCCGCGCCCGGGGCTGGGGCTGGCCGGGGTCACCCGCAACCGGTCTGcccccttg                                                                                                                             >  1:1859298/1‑246 (MQ=255)
                                                                                                                      ggggTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGcc                                                                                                                                                                                                                                                       <  2:4444724/130‑1 (MQ=255)
                                                                                                                      ggggTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCtgt                                                                                                                                                                                                                                                    <  2:1529654/133‑1 (MQ=255)
                                                                                                                      ggggTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGcac                                                                                                                                                                                                                                 <  1:4318651/151‑1 (MQ=255)
                                                                                                                       ggggtagtgggggctgggggcttttcttcTGTTGATGTAAGGGATGATGACGGCCCCGGGGCGGGCGGGTCCGGGGATGAGGCGGGCGGCTGGGGCGGATACCTGAATCGGGTACACACGTCCCCCGGCCTGGGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAAcc                                                                                                                                           <  2:3309442/214‑1 (MQ=255)
                                                                                                                        ggttatggggTGCGCCTGGCCATTCTTCGTTTGATGTAAGGGATGAGGCGGGCCCCCGGGCGGGCACGGTCGGGGAGGAGGCGGGCAGCCGGGGCGGATCCCTGAATCGGGTACACCCGTCCCCCGGCGTGGGGGCTGACATTGGGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAg                                                                                                                             <  1:1234415/242‑1 (MQ=255)
                                                                                                                          ttattTGGTGCTGCGGGCCTTTTTTCCGTTGATGTAAAGGATGATGAGGGCTCCCGTGCGGGCCCGGTCGGGGAGGATGCCGGCAGCCGGGGCGGATACCTGAATCGGGTACACCCGTCCCCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGaa                                                                                                                                                                                                           <  1:3732293/169‑1 (MQ=255)
                                                                                                                          ttatggggTGCGGCTGGCTATTCTTCCTTTGTTGTTAAGGATGATGACGCCTCCGGTGCCGGCCCGGCCCGGTATGATGCGGGGGGCGGGGGGGGATCCCGAAATCTGGAACCCCCGCCCCCGGGCGTGGGGGCTGACATGGGGACAGGGGGAGCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTAt                                                                                                                          <  1:4957282/245‑1 (MQ=255)
                                                                                                                              gtGGTGCTGCTGCCCTTTCTCCGGTTGATGTTAGGGATGAGGCGGGCCCCGGTGCGGGCCGGTGCGGGGATGATCCGGGCAGCGGGGGGGGACCCCTGAATCTGGTACACACGCCCTCCGGCCTGTGGGCTGACATTGTGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAAcc                                                                                                                                           <  1:3348185/229‑1 (MQ=255)
                                                                                                                                 gtgtcgcgggCCATTCTTCCGTTTGTTTTAAGGGTGGGGGCGGCCCCCCGGCCGGCACGGGCCGGGGATGGGCCGGGCGGCCGGGGGGGACCCCTGAATCGGGTCCACGCGTCCCCGGGCTGGGGGGCTAACTTTGGGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTg                                                                                                                     <  2:4054300/242‑1 (MQ=255)
                                                                                                                                 gTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTgg                                                                                                                                                                                                                                                  <  2:3514831/125‑1 (MQ=255)
                                                                                                                                              ttctccCGTTGATGTTAAGGATGATGACGGCCCCCGTGCCGGCACGGTCCGGGATGATCCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGc                                                                                                                                                                                                                      <  1:1361945/139‑1 (MQ=255)
                                                                                                                                                   ccGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAAc                                                                                                                                                                                                                          <  1:881445/130‑1 (MQ=255)
                                                                                                                                                         aTGTTAGGGATGAGGCCGGCCCCCGTGCCGGCCCGGTCCGGGAGGAGGCCGGCAGCTGGGGCGGATCCCTGAATCGGGTACACCCCTCCCCCGGCCTGGGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAg                                                                                                                                                                                                         <  2:2917387/141‑1 (MQ=255)
                                                                                                                                                                gtgtgatggCGGCCCCCGGACCGGCGCGGGCCGGGATGGGGGCGGGGGCGGGGGGGGATACAGGATGCGGGTCCACAGGCCCGCGGGCCTGGGGGCGGGCGTTGGGGCAGGGGGACCCCGCCACAGCATGAAAGCGTCGTCGGACGATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGtat                                                                                     <  1:1967197/247‑1 (MQ=255)
                                                                                                                                                                 gatggtGCCGGCCCCGGTGCCGGCGCGGTCCGGGATGATGCCGGCGGCGTGGGGGGATCCCTGATTCTGGTCCACACGCCCCCGGGCCTGGGGGCTGACATGGGGGCAGGGGAAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCt                                                                                                        <  1:48167/230‑1 (MQ=255)
                                                                                                                                                                 gaggatgaCGGCCCCCGTCCGGGCGGGGCCCGGGAGGATCCCGGCAGCCGGGGCGGATCCCTGAATCTGGTACACACGCCCTCCGGCCTGGGGGCTGACATGGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGAc                                                                                                                                                                                              <  1:3839216/141‑1 (MQ=255)
                                                                                                                                                                        ccccgcaccGTGCCGCCCGGGCCGGGGAGGATGCCGGCGGCCGGGGCGGATCCCTGAATCTGGTCCACCCTCCCTCCGGCCTGTGGGTTGACTTTGTGCCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTAt                                                                                                                          <  1:5476342/198‑1 (MQ=255)
                                                                                                                                                                         cGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGcc                                                                                                                                                                                                                     <  2:608475/113‑1 (MQ=255)
                                                                                                                                                                                 gTGCCGGCAGGGTCGGGGAGGAGGCCGGCAGCCGGGGGGGATCCCTGAATCTGGAACACACGTCCTCCGGCCTGGGGGCTGACATTGGGGCACGGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCggggctggggc                                                                                                                                                                <  1:4691749/158‑1 (MQ=255)
                                                                                                                                                                                  tgcggccccGGTCCGGGATGGCGCCGGCGCCCTGGGGGGATCCCTAACTCTGTCCCCCCCTCCCCCGGGCGTGTGGGCTGACTTTGTGGCCCGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGc                                                                                                         <  1:161655/204‑1 (MQ=255)
                                                                                                                                                                                          cgggcacGGGATGGGGCGGGCCGCGTGGGGGGATCCCTAAATTGGGCCCCCCCGTCCCCCGGGTGGTGGGTGGAATTTGGGGCACGGGCAACCCGCCACAGCATGAAAGCGTCATCGGACAATAGCGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATAc                                                                                   <  1:989860/221‑1 (MQ=255)
                                                                                                                                                                                            gTCGGGGATGATGCGGGCACCTGGGGCGTTTCCTTGATTCTGGCACCCCCGCCCTCCGGCCTGTGGGTTGACATTGTGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATg                                                                                                          <  1:5936287/201‑1 (MQ=255)
                                                                                                                                                                                            gTCGGGGATGAGGCCGGCGGCCGGGGGGAATCCCTGATTTTGGTCCACCCCCCCCCCGGCTTGTGGGCTAAATTTGTGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGccc                                                                                                                                  <  1:3272829/177‑1 (MQ=255)
                                                                                                                                                                                            gTCCGGGATGATGCCGGCAGCGGGGGGGGATCCCTGAATCTGGTACACCCGCCCCCCGGCCTGGGGGCTAACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGggctggct                                                                                                                                                           <  1:5497714/152‑1 (MQ=255)
                                                                                                                                                                                            gTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTgg                                                                                                                                                                                                                                                  <  1:371448/66‑1 (MQ=255)
                                                                                                                                                                                            gTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTgg                                                                                                                                                                                                                                                  <  1:5208722/66‑1 (MQ=255)
                                                                                                                                                                                            gTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACAt                                                                                                                                                                                                                                         <  2:3449961/74‑1 (MQ=255)
                                                                                                                                                                                            gTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGc                                                                                                                                                                                                                                   <  2:5779751/80‑1 (MQ=255)
                                                                                                                                                                                            gTCCGGGAGGATGCCGGCGGCCGGGGGGGTTCCCTGAATTTGGTACACCCGTCCTCCGCCCTGTGGGCTGACATTGTGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAg                                                                                                                             <  1:1359579/182‑1 (MQ=255)
                                                                                                                                                                                             gcccGGGATGATGCCGGCGGCGTGGGCGGATCCCTGAATCTGGTACACCCTCCCTCCGGCCTGGGGGCTGACTTTGTGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAAcc                                                                                                                                           <  2:2498860/166‑1 (MQ=255)
                                                                                                                                                                                             gcccGGGAGGATGCCGGCAGCTGGGGCGGATCCCTAAATCGGGCACCCCCGCCCCCCGGCCTGGGGGCTGACATTGTGGCACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTAt                                                                                                                          <  2:1060159/183‑1 (MQ=255)
                                                                                                                                                                                             gcccGGGAGGATGCCGGCAGCCTGGGCGGATCCCTGAATCTGGTCCCCCCGTCCCCCGGCCTGTGGGCTGACATTGTGGCACGGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCggggc                                                                                                                                                                      <  1:1158340/139‑1 (MQ=255)
                                                                                                                                                                                              gccggggatgggcccGGCGGCCGGGGGGGAGCCCTGAATCGGGTACACCCCCCCCCCGGCCTGGGGGCGGACTTTGGGGCGCGGGGCACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGc                                                                                                         <  2:4904347/189‑1 (MQ=255)
                                                                                                                                                                                              gccggggATGATGCGGGCAGCCGGGGGGGATCCCTAAATCTGGCACACCCGCCCCCCGGCGTGGGGGCTGACTTTGGGGCACGGGAACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGggctggct                                                                                                                                                           <  2:483202/148‑1 (MQ=255)
                                                                                                                                                                                              gccggggAGGATCCGGGCGGCTGGGGCGGATCCCTGAATCGGGTACACCCCCCCTCCGGCCTGGGGGCTGACATTGTGGCACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAAcc                                                                                                                                           <  2:4328287/164‑1 (MQ=255)
                                                                                                                                                                                              gccgggaaggggccgggCAGCCGGGGCGGATCCCTAAATCTGGTACCCCCGTCCCCCGCCCGGGGGGCTGACATGGGGGCACAGGAAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAAcc                                                                                                                                           <  2:5003428/153‑1 (MQ=255)
                                                                                                                                                                                               gttgtgtgtaggggcgggcccggggCGGACCCCTAAATCGGGTCCACCCGCCCCCCGGCCTGGGGGCGAATATTGGGCCACGGGAACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACgg                                                        <  1:3401378/228‑1 (MQ=255)
                                                                                                                                                                                                   ggATGATGCCGGCAGCCGGGGCGGATCCCTAAATCTGGTACCCCCGCCCCCCGGCCTGTGGGCTGACTTTGTGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGc                                                                                                         <  2:5045278/197‑1 (MQ=255)
                                                                                                                                                                                                    gatggtGCCGGCAGCCGGGGGGGATCCCTAATTCTGGTCCACACGTCCCCCGGCTTGTGGGCTGACATTGTGGCACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTg                                                                                                                                                 <  1:1082371/156‑1 (MQ=255)
                                                                                                                                                                                                      tggaggcGGGAGGCTGGGGGGATACCTGAAACCGGTAACCCCGCCCCCCGGGTGGGGGGTCAGAAATGGGGGCACGGGAACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGAc                                                 <  1:5241073/245‑1 (MQ=255)
                                                                                                                                                                                                          gCCGGCAGCGTGGGCGGATCCTGGAATCTGTCACACCCGTCCCCCGGCGGGTGGGCTGACATTGGGGCAAGGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGAc                                                 <  2:5082792/246‑1 (MQ=255)
                                                                                                                                                                                                                 ggcgggggggATCCCTGAATGTGGTACACACGTCCTCGGGCGTGGGGGTTGACTTTGGGGCAAAGGAAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCt                                                                                                                                                      <  1:3732420/134‑1 (MQ=255)
                                                                                                                                                                                                                 gcgggggCGGATCCCGAAATCGGGTACACCCGCCCCCCGGCCTGGGGGCTGAAATTGGGGCACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCggggctgggg                                                                                                                                                                 <  1:4235925/123‑1 (MQ=255)
                                                                                                                                                                                                                 gCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGcac                                                                                                                                                                                                                                 <  2:3157986/63‑1 (MQ=255)
                                                                                                                                                                                                                 gCCGGGGCGGATCCCTGAATCTGGTACACCCGCCCCCCGGCCTGGGGGCTGACATTGGGGCACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTAt                                                                                                                          <  2:2221421/166‑1 (MQ=255)
                                                                                                                                                                                                                    tgcgcttatcccTCAATCGGGCACACACGTCCCCCGGCGTGTGGGCTAACTTTGGGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTAtgtg                                  <  2:4380969/241‑1 (MQ=255)
                                                                                                                                                                                                                            tACCTGAATCTGGTACACACGTCCTCCGGCCtgtg                                                                                                                                                                                                                                                   >  1:2146673/1‑35 (MQ=255)
                                                                                                                                                                                                                             cccTGAATCTGGTACCCCCGTCCCCGGGCTTGTGGGCTGACTTTGGGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAAc                                                                                             <  1:3902132/182‑1 (MQ=255)
                                                                                                                                                                                                                             aCCTGAATCTGGTACACACGTCCTCCGGCCTGTgg                                                                                                                                                                                                                                                  >  2:3901479/1‑35 (MQ=255)
                                                                                                                                                                                                                             aCCTGAATCTGGTACACACGTCCTCCGGCCTGTgg                                                                                                                                                                                                                                                  >  1:1058994/1‑35 (MQ=255)
                                                                                                                                                                                                                              ctggAATCGGGAACACACGTCCTCCGGCTGGTGGGCTGACATTGGGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCggggctgggg                                                                                                                                                                 <  1:1412762/111‑1 (MQ=255)
                                                                                                                                                                                                                              ccTGATTCTGTTCCACCCCCCCCCCGCCTTGGGGGCTGACATGGGGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTg                                                                                                                     <  1:1233735/158‑1 (MQ=255)
                                                                                                                                                                                                                              ccTGAATCTGGTACACACGTCCTCCGGCCTGTggg                                                                                                                                                                                                                                                 >  1:376780/1‑35 (MQ=255)
                                                                                                                                                                                                                              ccTGAATCTGGTACACACGTCCTCCGGCCTGTggg                                                                                                                                                                                                                                                 >  2:486977/1‑35 (MQ=255)
                                                                                                                                                                                                                              ccTGAATCTGGTACACACGTCCTCCGGCCTGTggg                                                                                                                                                                                                                                                 >  2:1233483/1‑35 (MQ=255)
                                                                                                                                                                                                                              ccTGAATCGGATACACCCGTCCCCGGGCCTGGGGGCTGACTTTGTGGCCCGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCAc                                                                                                                       <  1:487081/156‑1 (MQ=255)
                                                                                                                                                                                                                                   ctCTGGCACACCCCCCCCCCGCCTTGGGGGCTGCCTTTGGGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAAcc                                                                                                               <  2:5362896/158‑1 (MQ=255)
                                                                                                                                                                                                                                   aTCTGGTACACACGTCCTCCGGTCGGTGGGGGGACAGTGGGGACGGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGc                   <  1:2034731/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTgg                                                                                                                                                                                                                                    <  2:2627246/40‑1 (MQ=255)
                                                                                                                                                                                                                                     cTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAAcc                                                                                                                                           >  1:5794496/1‑129 (MQ=255)
                                                                                                                                                                                                                                     cTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGc                   >  2:3276894/1‑249 (MQ=255)
                                                                                                                                                                                                                                     cTGGTACACACGTCCTCCGGCCTGTGGGCTGAAATTGTgg                                                                                                                                                                                                                                    >  1:2627684/1‑40 (MQ=255)
                                                                                                                                                                                                                                       ggTCCACACGCCCCCCGGCCTGGGGGCTGACTTTGGGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCggggctgggg                                                                                                                                                                 <  2:2861476/105‑1 (MQ=255)
                                                                                                                                                                                                                                       ggTACACCGGCCCCCCGGCTGGGGGGCGGACATGGGGGCAAGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGAt                                                                                                           <  1:1197587/159‑1 (MQ=255)
                                                                                                                                                                                                                                       ggTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCggggctgggg                                                                                                                                                                 >  1:2861932/1‑105 (MQ=255)
                                                                                                                                                                                                                                       ggTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCt                                                                                                                                                      >  2:2356934/1‑116 (MQ=255)
                                                                                                                                                                                                                                       ggTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTAt                                                                                                                          >  1:644515/1‑144 (MQ=255)
                                                                                                                                                                                                                                       ggTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTg                                                                                                                     >  1:3455668/1‑149 (MQ=255)
                                                                                                                                                                                                                                       ggTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGAt                                                                                                           >  2:1197341/1‑159 (MQ=255)
                                                                                                                                                                                                                                       ggTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCTGGCATGGCAGTg               >  2:288347/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTg               >  1:1703003/1‑251 (MQ=255)
                                                                                                                                                                                                                                          acacacGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAg                                                                                                                             >  1:1508975/1‑138 (MQ=255)
                                                                                                                                                                                                                                             cacGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTc                                                                                                                                                     >  1:4104833/1‑111 (MQ=255)
                                                                                                                                                                                                                                             cacGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGc                                                                                                                                    >  2:218884/1‑128 (MQ=255)
                                                                                                                                                                                                                                             cacGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTAt                                                                                                                          >  2:4493832/1‑138 (MQ=255)
                                                                                                                                                                                                                                             cacGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTAt                                                                                                                          >  1:291923/1‑138 (MQ=255)
                                                                                                                                                                                                                                             cacGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATAc                                                                                   >  2:3142406/1‑177 (MQ=255)
                                                                                                                                                                                                                                             cacGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCacaa          >  1:655907/1‑250 (MQ=255)
                                                                                                                                                                                                                                               cGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACt                                                                                                                                                                                    >  1:25043/1‑78 (MQ=255)
                                                                                                                                                                                                                                                  cctccGGCTTGGGGGCTGACATTGTGGCACGGGAACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGAc                                                 <  1:3816220/206‑1 (MQ=255)
                                                                                                                                                                                                                                                  cctccGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCt                                                                                                                                                      >  2:5750447/1‑105 (MQ=255)
                                                                                                                                                                                                                                                  cctccGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAg                                                                                        >  1:254768/1‑167 (MQ=255)
                                                                                                                                                                                                                                                  cctccGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATAcc                                                                                  >  1:5210353/1‑173 (MQ=255)
                                                                                                                                                                                                                                                  cctccGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGAc                                                 >  2:3815615/1‑206 (MQ=255)
                                                                                                                                                                                                                                                  cctccGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGGAAGGACAACGCCGGTTATGTGAAGTTCggg                         >  1:2848465/1‑230 (MQ=255)
                                                                                                                                                                                                                                                    tccGCCCTGTGGGTTGAAATTGGGGCACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATgtgcg  <  1:4559683/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
CGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCG  >  NC_000913/2073803‑2074294

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: