Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,074,0660TG37.5% 79.0 / 19.8 59V843G (GTG→GGG) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base T (21/14);  new base G (0/21);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 1.96e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCC‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGTTT  >  NC_000913/2073831‑2074303
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
gCCACGCCGGAAAGCAGCGGCATCTATGGATTCGTCCGTCGGGAGGGTGGCCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCC‑GGCCTGTGGGCTGA‑CATTGTgg                                                                                                                                                                                                                                              <  2:2561226/238‑1 (MQ=255)
gCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCC‑GGCCTGTGGGCTGA‑CATTGTgg                                                                                                                                                                                                                                              <  1:2417607/238‑1 (MQ=255)
   aCGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCC‑GGCCTGTGGGCTGA‑CATTGTGGcaca                                                                                                                                                                                                                                          <  1:3467269/239‑1 (MQ=255)
     gCCGGAAAGCAGCGGACGTTTTGGCTTCGTCCGTCTGGAGGGTGACCTGAGGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCC‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑gg                                                                                                                                                                                                                                       <  2:1346564/239‑1 (MQ=255)
                 cGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACCCACGTCCTCC‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCcaca                                                                                                                                                                                                                           >  1:1988406/1‑239 (MQ=255)
                                                                            ggTCTGTGGCCGGGGGGGTATATGGTGCTGCGGGCCTTTCTTCCGTTGATGTTAGGGATGATGCGGGCTCCCGGGCCGGCACGGTCGGGGAGGAGGCGGGCAGCCGGGGCGGATCCCTGAATCGGGTACACGCGCCCTCC‑GGCTTGTGGGCTGA‑CTTTGGGGCACG‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGcc                                                                                                                                                                                     <  2:318416/217‑1 (MQ=255)
                                                                                                    gTGCTGCTGACCTTTCGTCCGTTGATGAGAAGGATGACGCCGGCCCCAGCGCCGGGCGGGTCGGGGAGGTCCCCGGCAGCCGGGGGGTAACCCTAAATCTGGACCCCCCGCCCCCC‑GGCTGGGGGGCTGA‑CATTGGGGCACG‑GGAACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCtt                                                                                                                                         <  2:240419/238‑1 (MQ=255)
                                                                                                                              gTTAAGGATGAGGACGGCCCCCGGGCCGGCACGGCCGGGGATGATGCCGGCAGCCGGGGCGGATACCTGAATCTGGTACACACGCCCCCC‑GGCCGGGGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCa                                                                                                                                                                                                             <  2:3450741/144‑1 (MQ=255)
                                                                                                                                      ggggggcGGCTCGCGGGCCGGCACGGCCGGGGAGGGTGCCGGCGGCGGGGGCGAATCCCTGAATCGGGCCCACCCGTCCCCC‑GGCGGGGGGGCTGA‑CTTGGGGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACaa                                                                                                       <  2:406845/232‑1 (MQ=255)
                                                                                                                                           aCGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCC‑GGCCTGTGGGCTGA‑CATTGTGGcaca                                                                                                                                                                                                                                          <  2:133331/103‑1 (MQ=255)
                                                                                                                                              gcccccGCGCCGGCACGGTCCGGGATGATGCCGGCACCCTGGGCGGATACCTGAATCTGGTCCACACGTCCTCC‑GGCCTGTGGGCTGA‑CATTGCGGCACA‑GGGAACCCGCCACAGCATGAAAGCGACATCGGACACTAACGACTTCCGCGCCCGGGGCTGGGGCCGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGGATCACTTACCACCTGATGCTGGAGTCACCACt                                                                                                     >  2:3307795/4‑232 (MQ=255)
                                                                                                                                                           cccggccggggATGATGCGGCCACGCGGGGGGAGTACCGGAGTCTGGACCAGAC‑TGCTCCGGGCCGGGGGGCGGA‑CATTGGGGAACA‑GGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGAc                                                                                 <  2:1448124/231‑1 (MQ=255)
                                                                                                                                                                tcggggATGATGCGGGCAGCCTGGGCGGATCCCTGAATCTGGTACACACGTCCTCC‑GGCCTGGGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTg                                                                                                                                              <  1:1403610/170‑1 (MQ=255)
                                                                                                                                                                       tgaggCGGGCAGCCTGGGCGGATACCTGAATCGGGTACACCCGTCCTCC‑GGCCTGGGGGCTGA‑CTTTGTGGCCCA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGAt                                                                      <  1:982981/238‑1 (MQ=255)
                                                                                                                                                                        ggggccggccgCGGGGGCGGATACCTGAAGCTGGTACCCACGTCCTCC‑GGCTTGTGGGCTGA‑CATGGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATg                                                                     <  1:687882/235‑1 (MQ=255)
                                                                                                                                                                            ccGGCAGCCTGGGCGGATACCTGAATCTGGTACCCACGTCCTCC‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAATGGAAACCGGTCGGCCCTCCAGCATCACTGACCACCCGATGCTGGAGCCCCAACTGCAGTATACCTCGCAGGGACATTCCATGGCAGACgat                                                                >  2:407643/1‑237 (MQ=255)
                                                                                                                                                                             cggcggCCGGGGGGGATACCGGAATCGGGTACACACGTCCCCC‑GGCCTGTGGGCTGA‑CATTGGGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGt                                                                <  1:2474505/238‑1 (MQ=255)
                                                                                                                                                                             cGGAAGCCGGGGCGGATACCTGATTTGGGTCCACACGTCCTCC‑GGCTTGTGGGCTGA‑CTTTGGGGCACA‑GGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGt                                                                <  1:2406601/238‑1 (MQ=255)
                                                                                                                                                                                    cggggCGGATACCAGAATCTGGTACCCCCGTCCTCC‑GGCCGGTGGGCTGA‑CATTGGGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCCTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGAc                                                          <  1:2234035/235‑1 (MQ=255)
                                                                                                                                                                                        ggggATCCCTGACTTTGATCCACCCGTCCTCG‑GGCCTGGGGGCGAA‑CATTGGGGCACG‑GGAAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGc                                                     <  1:1266471/236‑1 (MQ=255)
                                                                                                                                                                                        ggggACACCTGAACCTGTAACCCACTCCCCCC‑GGCGTGTGGGCTGA‑CATTGGGCAACA‑GGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGc                                                     <  1:3132987/236‑1 (MQ=255)
                                                                                                                                                                                        gCGGATACCTAAATCTGGTACACCCGCCCCCG‑GGCCTGTGGGCTAA‑CTTTGGGGCACG‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGc                                                     <  1:1967024/238‑1 (MQ=255)
                                                                                                                                                                                              cccggAATCTGGTACACCCGCCCTCG‑GGCCTGGGGGCTGA‑CTTGGTGGCACG‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATg                                                                     <  1:2729814/212‑1 (MQ=255)
                                                                                                                                                                                                 ggattgtGGTCCACCCGCCCCCG‑GGCCTGTGGTTTGA‑CTTGGGGGCACA‑GGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTAtgtg                                           <  1:1054021/233‑1 (MQ=255)
                                                                                                                                                                                                     gtggacacacccGCCCCCC‑GGCGTGGGGTCAGA‑ATTTGGGCCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGt                                       <  2:174268/233‑1 (MQ=255)
                                                                                                                                                                                                      cTGGTAAACACGCCCTCC‑GGCCTGTGGGCTGA‑CATTCTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGCCCTCGGACAATAACGACTTCCGCGCCCAGTGCTGGGGCTGGGTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCACCAACTGCAGTATACCAGGCAGGGACTTCCCCTGGATGACGGTAAGGACAACGCTGATTATGTGAAGt                                       >  2:877994/1‑238 (MQ=255)
                                                                                                                                                                                                          tACACACGACCTCC‑GGCCTGTGGGCTGA‑CATTGCGGCACA‑GGGAACCCGCCACAGCATGAAATCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCGGGGGCTGGCTGGGCTCACCGGAAACCGGTCTTCCCTCCAGTATCACTGACAACATGAGGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTACGGACAACGCCGGTCATGGGAAGTTCggg                                  >  2:1430549/1‑239 (MQ=255)
                                                                                                                                                                                                            cccccgcccccc‑GGCTTGGGGGCTGA‑CTTGGGGGCACG‑GGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCa                                <  1:166289/229‑1 (MQ=255)
                                                                                                                                                                                                            ccaccGTCCGCC‑GG‑CTGGGGGGTGAGCTTTGGGGCGCG‑GGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCa                                <  2:3310691/235‑1 (MQ=255)
                                                                                                                                                                                                            cacccGCCCTCC‑GGCCGGTGGGCGGA‑AATTGGGCCCAA‑GGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCa                                <  1:535815/239‑1 (MQ=255)
                                                                                                                                                                                                               cggTCCCCC‑GGCCAGTGGGCTAA‑CATGAGGGCACG‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATgg                             <  1:3252153/237‑1 (MQ=255)
                                                                                                                                                                                                               ccgcccccc‑GGCCTGTGGGCTGA‑CATGGGGCCACA‑GGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATgg                             <  1:242107/232‑1 (MQ=255)
                                                                                                                                                                                                               acGTCCTTC‑GGACAGTGGGCCGA‑AATTGCGGAACA‑GGGAACCCGACACAGCACGAAAGCGTCATCGGACAATAACGACTTCCGCGACCGGGGCTGGGACTGGCTGGGCTCACTGGTAACCGGCCGCCCGTTCAGTATCACTGACAACATGATGCTGGAGCCACAAACGCCGTATACCTGGCCGGGACTTTCCCTGGATGCAGGTAAGGAACACGCCGGTTATGTGCTGTTCgggaaag                              >  2:342881/1‑234 (MQ=255)
                                                                                                                                                                                                                cctccccc‑GGCCGGGGGGCTGA‑CTTGGGGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGc                            <  1:754077/237‑1 (MQ=255)
                                                                                                                                                                                                                 gTCCTCC‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGa                                                                       >  1:2562448/1‑195 (MQ=255)
                                                                                                                                                                                                                 gTCCTCC‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCa                           >  1:1877286/1‑239 (MQ=255)
                                                                                                                                                                                                                 gTCCTCC‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCa                           >  2:2006741/1‑239 (MQ=255)
                                                                                                                                                                                                                 gTCCTCC‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGATGGGGCTGGCTGGGCTCACTGGaaa                                                                                                                                                      >  1:1055853/1‑116 (MQ=255)
                                                                                                                                                                                                                   cctcc‑GGCTGGGGGGCGGA‑CATTGTGGCACA‑GGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAg                          <  1:2434749/238‑1 (MQ=255)
                                                                                                                                                                                                                   cctcc‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGt                         >  1:3272097/1‑239 (MQ=255)
                                                                                                                                                                                                                    ctcc‑GGCTGGTGGGTCGA‑‑ATTGGTGCACAGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTg                        <  1:1656526/239‑1 (MQ=255)
                                                                                                                                                                                                                    ctcc‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGt                         >  1:2915146/1‑238 (MQ=255)
                                                                                                                                                                                                                      cc‑ggccggGGGGCTGA‑CATTGTGGCACG‑GGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGca                      <  1:1869066/239‑1 (MQ=255)
                                                                                                                                                                                                                      cc‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGc                                                                                                            >  1:1729546/1‑153 (MQ=255)
                                                                                                                                                                                                                          gCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATACCGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGGGAAGTACGGGCATGGCCGTGCacaa                   >  2:3496545/1‑239 (MQ=255)
                                                                                                                                                                                                                          gCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACaa                                                                                                       >  2:1840739/1‑155 (MQ=255)
                                                                                                                                                                                                                          gCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGcaca                    >  1:3272232/1‑238 (MQ=255)
                                                                                                                                                                                                                          gCCTGTGGGCTGA‑CATTGGGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACaa                                                                                                       <  1:1840739/155‑1 (MQ=255)
                                                                                                                                                                                                                            cTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATTGCAGTGCacaaca                 >  1:391473/1‑239 (MQ=255)
                                                                                                                                                                                                                              gtgGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGt                                                 >  2:2080432/1‑205 (MQ=255)
                                                                                                                                                                                                                              gtgGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGt                                                 <  1:2080432/205‑1 (MQ=255)
                                                                                                                                                                                                                               tgGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACAtgt              >  2:25753/1‑239 (MQ=255)
                                                                                                                                                                                                                               tgGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGCTCGGGCATGGCAGTGCACAACAtgt              >  1:1763865/1‑239 (MQ=255)
                                                                                                                                                                                                                                     gA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATgtgcgtg         >  1:1456372/1‑238 (MQ=255)
                                                                                                                                                                                                                                     gA‑CATGGGGGCACA‑GGAACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATgtgcgtgc        <  1:1086869/239‑1 (MQ=255)
                                                                                                                                                                                                                                          ttGTGGCACA‑GGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGggctggggctgg                                                                                                                                                                      <  2:284282/77‑1 (MQ=255)
                                                                                                                                                                                                                                          ttGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGggctggggctgg                                                                                                                                                                      >  1:284282/1‑77 (MQ=255)
                                                                                                                                                                                                                                          ttGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACAtgt              >  1:1068940/1‑229 (MQ=255)
                                                                                                                                                                                                                                          ttGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGGGCACAACATGTGCGGTCCGGt    >  2:689702/1‑239 (MQ=255)
                                                                                                                                                                                                                                             tgGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGttt  >  1:1803376/1‑238 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
GCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCC‑GGCCTGTGGGCTGA‑CATTGTGGCACA‑GGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGTTT  >  NC_000913/2073831‑2074303

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