Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,048,596 G→A G32G (GGC→GGT gfcC ← capsule biosynthesis GfcC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,048,5960GA94.5% 179.9 / NA 55G32G (GGC→GGTgfcCcapsule biosynthesis GfcC family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/31);  minor base C (2/1);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTG  >  NC_000913/1048350‑1048837
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
aCCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAg                                                                                                                                                                                                                                               <  2:358918/251‑1 (MQ=255)
    gTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGGCACCAACTGGCCAACATTTTCCACGGGGCCAACGGAAAGCGTCTGGTGGTCCCCAGGc                                                                                                                                                                                                                                             >  2:11205/1‑249 (MQ=255)
         attgagattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGCCAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATaa                                                                                                                                                                                                                                       >  2:124517/1‑250 (MQ=255)
         attgagattgagTAACTGCTGGCGCACGGATTTAATAGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCAACAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGCGCCAACGGAAAGTGTCTGTTGTTCACCCGGCAGATaa                                                                                                                                                                                                                                       >  1:381050/1‑250 (MQ=255)
               attgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGNGTCTGTTGTTCACCAGGCAGATAAATAGTc                                                                                                                                                                                                                                 >  2:217333/1‑250 (MQ=255)
                 tgagTACCTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCAc                                                                                                                                                                                                                               <  2:266947/250‑1 (MQ=255)
                          gctgGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCtt                                                                                                                                                                                                                       <  1:210563/249‑1 (MQ=255)
                             ggCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTgcg                                                                                                                                                                                                                    <  1:213099/249‑1 (MQ=255)
                                        ttAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGACGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATggg                                                                                                                                                                                                        <  2:374436/250‑1 (MQ=255)
                                                       acaTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAg                                                                                                                                                                                        <  1:38437/251‑1 (MQ=255)
                                                        cctctgcttcAGCTTCTGCGTCCCAGGAAGCAAGCCGTGCCATGACATGTTGGTAGTCTGTAAGCGCTTTAGCTGTCGCCGCAGAGTCGGTCAGCAAGGCCCCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGt                                                                                                                                                                                       <  2:20882/243‑1 (MQ=255)
                                                         aTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAATCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGTCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTTTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTa                                                                                                                                                                                      <  1:206344/251‑1 (MQ=255)
                                                         aTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACAGGTTGTTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTGCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTa                                                                                                                                                                                      <  2:147300/251‑1 (MQ=255)
                                                                     ttCTGCTTGCCAGGAAGCAAGCGTTGCCATGTCATGTTGGTAGTCTTTAAGCGCTTTAGCTGCCGGCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAGCTGGACAGCATTTTCCACGGGGCCAACGGAAAGGGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGGAAAGTaaggggtcaatg                                                                                                                                                                           <  2:81135/250‑12 (MQ=255)
                                                                             cccAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGa                                                                                                                                                                   <  1:362703/250‑1 (MQ=255)
                                                                                        agcTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGTGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTAtt                                                                                                                                                       <  1:259157/251‑1 (MQ=255)
                                                                                           ggtgCCATGACATGTTGGTAGGCTTTAAGCGCTGTAGCTGTCGCCGCAGTGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGTTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATTCGACTGTAATTTATTcgt                                                                                                                                                    <  2:413824/250‑3 (MQ=255)
                                                                                            gggCCAGGACATGTTGGTAGTCTTTAAGCGNTTTAGCTTTCGCCGCTGAGTCGGTCAGCAAGGCGCCAGGCCACCAGAGGCGATCGCGCAGTTGCGGGTGGGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACTGCTGTAATTTATTCAt                                                                                                                                                    <  2:224538/248‑1 (MQ=255)
                                                                                             ggCCATGACATGTTGGTAGTCTTTAAGCGCTNTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATgg                                                                                                                                                  <  2:233763/250‑1 (MQ=255)
                                                                                              gCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATggg                                                                                                                                                 >  1:130486/1‑250 (MQ=255)
                                                                                                   gACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCggg                                                                                                                                            >  1:177308/1‑251 (MQ=255)
                                                                                                   gACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGACACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAAAAGACACCGTTCCTTGAGCAAAGGCATGGGGTGTCATTACGTAAAGTACGGTGGCACTGAAATACGACTGTAATTTATTCATGGTGCggg                                                                                                                                            >  2:55045/1‑251 (MQ=255)
                                                                                                     cATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGCGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCCACTGGACAACATTTTCCACGGGGCCAACGGCAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGGGTCCGTACGTAAAGAACGCTGGCAAAGACATACGACTGCAATTTATTCATGGTGCGGGtt                                                                                                                                          >  2:60905/1‑251 (MQ=255)
                                                                                                      aTGTTGGTAGTCTTTAAGCGCTTTAGCTTTCNCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGttt                                                                                                                                         <  2:215312/251‑1 (MQ=255)
                                                                                                      aTGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAGCGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGttt                                                                                                                                         <  1:312520/251‑1 (MQ=255)
                                                                                                                tCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGACAGCAAGGCACCAGGCCATCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAAAATTTTCCACGGGGCCAACGGAAAGCGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGCTCCTTGCGCAAAAGCATGGGATGTCATTACGTAAAGAACGCTGGCAATGAAATACGACTGTAATTTATTCATGGAGCGGGTTTAAGAAACGTc                                                                                                                               >  1:435729/1‑251 (MQ=255)
                                                                                                                         gcTTTAGCTTTCTCCGCAGAGTCGGTCAGCAAGGCCCCAGGCCTCCAGAGGCGATCGCGCAGTTGCGGTTGTGTTACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACg                                                                                                                      <  1:271811/251‑1 (MQ=255)
                                                                                                                             tAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGa                                                                                                                   >  1:302054/1‑250 (MQ=255)
                                                                                                                                  ttCGCCGCAGAGTCGGTCAGCAAGGCACCCGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGc                                                                                                                                                                                                             >  1:138914/1‑155 (MQ=255)
                                                                                                                                  ttCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGc                                                                                                                                                                                                             <  2:138914/155‑1 (MQ=255)
                                                                                                                                        gcAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTgcg                                                                                                       <  1:339547/251‑1 (MQ=255)
                                                                                                                                                     gCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGcc                                                                                          >  1:268679/1‑251 (MQ=255)
                                                                                                                                                           cACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACttt                                                                                    <  1:326391/251‑1 (MQ=255)
                                                                                                                                                               aGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGAGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACgggg                                                                                                                   >  1:449414/1‑215 (MQ=255)
                                                                                                                                                               aGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTGCCACGGTGCCAACGGAAAGTGTCTGTTGTGCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGa                                                                                                                   <  2:449414/216‑1 (MQ=255)
                                                                                                                                                               aGGCCACCAGAGGCGATCGCGCAGGTGCGGGTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCATATAAATAGGCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACACTTGCCGAATTTGACg                                                                                <  2:166296/251‑1 (MQ=255)
                                                                                                                                                                        gagGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCCCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGGGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGCCGAACCTGGcc                                                                       >  2:381420/1‑251 (MQ=255)
                                                                                                                                                                        gagGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGcc                                                                       <  1:358080/251‑1 (MQ=255)
                                                                                                                                                                            cGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAATTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATACCCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTa                                                                   <  1:257231/251‑1 (MQ=255)
                                                                                                                                                                            cGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCAACGGGGCCAACGGAAAGTGTCTGATGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGCCTTTGACGAACCTGGCCAGt                                                                    >  2:410737/1‑250 (MQ=255)
                                                                                                                                                                                   gcAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTCGACGAACCTGGCCAGTAGCCGAAt                                                            >  2:108271/1‑251 (MQ=255)
                                                                                                                                                                                    cAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGGTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATc                                                           <  2:438423/251‑1 (MQ=255)
                                                                                                                                                                                                     ccaacgtgACAACAGTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCATAATGTGtttt                                          <  1:118303/244‑1 (MQ=255)
                                                                                                                                                                                                     ccAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGtttt                                          >  1:297719/1‑251 (MQ=255)
                                                                                                                                                                                                     ccAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGtttt                                          <  1:6022/251‑1 (MQ=255)
                                                                                                                                                                                                     ccAACTGGACAACATTTTCCACGGGGCCAACGGAAAGGGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGtttt                                          >  1:200320/1‑251 (MQ=255)
                                                                                                                                                                                                        aCTGGACAACATTTTCCACGGGGCCAACGGCAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGAGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGGTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGCCTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGcc                                       >  2:64939/1‑251 (MQ=255)
                                                                                                                                                                                                            gACAACATTTTCCACGGGGCCCACGGAAAGTGTCTGTTGTTCACCAGGCAGATAACTAGGCATCGTTCCTTGCGCAAAGGCATGGGGTGTAATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCGTGGTGCGGGTTTAAGAAACGTCATGTCTACGGGAATAACCCCTGCGCCTAACAGTTGCCGACTTTGGCGAACCTGGCCAGTAGCCGAATAGACACAGAATGTGTTTTGCCACTg                                   <  2:34855/251‑1 (MQ=255)
                                                                                                                                                                                                                 caTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgt                              <  1:97836‑M1/251‑6 (MQ=255)
                                                                                                                                                                                                                                 aaCGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaacc                 >  2:347610‑M1/1‑230 (MQ=255)
                                                                                                                                                                                                                                      aaaGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagtt         >  1:351814‑M1/1‑225 (MQ=255)
                                                                                                                                                                                                                                            tCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCCTGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgttgactgagtaaccagacagttgatgtg   >  2:269466‑M1/1‑219 (MQ=255)
                                                                                                                                                                                                                                              tgttgtTCCCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATGCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGattgtggactgagtaaccagacagttgatgtgc  <  2:161038‑M1/250‑34 (MQ=255)
                                                                                                                                                                                                                                                 tgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACtt                                                                                     <  1:201846/164‑1 (MQ=255)
                                                                                                                                                                                                                                                 tgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACtt                                                                                     >  2:201846/1‑164 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
ACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTGGCAATGAAATACGACTGTAATTTATTCATGGTGCGGGTTTAAGAAACGTCATTTCTACGGGAATAACCCCTGCGCCTAACATTTGCCGACTTTGACGAACCTGGCCAGTAGCCGAATCGACCCAGAATGTGTTTTGCCACTGTTTTTCTGGCCGCGTGGAAGTCACCTCTTCCTG  >  NC_000913/1048350‑1048837

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