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.

GCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGTTT  >  NC_000913/2073831‑2074303
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
gCCACGCCGGAAAGCAGCGGCATCTATGGATTCGTCCGTCGGGAGGGTGGCCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTgg                                                                                                                                                                                                                                              <  2:2561226/238‑1 (MQ=255)
gCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTgg                                                                                                                                                                                                                                              <  1:2417607/238‑1 (MQ=255)
   aCGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGcaca                                                                                                                                                                                                                                          <  1:3467269/239‑1 (MQ=255)
     gCCGGAAAGCAGCGGACGTTTTGGCTTCGTCCGTCTGGAGGGTGACCTGAGGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAgg                                                                                                                                                                                                                                       <  2:1346564/239‑1 (MQ=255)
                 cGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACCCACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCcaca                                                                                                                                                                                                                           >  1:1988406/1‑239 (MQ=255)
                                                                            ggTCTGTGGCCGGGGGGGTATATGGTGCTGCGGGCCTTTCTTCCGTTGATGTTAGGGATGATGCGGGCTCCCGGGCCGGCACGGTCGGGGAGGAGGCGGGCAGCCGGGGCGGATCCCTGAATCGGGTACACGCGCCCTCCGGCTTGTGGGCTGACTTTGGGGCACGGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGcc                                                                                                                                                                                     <  2:318416/217‑1 (MQ=255)
                                                                                                    gTGCTGCTGACCTTTCGTCCGTTGATGAGAAGGATGACGCCGGCCCCAGCGCCGGGCGGGTCGGGGAGGTCCCCGGCAGCCGGGGGGTAACCCTAAATCTGGACCCCCCGCCCCCCGGCTGGGGGGCTGACATTGGGGCACGGGAACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCtt                                                                                                                                         <  2:240419/238‑1 (MQ=255)
                                                                                                                              gTTAAGGATGAGGACGGCCCCCGGGCCGGCACGGCCGGGGATGATGCCGGCAGCCGGGGCGGATACCTGAATCTGGTACACACGCCCCCCGGCCGGGGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCa                                                                                                                                                                                                             <  2:3450741/144‑1 (MQ=255)
                                                                                                                                      ggggggcGGCTCGCGGGCCGGCACGGCCGGGGAGGGTGCCGGCGGCGGGGGCGAATCCCTGAATCGGGCCCACCCGTCCCCCGGCGGGGGGGCTGACTTGGGGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACaa                                                                                                       <  2:406845/232‑1 (MQ=255)
                                                                                                                                           aCGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGcaca                                                                                                                                                                                                                                          <  2:133331/103‑1 (MQ=255)
                                                                                                                                              gcccccGCGCCGGCACGGTCCGGGATGATGCCGGCACCCTGGGCGGATACCTGAATCTGGTCCACACGTCCTCCGGCCTGTGGGCTGACATTGCGGCACAGGGAACCCGCCACAGCATGAAAGCGACATCGGACACTAACGACTTCCGCGCCCGGGGCTGGGGCCGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGGATCACTTACCACCTGATGCTGGAGTCACCACt                                                                                                     >  2:3307795/4‑232 (MQ=255)
                                                                                                                                                           cccggccggggATGATGCGGCCACGCGGGGGGAGTACCGGAGTCTGGACCAGACTGCTCCGGGCCGGGGGGCGGACATTGGGGAACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGAc                                                                                 <  2:1448124/231‑1 (MQ=255)
                                                                                                                                                                tcggggATGATGCGGGCAGCCTGGGCGGATCCCTGAATCTGGTACACACGTCCTCCGGCCTGGGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTg                                                                                                                                              <  1:1403610/170‑1 (MQ=255)
                                                                                                                                                                       tgaggCGGGCAGCCTGGGCGGATACCTGAATCGGGTACACCCGTCCTCCGGCCTGGGGGCTGACTTTGTGGCCCAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGAt                                                                      <  1:982981/238‑1 (MQ=255)
                                                                                                                                                                        ggggccggccgCGGGGGCGGATACCTGAAGCTGGTACCCACGTCCTCCGGCTTGTGGGCTGACATGGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATg                                                                     <  1:687882/235‑1 (MQ=255)
                                                                                                                                                                            ccGGCAGCCTGGGCGGATACCTGAATCTGGTACCCACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAATGGAAACCGGTCGGCCCTCCAGCATCACTGACCACCCGATGCTGGAGCCCCAACTGCAGTATACCTCGCAGGGACATTCCATGGCAGACgat                                                                >  2:407643/1‑237 (MQ=255)
                                                                                                                                                                             cggcggCCGGGGGGGATACCGGAATCGGGTACACACGTCCCCCGGCCTGTGGGCTGACATTGGGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGt                                                                <  1:2474505/238‑1 (MQ=255)
                                                                                                                                                                             cGGAAGCCGGGGCGGATACCTGATTTGGGTCCACACGTCCTCCGGCTTGTGGGCTGACTTTGGGGCACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGt                                                                <  1:2406601/238‑1 (MQ=255)
                                                                                                                                                                                    cggggCGGATACCAGAATCTGGTACCCCCGTCCTCCGGCCGGTGGGCTGACATTGGGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCCTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGAc                                                          <  1:2234035/235‑1 (MQ=255)
                                                                                                                                                                                        ggggATCCCTGACTTTGATCCACCCGTCCTCGGGCCTGGGGGCGAACATTGGGGCACGGGAAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGc                                                     <  1:1266471/236‑1 (MQ=255)
                                                                                                                                                                                        ggggACACCTGAACCTGTAACCCACTCCCCCCGGCGTGTGGGCTGACATTGGGCAACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGc                                                     <  1:3132987/236‑1 (MQ=255)
                                                                                                                                                                                        gCGGATACCTAAATCTGGTACACCCGCCCCCGGGCCTGTGGGCTAACTTTGGGGCACGGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGc                                                     <  1:1967024/238‑1 (MQ=255)
                                                                                                                                                                                              cccggAATCTGGTACACCCGCCCTCGGGCCTGGGGGCTGACTTGGTGGCACGGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATg                                                                     <  1:2729814/212‑1 (MQ=255)
                                                                                                                                                                                                 ggattgtGGTCCACCCGCCCCCGGGCCTGTGGTTTGACTTGGGGGCACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTAtgtg                                           <  1:1054021/233‑1 (MQ=255)
                                                                                                                                                                                                     gtggacacacccGCCCCCCGGCGTGGGGTCAGAATTTGGGCCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGt                                       <  2:174268/233‑1 (MQ=255)
                                                                                                                                                                                                      cTGGTAAACACGCCCTCCGGCCTGTGGGCTGACATTCTGGCACAGGGAACCCGCCACAGCATGAAAGCGCCCTCGGACAATAACGACTTCCGCGCCCAGTGCTGGGGCTGGGTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCACCAACTGCAGTATACCAGGCAGGGACTTCCCCTGGATGACGGTAAGGACAACGCTGATTATGTGAAGt                                       >  2:877994/1‑238 (MQ=255)
                                                                                                                                                                                                          tACACACGACCTCCGGCCTGTGGGCTGACATTGCGGCACAGGGAACCCGCCACAGCATGAAATCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCGGGGGCTGGCTGGGCTCACCGGAAACCGGTCTTCCCTCCAGTATCACTGACAACATGAGGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTACGGACAACGCCGGTCATGGGAAGTTCggg                                  >  2:1430549/1‑239 (MQ=255)
                                                                                                                                                                                                            cccccgccccccGGCTTGGGGGCTGACTTGGGGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCa                                <  1:166289/229‑1 (MQ=255)
                                                                                                                                                                                                            ccaccGTCCGCCGGCTGGGGGGTGAGCTTTGGGGCGCGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCa                                <  2:3310691/235‑1 (MQ=255)
                                                                                                                                                                                                            cacccGCCCTCCGGCCGGTGGGCGGAAATTGGGCCCAAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCa                                <  1:535815/239‑1 (MQ=255)
                                                                                                                                                                                                               cggTCCCCCGGCCAGTGGGCTAACATGAGGGCACGGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATgg                             <  1:3252153/237‑1 (MQ=255)
                                                                                                                                                                                                               ccgccccccGGCCTGTGGGCTGACATGGGGCCACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATgg                             <  1:242107/232‑1 (MQ=255)
                                                                                                                                                                                                               acGTCCTTCGGACAGTGGGCCGAAATTGCGGAACAGGGAACCCGACACAGCACGAAAGCGTCATCGGACAATAACGACTTCCGCGACCGGGGCTGGGACTGGCTGGGCTCACTGGTAACCGGCCGCCCGTTCAGTATCACTGACAACATGATGCTGGAGCCACAAACGCCGTATACCTGGCCGGGACTTTCCCTGGATGCAGGTAAGGAACACGCCGGTTATGTGCTGTTCgggaaag                              >  2:342881/1‑234 (MQ=255)
                                                                                                                                                                                                                cctcccccGGCCGGGGGGCTGACTTGGGGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGc                            <  1:754077/237‑1 (MQ=255)
                                                                                                                                                                                                                 gTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGa                                                                       >  1:2562448/1‑195 (MQ=255)
                                                                                                                                                                                                                 gTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCa                           >  1:1877286/1‑239 (MQ=255)
                                                                                                                                                                                                                 gTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCa                           >  2:2006741/1‑239 (MQ=255)
                                                                                                                                                                                                                 gTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGATGGGGCTGGCTGGGCTCACTGGaaa                                                                                                                                                      >  1:1055853/1‑116 (MQ=255)
                                                                                                                                                                                                                   cctccGGCTGGGGGGCGGACATTGTGGCACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAg                          <  1:2434749/238‑1 (MQ=255)
                                                                                                                                                                                                                   cctccGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGt                         >  1:3272097/1‑239 (MQ=255)
                                                                                                                                                                                                                    ctccGGCTGGTGGGTCGAATTGGTGCACAGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTg                        <  1:1656526/239‑1 (MQ=255)
                                                                                                                                                                                                                    ctccGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGt                         >  1:2915146/1‑238 (MQ=255)
                                                                                                                                                                                                                      ccggccggGGGGCTGACATTGTGGCACGGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGca                      <  1:1869066/239‑1 (MQ=255)
                                                                                                                                                                                                                      ccGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGc                                                                                                            >  1:1729546/1‑153 (MQ=255)
                                                                                                                                                                                                                          gCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATACCGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGGGAAGTACGGGCATGGCCGTGCacaa                   >  2:3496545/1‑239 (MQ=255)
                                                                                                                                                                                                                          gCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACaa                                                                                                       >  2:1840739/1‑155 (MQ=255)
                                                                                                                                                                                                                          gCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGcaca                    >  1:3272232/1‑238 (MQ=255)
                                                                                                                                                                                                                          gCCTGTGGGCTGACATTGGGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACaa                                                                                                       <  1:1840739/155‑1 (MQ=255)
                                                                                                                                                                                                                            cTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATTGCAGTGCacaaca                 >  1:391473/1‑239 (MQ=255)
                                                                                                                                                                                                                              gtgGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGt                                                 >  2:2080432/1‑205 (MQ=255)
                                                                                                                                                                                                                              gtgGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGt                                                 <  1:2080432/205‑1 (MQ=255)
                                                                                                                                                                                                                               tgGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACAtgt              >  2:25753/1‑239 (MQ=255)
                                                                                                                                                                                                                               tgGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGCTCGGGCATGGCAGTGCACAACAtgt              >  1:1763865/1‑239 (MQ=255)
                                                                                                                                                                                                                                     gACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATgtgcgtg         >  1:1456372/1‑238 (MQ=255)
                                                                                                                                                                                                                                     gACATGGGGGCACAGGAACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATgtgcgtgc        <  1:1086869/239‑1 (MQ=255)
                                                                                                                                                                                                                                          ttGTGGCACAGGGACCCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGggctggggctgg                                                                                                                                                                      <  2:284282/77‑1 (MQ=255)
                                                                                                                                                                                                                                          ttGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGggctggggctgg                                                                                                                                                                      >  1:284282/1‑77 (MQ=255)
                                                                                                                                                                                                                                          ttGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACAtgt              >  1:1068940/1‑229 (MQ=255)
                                                                                                                                                                                                                                          ttGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGGGCACAACATGTGCGGTCCGGt    >  2:689702/1‑239 (MQ=255)
                                                                                                                                                                                                                                             tgGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGttt  >  1:1803376/1‑238 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              
GCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAGCATGAAAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCACTGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCACAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTATGTGAAGTTCGGGCATGGCAGTGCACAACATGTGCGTGCCGGTTT  >  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: