Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,160 G→C 5.9% intergenic (‑143/+28) yedQ ← / ← yodC putative membrane‑anchored diguanylate cyclase/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,028,1600GC 5.9% 285.3 / 7.9 85intergenic (‑143/+28)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (37/41);  new base C (2/3);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTT  >  NC_000913/2027887‑2028427
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                           
tCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGa                                                                                                                                                                                                                                                                            >  2:143312/1‑275 (MQ=255)
 cgagCCAGGTGGGCAAAGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTGTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCAGAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                                                                                                                                                                                           <  2:411188/272‑1 (MQ=255)
     ccAGGGGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc                                                                                                                                                                                                                                                                          <  1:86162/272‑1 (MQ=255)
      cAGGTGGGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                                                                                                                                                                                                      <  2:170725/275‑1 (MQ=255)
                  gTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCAACTATCATTACTTGCAGGCATTAGCCCACCAGGCCGGATGTGATCCGGCCTGGTGTTaa                                                                                                                                                                                                                                                            >  2:428336/1‑273 (MQ=255)
                                                  aTTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTCTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGTTTCTGGCATAAACTGGCGATACAGTACCGCAAAACACCGATAATTTTATCATGGTTGCATTATAAGACGCTACTATCCTTCATGGCAGGCATAAGCACCCAAGTCCGGAGGTGATCCGGCCTGCGGTAAAAGATACCCTTTACAACTTTTTCCgaaaacgcatg                                                                                                                                                                                                                          >  4:312097/1‑265 (MQ=255)
                                                                                 gCGAGTTTTTTCAACCAGCTCTGGTTTGCCATTTTTGTCGCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                                                                                                                                                                                           <  2:416020/195‑1 (MQ=255)
                                                                                   gAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTTTCTCGTGCTGCACCCTGACGACTCTCAACCTTCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAACACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTTATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                                                                                                                                                                                         <  2:312101/275‑1 (MQ=255)
                                                                                        ttttCGACCGACTCCGGTTGTCCAATTTTGTGTCGTGATTACCGCTGACGAATCCCAACCTGCGTTTGTATGATTCGGGCAGACCCGGGCGACAGCGTACCGCAAAACATCGATAATGTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt                                                                                                                                                                                     <  4:450557/274‑1 (MQ=255)
                                                                                                        atatgtattttttCTCACGTGCTCTCCCNNNACGAAACACAGCCGACTTTTTTATGTGTATGGCAGACCCTGGCGATAGCGTACCGCAAAAAAGCGATAATTGTATAATGGTGGCAAGATAAGACGCGACTATCATTCATTGCAGTCAGTAGCACACCAGGCCGGATGTGATCCGGCCTGCTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCAGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGt                                                                                                                                                                     <  3:147185/266‑1 (MQ=255)
                                                                                                             ccATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                                                                                                                                                                                             >  4:438996/1‑245 (MQ=255)
                                                                                                             ccATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCAATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCg                                                                                                                                                                >  2:103205/1‑274 (MQ=255)
                                                                                                                 ttttGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                                                                                                                                                                                         <  2:439001/245‑1 (MQ=255)
                                                                                                                          gTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGAAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGaca                                                                                                                                                  <  2:407500/275‑1 (MQ=255)
                                                                                                                           tgctgcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAATTCTTCCGCATAACGACTCCCCTCCCCCGGACCCAGCTCGGTTTCATGAAAAGCTTCCCCTTTGACCCCGTAACCGTCATAcctttctctt                                                                                                                                                 >  1:147185/1‑267 (MQ=255)
                                                                                                                               gcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                                                                                                                                                                                           >  6:317748/1‑149 (MQ=255)
                                                                                                                                cACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc                                                                                                                                                                                                                                                                          <  6:21749/149‑1 (MQ=255)
                                                                                                                                  cccTGACGCCTCCCCACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCa                                                                                                                                          <  2:47214/275‑1 (MQ=255)
                                                                                                                                  cccTGACGACTCCCAACCTTCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCa                                                                                                                                          <  2:19720/275‑1 (MQ=255)
                                                                                                                                   ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGa                                                                                                                                                                                                                                                                            >  5:189126/1‑144 (MQ=255)
                                                                                                                                   ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGa                                                                                                                                                                                                                                                                            <  6:189126/144‑1 (MQ=255)
                                                                                                                                   ccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCgg                                                                                                                                                                                                                                                                       <  5:207653/149‑1 (MQ=255)
                                                                                                                                     tGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACTCTTTCACACTTCTTCCGCAGAATGTCTCCCCTCCCCAGGAACAAGCTCGGTTTCATGAAAAGTTTccc                                                                                                                                                                                 >  2:161134/1‑233 (MQ=255)
                                                                                                                                        cgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                                                                                                                                                                                                  <  5:76276/149‑1 (MQ=255)
                                                                                                                                         gacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                                                                                                                                                                                                 >  6:385435/1‑149 (MQ=255)
                                                                                                                                          acTCCCAACCTGCTTTTTTATGATTATGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATACTTTTATAATGGTGGCATGCTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGCACCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGACCGCCTCCCCCCCCCCGGAACAAGCTCGGTTTCCTGAAAAGCTTCCCGCTTGACCCCGTAACCGCCAAACCAACAACCCCCAAcccaacc                                                                                                                                      >  4:497663/1‑266 (MQ=255)
                                                                                                                                            tCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                                                                                                                                                                                              >  7:10010/1‑148 (MQ=255)
                                                                                                                                            tCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt                                                                                                                                                                                                                                                              >  6:302696/1‑149 (MQ=255)
                                                                                                                                            tCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGa                                                                                                                                >  2:149195/1‑275 (MQ=255)
                                                                                                                                              ccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa                                                                                                                                                                                                                                                            <  7:4315/149‑1 (MQ=255)
                                                                                                                                              ccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCt                                                                                                                                   <  2:497668/270‑1 (MQ=255)
                                                                                                                                                aaCCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtgctaa                                                                                                                                                                                                                                                          <  8:10010/149‑6 (MQ=39)
                                                                                                                                                aaCCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                                                                                                                                                                                          >  5:32479/1‑149 (MQ=255)
                                                                                                                                                  ccTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGAt                                                                                                                                                                                                                                                        >  7:376264/1‑149 (MQ=255)
                                                                                                                                                    tGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                                                                                                                                                                                                    <  5:321649/135‑1 (MQ=255)
                                                                                                                                                    tGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                                                                                                                                                                                                    >  6:321649/1‑135 (MQ=255)
                                                                                                                                                          tttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCCTGAAAAGCTTCCCGCTTGACCCCGTCACCGTCATACCAACGACCCCCCACCAACCCGCTGGAATATCCACGGACa                                                                                                                   >  4:81015/1‑274 (MQ=255)
                                                                                                                                                           ttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTca                                                                                                                                                                                                                                               <  6:59954/149‑1 (MQ=255)
                                                                                                                                                                        aTAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGgccg                                                                                                    <  2:81016/275‑1 (MQ=255)
                                                                                                                                                                        aTAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGcc                                                                                                                                                                                                                                                                     >  5:97603/1‑114 (MQ=255)
                                                                                                                                                                        aTAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGcc                                                                                                                                                                                                                                                                     <  6:97603/114‑1 (MQ=255)
                                                                                                                                                                           aCCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTTCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAAAAGCTCGGTTTCATGAAAAGATCACCGCTTGACCCAGTAACCGTCATACCAACGACACTAAACAATACCCCGGGCGTATACAGTGCAAATCATCCGAGGgcccc                                                                                                   >  4:449560/1‑271 (MQ=255)
                                                                                                                                                                            ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACActt                                                                                                                                                                                                                                          >  4:250437/1‑137 (MQ=255)
                                                                                                                                                                            ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                                                                                                                                                                                              <  8:376264/149‑1 (MQ=255)
                                                                                                                                                                            ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                                                                                                                                                                                              <  5:385435/149‑1 (MQ=255)
                                                                                                                                                                            ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                                                                                                                                                                                              >  8:281390/1‑149 (MQ=255)
                                                                                                                                                                             cTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                                                                                                                                                                                                  <  6:24701/112‑1 (MQ=255)
                                                                                                                                                                             cTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                                                                                                                                                                                                  >  5:24701/1‑112 (MQ=255)
                                                                                                                                                                             cTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                                                                                                                                                                                                 >  7:108827/1‑113 (MQ=255)
                                                                                                                                                                             cTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                                                                                                                                                                                                 <  8:108827/113‑1 (MQ=255)
                                                                                                                                                                                gCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACttcttc                                                                                                                                                                                                                                      <  2:250441/137‑1 (MQ=255)
                                                                                                                                                                                gCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCNCAGAACGACTCCCCTCCCCCGGAACAAGCCCGGTTTCATGcaa                                                                                                                                                                                         >  3:95284/1‑182 (MQ=255)
                                                                                                                                                                                   aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAACCGGCCTGATGTTAATGATAc                                                                                                                                                                                                                                                      <  8:137461/118‑1 (MQ=255)
                                                                                                                                                                                   aTAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAACCGGCCTGATGTTAATGATAc                                                                                                                                                                                                                                                      >  7:137461/1‑118 (MQ=255)
                                                                                                                                                                                    tAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt                                                                                                                                                                                     <  1:95284/182‑1 (MQ=255)
                                                                                                                                                                                     aGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAAc                                                                                                                                                      >  4:131075/1‑212 (MQ=255)
                                                                                                                                                                                       cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAct                                                                                                                                                                                                                          >  1:23212/1‑142 (MQ=255)
                                                                                                                                                                                         tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGaca                                                                                                                                                  <  2:131076/212‑1 (MQ=255)
                                                                                                                                                                                          aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTa                                                                                                                                                                                                                                                             >  3:23234/1‑104 (MQ=255)
                                                                                                                                                                                           ccGNAAAACANCGANAATTTTANAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                                                                                                                                                                                                      <  3:23212/142‑1 (MQ=255)
                                                                                                                                                                                              cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGa                                                                                                                                                                                                                                                         <  1:23234/104‑1 (MQ=255)
                                                                                                                                                                                              cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAg                                                                                                                               >  2:351255/1‑226 (MQ=255)
                                                                                                                                                                                                aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAc                                                                                                                                                                                                                                                      <  8:237928/105‑1 (MQ=255)
                                                                                                                                                                                                aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAc                                                                                                                                                                                                                                                      >  7:237928/1‑105 (MQ=255)
                                                                                                                                                                                                 aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATCACGCCGGAACAGCCAGTGCCAATCATCCGCGGGCCCCCCCCTTTCACCGTACCTTCCtacc                                                                           >  4:43860/1‑272 (MQ=255)
                                                                                                                                                                                                         tAATTTTATAATGGTGGACTTATCAGACGCTACGATCATTCATTGCAGGCAGTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATa                                                                   <  2:500033/275‑1 (MQ=255)
                                                                                                                                                                                                         gacagttttaATGTGGGCATGTTAAGACGCGACGATCATTCACTTCAGCCATCAGCCCACCCGGCCGGATGTTGTCCNGCCTGATGTTACTGATACCCTTTCCCACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGTTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATa                                                                   <  3:121590/267‑1 (MQ=255)
                                                                                                                                                                                                                      atggCATGATAAGACGCGACAATCAGTAATTGCAGGAATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACCCGTCTACCGCAGAACGACTCCCCTCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGc                                                     <  2:351035/274‑1 (MQ=255)
                                                                                                                                                                                                                       tggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTc                                                                                                                                                                                   >  7:61433/1‑149 (MQ=255)
                                                                                                                                                                                                                         gCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAAc                                                                       >  1:121590/1‑255 (MQ=255)
                                                                                                                                                                                                                                aaGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTCACCATAAAGCCCCTACTCGCCTCCttttt                                            >  2:326457/1‑275 (MQ=255)
                                                                                                                                                                                                                                         cTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACACTCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATACAGCTCATACCCGCCTCCTTTTTTCGAGTGaa                                   >  1:87452/1‑275 (MQ=255)
                                                                                                                                                                                                                                             cattcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                                                                                                                                             <  7:203673/149‑1 (MQ=255)
                                                                                                                                                                                                                                             cattcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCa                                                                                                                                                             <  5:302696/149‑1 (MQ=255)
                                                                                                                                                                                                                                              attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa                                                                                                                                                                                        >  6:243198/1‑121 (MQ=255)
                                                                                                                                                                                                                                              attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa                                                                                                                                                                                        <  5:243198/121‑1 (MQ=255)
                                                                                                                                                                                                                                                        ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGa                                                                                                                                                    >  1:93749/1‑147 (MQ=255)
                                                                                                                                                                                                                                                           aTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc                                                                                                                                               >  7:176873/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            ttAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACt                                                                                                                                                <  3:93749/147‑1 (MQ=255)
                                                                                                                                                                                                                                                              aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAAc                                                                                                                                            <  8:323614/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   acCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATAc                                                                                                                                       <  7:281390/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGgccg                                                                                                    <  1:60654/178‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACt                                                                          <  2:488422/204‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGa                                                                                                                                <  8:445020/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAAtttt  <  2:41766/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                               ttgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCTCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATc                                                                                                                           >  8:212433/2‑149 (MQ=255)
                                                                                                                                                                                                                                                                                 gATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCa                                                                                                                         <  5:317748/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                           
TCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTT  >  NC_000913/2027887‑2028427

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