Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,074,021 A→C 13.3% Y828S (TAC→TCC)  flu → CP4‑44 prophage, antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,074,0210AC13.3% 169.8 / 3.3 66Y828S (TAC→TCC) fluCP4‑44 prophage, antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base A (47/8);  new base C (3/8);  total (50/16)
Fisher's exact test for biased strand distribution p-value = 2.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTA  >  NC_000913/2073790‑2074258
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
gtcccgtcGGTCCCGTTAAAAATGGGGGTTTTGCCCGTGGGGCCCGGCCGGAAAGCAGCGGCAGGTAGGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAAt                                                                                                                                                                                                                                            <  1:607726/233‑1 (MQ=255)
 gcctccgcgTCACGATAACCATGGCGGTATTGCCCGTGGGGCCAGGCCGAAAAGGAGCGGCAGGTATGGCTTGGTCCGTCGGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGGCCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAAt                                                                                                                                                                                                                                            <  1:261346/230‑1 (MQ=255)
           cACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGGATATGGTGCTGCTGGCCATTCTTCCGTTTATGTTAAGGATGATGACGGCTCCCGTGCCGGCCCGGGCCGGGGTGATGCCGGCAGCCTGGGC‑GGGACCCTGAAACTGGTaccac                                                                                                                                                                                                                                  >  2:1216709/1‑235 (MQ=255)
                  aCAATGGCGGTATTGCCCGTGGGGCAACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGtcctcc                                                                                                                                                                                                                          <  2:1150821/239‑1 (MQ=255)
                       gggggttgtGCCCGGGTGGCCACGCCGGAAAGCAGCGTGAGCTATGGATTCGTCCGTCGTGAGGGTGACCTGATGAGCACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATTACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCt                                                                                                                                                                                                                     <  2:1229799/231‑1 (MQ=255)
                            tATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTgg                                                                                                                                                                                                                 <  2:430521/238‑1 (MQ=255)
                                gCCCGGGGGGCCACGCCGGAAAGCGGCAGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGa                                                                                                                                                                                                            <  2:614428/239‑1 (MQ=255)
                                     tGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGACCACACGTCCTCCGGCCTGTGGGCTGACATtg                                                                                                                                                                                                       >  2:466273/1‑239 (MQ=255)
                                                  gAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCCCAGGGaa                                                                                                                                                                                           >  2:362717/1‑238 (MQ=255)
                                                                          gtccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGCCATTTTGGCCCAGGGAACCCGCCCCAG‑CATGA‑AAGCGTAATc                                                                                                                                                                >  2:253200/1‑239 (MQ=255)
                                                                          gtccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAACGCCATc                                                                                                                                                                >  2:1151513/1‑239 (MQ=255)
                                                                          gtccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGGGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCacc                                                                                                                                                                >  2:1297291/1‑237 (MQ=255)
                                                                          gtccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCATCCGGCCTGTGGGCGGACATTGTGGCACAGGGAACCCGCCCCAG‑CATGA‑ACGCGGCATc                                                                                                                                                                >  2:750738/1‑239 (MQ=255)
                                                                                     gggTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACCCGTCCTCCGGCCTTTGGGCTGACCTGGTGGCACCGGAAACCCGCCACAG‑CATGA‑AAACGTCACCGGACAATAACg                                                                                                                                                     >  2:1246747/1‑239 (MQ=255)
                                                                                        tGACCTGATGCGTACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGCTGATGACGGCTCCCGTGCCGGCTCGGTCCGGCATGATGCCGGCAGCCTGTGC‑GGCTCCCTGAATCTGGTCCACGAGTCCTCCGGTCGGCGGGTCGGCATTCAGGCACCGGGGAACCGCCCCCG‑CCTGA‑AAGCGGCCTCgcaaaaaaaaaca                                                                                                                                                   >  1:1100817/1‑226 (MQ=255)
                                                                                        tGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGAGGATGACGGCTCCCGTGCCGGCACGGGCCGGCGTTATGACGGCAGACTGGGC‑GGCTACCTGAATCGGGAACACACATCCTCCGGCCTGCAGGCTGACCTTGTGGCACAGGGCCCCCCCCCCCC‑CATGA‑AAGCAGCCTCGGACAATAACGaa                                                                                                                                                   >  1:942765/1‑237 (MQ=255)
                                                                                                tgaGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTcccggc                                                                                                                                           >  2:1096356/1‑234 (MQ=255)
                                                                                                         agGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGGTGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCGGAATCTGGAACACACGTCCTCCGGCCTGGGGGGTGCCATTGTGGCACAGGGAACCCGCCCACG‑AATGA‑AAGAGTCATCGGACAAAAACAAATTCCCCGCCCggggcggg                                                                                                                                 >  1:1099792/1‑239 (MQ=255)
                                                                                                         agGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCGGACATTGTGGCACAGGGAACCCGCCACAG‑CATGa                                                                                                                                                                           >  2:1011254/1‑198 (MQ=255)
                                                                                                              gCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGCACTCGGACCATAACGACTTCCGCGCCCGgggctggggctg                                                                                                                            >  2:1217050/1‑239 (MQ=255)
                                                                                                                         tgtgACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCCGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGGACACACGTCCTCCGGCCTGCGGGCCGACATTGTGGCCCACGGAACCCGCCACCG‑CCGGC‑AAACGGCATCGGACCATAACGAACTCCCCGCCCggggcccggaccggtcgggccgac                                                                                                                 >  1:1056089/1‑220 (MQ=255)
                                                                                                                                 cgGGGGTATATGGTGTTGCTGGCCATTCTTCCGTTGATGGTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGCATCTGGTCCCCACGCCATCCGGCCTGTGGGCTGACCTTTTGGCACAGGGAAACCGCCACA‑‑CACGACAAGCGGAATCGGCCAATACCGCCATCCGCGCCCcgggctggtgcgtcctggcgccgcctgaaaa                                                                                                          >  1:1199810/1‑218 (MQ=255)
                                                                                                                                 cgGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGa                                                                                                                                                                                                                                              >  2:296502/1‑108 (MQ=255)
                                                                                                                                 cgGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGa                                                                                                                                                                                                                                              <  1:296502/108‑1 (MQ=255)
                                                                                                                                 cgGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCCCGCCCGGGGCTGGGGCTGGCTGGGGTCAC‑TGGAAAcc                                                                                                        >  2:519149/1‑239 (MQ=255)
                                                                                                                                   ggggTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCCCAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCg                                                                                                       >  2:857574/1‑238 (MQ=255)
                                                                                                                                              gactgctgGCCAGTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGTAAGCCTGGGC‑TGTTCACTGAATCTGGTACCCCCGTCCTCCGGCCTGTAGCGTGGCCTTGGTGCCCAAGGAACCC‑CCACCGTCCTGA‑CAGCGTCATTGGACAATACCGACTTCCGCCACCGGGGCTGGGCCGTCCTGGGCCCAC‑tcaaaaagagaccgccctgca                                                                                           >  1:886275/3‑219 (MQ=255)
                                                                                                                                                ctgctgGCCATTCTTCCGTTGCTGTTAAGGATGCTGACGGCTCCCGTGCAGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGCTACCTGATTCTGGTACTCCCGTCTTCCGGCCTGGGGGCTGACATTGTGTCACCGGAAAACCGCCACCG‑CATGT‑CAGGTCCATTGGCCACTACAGCCTTCCGCGCCCGGGGCTGGGCCGGGCCGGGCACACGTGGaacaccatacacctcccga                                                                                          >  1:200571/1‑222 (MQ=255)
                                                                                                                                                ctgctgGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGc                                                                                                                                                                                                               >  2:965833/1‑124 (MQ=255)
                                                                                                                                                ctgctgGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGGAACCGGGCTGCCCTTCAGt                                                                                         >  2:744703/1‑239 (MQ=255)
                                                                                                                                                ctgctgGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGt                                                                                         >  2:61627/1‑239 (MQ=255)
                                                                                                                                                ctgctgGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCGCCCGTGCCGGAACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAGTCTGGCACACACGGCCTACGGACTGGGGGCTGCCATTGGGGCACAGGGAACCCGCCACAC‑CCTTA‑ACAGGCCATCGGACAATAGCGACTTCCGAGCCCTGGGCCGGGGCAGGCCGGGGGCCG‑TGGAAACCGGGCGGCCCtgtagg                                                                                         >  1:1325429/1‑234 (MQ=255)
                                                                                                                                                ctgctgGCCAGTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGc                                                                                                                                                                                                               <  1:965833/124‑1 (MQ=255)
                                                                                                                                                ctgctgGCAATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCCCAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGGCAATAACGACTTCCGAGACCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAAccctgtctcttataac                                                                                          >  2:321834/1‑224 (MQ=255)
                                                                                                                                                   ctgGCCATTCTTCCGTTGATGTTAAGGATGCTGACGGCTCCCGTGCCGGCACGGTACGGGATGATTCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGCAATTGTGGCACAGTGGAACTGCCACAG‑AATCA‑AACCGTGAACGGACAATACCGAATTCCCCGCCCGCGGCCCGGGCTGCGCGGGCTCCC‑TTGAAAcccatcgccccctcagaaca                                                                                      >  1:1280926/1‑221 (MQ=255)
                                                                                                                                                   ctgGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACCGGGAACCCGCCCCAG‑CCTGA‑AAGCGGCACCGGCCAATAACAACTTCCCCGCCCGGGGCCGGGGCGGGCCGGGCCCAC‑TGAAAACCggccctccccccagagt                                                                                       >  1:1152584/1‑223 (MQ=255)
                                                                                                                                                   ctgGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CCTGA‑AAGCGTCATCGGCCAATAACGACTTCCGCGCCCCGGCCTGGGGCGGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTCCAGGAt                                                                                       >  1:894823/1‑238 (MQ=255)
                                                                                                                                                   ctgGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGCCAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCCCAC‑TGGAAACCGGTTTGCCCTTCAGTAt                                                                                       >  2:568496/1‑238 (MQ=255)
                                                                                                                                                   ctgGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCCCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAgatc                                                                                       >  2:476563/1‑235 (MQ=255)
                                                                                                                                                   ctgGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGGCATCGGACCATAACGACTTCCGCGCCCGGGGCTGGGCCTGGCTGGGCCCAC‑TGGAAACCGGTCTGCCCTTCAGTATc                                                                                      >  2:566596/1‑239 (MQ=255)
                                                                                                                                                      gCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGGCATTTTGGCACAGGGAACCCGCCACAG‑CCTGA‑AAGCGTCCTCGGACCATAACGACTTCCGCCCCCGGGGCTGGGCATGCCTGGGCTCAC‑TGGAACCCGGTCGCCCCTTccggtcaacc                                                                                   >  1:967356/1‑230 (MQ=255)
                                                                                                                                                        cATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCGGGCATTTGGGCACAGGGAACCCGCCACAG‑CATGA‑AACCGGCATCGGACAATACCGCCTTCCACGCCCGGGGCTGGGGCTGGCTGGGCCCCC‑CGGAAACCGGTCTGcccccaaggacacctga                                                                                 >  1:35425/1‑225 (MQ=255)
                                                                                                                                                                             gtggatggggccccccggcccGGCGCGGTCGGGGAGGGCGCCGGCAGCGGGGGC‑GGATCCATGAATCTGGTACACACGTCCCCGGGCTTGTGGGATGACATTGTGCCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCca‑c                                                            <  2:1099792/219‑1 (MQ=255)
                                                                                                                                                                                gatgaCGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑GGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATCCTTGAGCCC‑CAAcg                                                        >  2:920844/1‑238 (MQ=255)
                                                                                                                                                                                 atgaCGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGCCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGCCAACCTGATGCTGGAGCCA‑CAACTg                                                       >  2:183881/1‑239 (MQ=255)
                                                                                                                                                                                  tgaCGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCCGGGGCTGGCTGGGCTCAC‑TGGAAACCGTGCTGCCCTTCAGAATCACTGACAACCTGATACTGGAGCCA‑CAACCGc                                                      >  1:1157307/1‑239 (MQ=255)
                                                                                                                                                                                    ggggCTCCCGGGCCGGCCGGGTGCGGGATGATGGCGGCACGCGGGGC‑GGATCCCGGAACCTGGTCCACCGTTCCCCGGGCGGGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCa                                                     <  2:481866/236‑1 (MQ=255)
                                                                                                                                                                                        cTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACCCACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCC‑CAACTGAGGtat                                                 >  2:1084984/1‑238 (MQ=255)
                                                                                                                                                                                         ccccGTGCCGGCACGGTCCGGGATGAGGCCGGCAGCCGGGGG‑GGATCCCTGAATCTGGTACACCCGTCCCCCGGCCTGTGGGCTGACTTGGGGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCAGTATAc                                               <  2:1152584/238‑1 (MQ=255)
                                                                                                                                                                                               gccggccgGGTCGGGGATGATGCCGGCACCCTGGGC‑GAATCCCTGAATCTGGACCCCACGCCCGCGGTCTTGGGGGTTAAAATGGGGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCAGTATACCTGGc                                          <  2:35425/238‑1 (MQ=255)
                                                                                                                                                                                               gCCGGCACGGGCGGGAATGGGGCCGACAG‑CTGGGCGGGATCCCTGATCCGTGTACACAGGTCCCCCGGCTTGGGGGCTGACATTGTGGCACAGGGACCCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCAGTATACCTGGc                                          <  2:843427/238‑1 (MQ=255)
                                                                                                                                                                                                   gCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTCCCGTATACCTGGCAggg                                      >  2:427716/1‑238 (MQ=255)
                                                                                                                                                                                                              gatgatGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGCCTGCCCTTCGGTATCACTGGCAACCTGATGCCGGAGCAC‑CAACTGCGGGATACCTGGCAGGGGCTTTCCCTGGa                          >  1:768202/1‑239 (MQ=255)
                                                                                                                                                                                                               aggaggCCGGGACCCGGGGC‑GAATCCCGGAACCGGGAACAACCGTCCTAGGGCCTGTGGTGTGACTTGGGGGCAGGGGGACCCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCAGTATACCTGGCAGGGACTTTCCCTGGa                          <  2:138118/233‑1 (MQ=255)
                                                                                                                                                                                                                  atgCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCCGTATCCCTGACAACCTGATGCTGGAGCCA‑CAACTCCAGGTAACCTGGCAGGGCCTTTCCCCGGGTGag                      >  1:1251795/1‑238 (MQ=255)
                                                                                                                                                                                                                   tgCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCCC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCCGATGCTGGAGCCA‑CAACTGCAGTATACCTGCCAGGGACTTTCCCTGGATGACg                     >  2:549897/1‑239 (MQ=255)
                                                                                                                                                                                                                   tgCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCCCGGACAACCTGATGCTGGAGCCA‑CAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACg                     >  2:996829/1‑239 (MQ=255)
                                                                                                                                                                                                                   tgCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCAGtata                                                >  2:717872/1‑212 (MQ=255)
                                                                                                                                                                                                                   tgCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCGGGGCTCAC‑TGAAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCCGTATACCTGGCAAGGACTTACCCAgggggagg                     >  1:1200173/1‑233 (MQ=255)
                                                                                                                                                                                                                   tgCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTACCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTGTGCCCTTCAGTCTCCCTGGCAACCTGCTGCTCGGGCGC‑CAAACGCAGGAGACCTGGCCGGGACGTTCCCCGGATGACg                     >  1:115343/1‑239 (MQ=255)
                                                                                                                                                                                                                   tgCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCACCCCAG‑CATGA‑AAGCGTCATCGGAAAATAACGACTTCCGCGCCCGGGGCTCGGACTGGCTGGGCCCAC‑TGGCAACCGGTCTGACCTTCCGGATTCCTGACAACCTGATGCCGACGGCC‑CCACGCCAAGATACCTGGCCGGGCCTTTACCTgtgggcgg                     >  1:89628/1‑232 (MQ=255)
                                                                                                                                                                                                                   tgCCGGCAGCCGGGGC‑GGATACCTGAATCGGGTCCACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCGGGGGCTGGCGGGGCTCAC‑TGGAACACGGTCTGCCCTTCAGTATACCTGCCAACCTGATGCTGG‑GCCACCCACGCCACTATACCTGCCAGGCGCTTTCCCCGGATGaag                     >  1:933513/1‑237 (MQ=255)
                                                                                                                                                                                                                        gCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCGGGCTGGGCTCAC‑TGGAAACAGGTCTGCCCTTCAGTATCAATGCCAACCTGTTGCTGGAGCCA‑CAACTCCCGAATACCGTGCAGGGGACTTCCCCGGATGAAGGAAAg                >  1:297838/1‑239 (MQ=255)
                                                                                                                                                                                                                               gggg‑ggATCCCTAAACTTGGTACCCCCGTCCCCGGGCCTGGGGGCTGACTTGGGGGCACAGGAACCCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAAc          <  2:1251795/238‑1 (MQ=255)
                                                                                                                                                                                                                                 gC‑GGATCCCGTAATTTGGTACACCCGCCGTCGGGCTGGGGGGGTGACTTGGTGGCACGGGGACCCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGc        <  2:942765/238‑1 (MQ=255)
                                                                                                                                                                                                                                       tACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCAGTATACCTGGCAGGGACTTTCCCTGGCGGAGGGAAAGGACAAAACCGGGTa  >  1:1196786/1‑239 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                 
GTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGC‑GGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACATTGTGGCACAGGGAACCCGCCACAG‑CATGA‑AAGCGTCATCGGACAATAACGACTTCCGCGCCCGGGGCTGGGGCTGGCTGGGCTCAC‑TGGAAACCGGTCTGCCCTTCAGTATCACTGACAACCTGATGCTGGAGCCA‑CAACTGCAGTATACCTGGCAGGGACTTTCCCTGGATGACGGTAAGGACAACGCCGGTTA  >  NC_000913/2073790‑2074258

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