Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,160 G→C 11.5% 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,1600GC11.5% 352.9 / 30.6 121intergenic (‑143/+28)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (50/56);  new base C (7/7);  total (58/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 122 total aligned reads displayed.

CTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAG  >  NC_000913/2028012‑2028309
                                                                                                                                                     |                                                                                                                                                     
ctgcACCCTGCCGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGAGAATTTTATAATGGTGGCATGATAAGACGCGACTATCATGCATTGCAGGCATTAGCACACCAGGCCGGATGTGa                                                                                                                                                      >  2:1515797/1‑150 (MQ=255)
ctgcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGa                                                                                                                                                      >  2:603857/1‑150 (MQ=255)
  gcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGCCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATc                                                                                                                                                    >  2:1613222/1‑150 (MQ=255)
    aCCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGTGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTCCTCCg                                                                                                                                                  >  1:1192572/1‑150 (MQ=255)
       cTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                                                                               <  2:544092/150‑1 (MQ=255)
        tGACGACTCCCACCCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                                                                              <  2:338990/150‑1 (MQ=255)
          ccgacTCCCAACCTGCTTTTTTATGAATCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                                                                            <  2:1571630/149‑1 (MQ=255)
          acgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCTTACCGCAAAACAGCGATAATTTTATAATGGTGGGATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                                                                            <  1:953531/150‑2 (MQ=255)
          acgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                                                                            <  1:1579207/150‑1 (MQ=255)
          acgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                                                                            >  1:1453365/1‑150 (MQ=255)
          acgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                                                                            <  2:459248/150‑2 (MQ=255)
           cgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                                                                           <  1:535298/150‑1 (MQ=255)
           cgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGGACCGAAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                                                                           <  1:911806/150‑1 (MQ=255)
            gacTCACAACCTGCTTTTTTATGATTCGGGCATAACCTGGCGATAGCGTACCGCAATACAGCGATAATTTTATAATGGTGGCATGAAAAGACGCGACTATCAATCATTGCACGCATTAGCACACCAGGCCGGATCACATCCGGCCTGGTg                                                                                                                                          <  1:311015/150‑1 (MQ=255)
             acTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt                                                                                                                                         <  2:1464360/150‑1 (MQ=255)
              cTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt                                                                                                                                        <  1:1103153/150‑1 (MQ=255)
              atCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGTAAAACAGCGATAATTTTATAATGGTGGTATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                                                                        <  1:401269/149‑2 (MQ=255)
                 ccAACCTGCTTTTTTTATGATTCTGGCATACCCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtgct                                                                                                                                      <  1:937405/150‑4 (MQ=255)
                 ccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                                                                     <  2:1314285/150‑1 (MQ=255)
                  cAACCTTCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                                                                    <  1:868977/150‑1 (MQ=255)
                         cTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCt                                                                                                                             <  2:605810/150‑1 (MQ=255)
                           ttttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCtt                                                                                                                            >  1:1574002/1‑150 (MQ=255)
                            tttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAACCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCttt                                                                                                                           <  1:1584883/150‑1 (MQ=255)
                            tttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAACGGTTGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCttt                                                                                                                           <  1:918256/150‑1 (MQ=255)
                             ttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                                                                          <  2:1522518/150‑1 (MQ=255)
                             ttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                                                                          <  1:892348/150‑1 (MQ=255)
                                tATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCACGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcaca                                                                                                                       >  2:1086000/1‑150 (MQ=255)
                                  tGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct                                                                                                                     >  1:1300661/1‑150 (MQ=255)
                                  tGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct                                                                                                                     >  1:1280658/1‑150 (MQ=255)
                                        tGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATATTTTTATAATGGTGGCATGATAAGACGCGTCTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGTTCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTcc                                                                                                               <  1:300273/150‑1 (MQ=255)
                                         ggCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCg                                                                                                              <  1:466915/150‑1 (MQ=255)
                                         ggCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCg                                                                                                              <  1:1308858/150‑1 (MQ=255)
                                          gCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc                                                                                                             <  1:828627/150‑1 (MQ=255)
                                           cATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCa                                                                                                            >  2:1290656/1‑150 (MQ=255)
                                               aCCTGGCGATATCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                                                                        <  2:1052360/150‑1 (MQ=255)
                                                ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATATGACGCGCCTATCATTCATTGCAGGCATTAGCACCCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTACACACTTCTTCCGCAGAACg                                                                                                       >  1:597724/1‑150 (MQ=255)
                                                  tGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAc                                                                                                     <  2:813777/150‑1 (MQ=255)
                                                     cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                                                                                  <  2:1212437/150‑1 (MQ=255)
                                                     cGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTCATGTTAATGATACCCTTTCACACGTCTTCCGCAGAACGActcc                                                                                                  >  1:1437576/1‑150 (MQ=255)
                                                             tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTccccc                                                                                          <  2:1077873/150‑1 (MQ=255)
                                                               ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCgg                                                                                        >  1:183150/1‑150 (MQ=255)
                                                               ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCgg                                                                                        >  1:1370591/1‑150 (MQ=255)
                                                                cGCAAAAAAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGa                                                                                       <  2:1553560/150‑1 (MQ=255)
                                                                 gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACGCTTTCACACTTCTTCCGCAGAACGACTCGCCTCCCCCGGaa                                                                                      >  1:1577717/1‑150 (MQ=255)
                                                                 gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGAGATCCGGCATGATGTTAATGATACCCTTTCACACTTCTTCCGCAGACCGACTCCACTCCCCCGGaa                                                                                      >  1:880138/1‑150 (MQ=255)
                                                                  cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGACGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCTCCTCCCCCGGAAc                                                                                     <  2:1157039/150‑1 (MQ=255)
                                                                     aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGCCTCCCCTCCCCCGGAACaac                                                                                  >  2:901217/1‑149 (MQ=255)
                                                                     aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAg                                                                                  <  1:1102049/150‑1 (MQ=255)
                                                                     aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAg                                                                                  >  1:28312/1‑150 (MQ=255)
                                                                       cAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                                                                                >  1:1728797/1‑150 (MQ=255)
                                                                       cAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAAGATACACTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGTAACAAGca                                                                                >  2:152322/1‑149 (MQ=255)
                                                                         gCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg                                                                              <  2:757930/150‑1 (MQ=255)
                                                                           gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGt                                                                            <  2:1292405/150‑1 (MQ=255)
                                                                           gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGt                                                                            >  1:1681763/1‑150 (MQ=255)
                                                                             aaaTTTTATAATGGGGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGttt                                                                          <  1:1461600/149‑1 (MQ=255)
                                                                                ttttATAATGGGGGCAGGTTAAGACTCGACTATCATTCATTGCAGGATTTAGCACACCAGGCCGGTTGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGATCAAGCTCGGTTTCAt                                                                       <  1:152322/150‑1 (MQ=255)
                                                                                  ttAAAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGa                                                                     <  1:879557/150‑1 (MQ=255)
                                                                                   tataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGACCGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                                                                    >  1:1182069/1‑150 (MQ=255)
                                                                                    ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTATCACACCAGGCCGGATGTGAACCGGCCTGATGTTAATGATACCCTATCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGTATACATGaaa                                                                   >  2:339021/1‑150 (MQ=255)
                                                                                          gtggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTc                                                             <  1:641949/150‑1 (MQ=255)
                                                                                                 gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg                                                      <  2:1728797/150‑1 (MQ=255)
                                                                                                 gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg                                                      <  2:651671/150‑1 (MQ=255)
                                                                                                 gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg                                                      <  2:1384289/150‑1 (MQ=255)
                                                                                                  aTAAGACTCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCGCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa                                                     <  1:1021921/150‑1 (MQ=255)
                                                                                                     aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAccc                                                  <  1:1444911/150‑1 (MQ=255)
                                                                                                     aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAccc                                                  >  1:723809/1‑150 (MQ=255)
                                                                                                     aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAccc                                                  >  2:1648600/1‑150 (MQ=255)
                                                                                                     aGACGCGACTATCAGTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGcccc                                                  >  1:1354180/1‑150 (MQ=255)
                                                                                                      gACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcccc                                                 >  2:1655812/1‑150 (MQ=255)
                                                                                                        cgcgAATATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGt                                               <  2:199343/150‑1 (MQ=255)
                                                                                                          cgACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTaa                                             >  1:1101739/1‑150 (MQ=255)
                                                                                                            aCTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGGGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAAcc                                           <  2:872221/150‑1 (MQ=255)
                                                                                                             cTATCATTCCTTGGCGGCATTTGCACCCCAGGCCGGTTGTGATGCGGCCTGATGTTAATGATACCCTTTCACACATCTTCCGGAGAACGACTCCCCTCCCCCGGAACAAGCTCTGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCg                                          <  1:1083795/150‑1 (MQ=255)
                                                                                                                  attcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATa                                     >  1:1558360/1‑150 (MQ=255)
                                                                                                                     cattGCAGGCATTAGCACACCAGGCCGGATTTGATCCGGCCTTATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATTAACAGATTCCCGCTTGACCCCGTACCCGTCATAccc                                  >  1:307328/1‑149 (MQ=255)
                                                                                                                     cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCa                                  <  2:28312/150‑1 (MQ=255)
                                                                                                                         gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGa                              >  2:446231/1‑150 (MQ=255)
                                                                                                                          cAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGac                             >  1:1708161/1‑150 (MQ=255)
                                                                                                                          cAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGac                             >  1:595970/1‑150 (MQ=255)
                                                                                                                           aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATCCCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATAAAAAGCTTCCCGCTTGCCCCCGTAACCGTCATACCCACGacc                            >  1:73619/1‑149 (MQ=255)
                                                                                                                            ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGacac                           >  2:1588181/1‑150 (MQ=255)
                                                                                                                            ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGacac                           <  1:943231/150‑1 (MQ=255)
                                                                                                                            ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGacac                           <  1:645646/150‑1 (MQ=255)
                                                                                                                            ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAGGTTTATGATACCCTTTCACACTTCTTCCGCAGACCGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAATGCTGCCCGCTTGAGCGAGTAACCGTCATCCCAGCGacac                           >  2:1430890/1‑150 (MQ=255)
                                                                                                                              cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc                         >  2:1433260/1‑150 (MQ=255)
                                                                                                                              cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc                         <  1:465255/150‑1 (MQ=255)
                                                                                                                              cATTAGCACACCAGGCCGGATGTGATCCGGCCGGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTGCATGAAAAGCTTCCCGCTTGACCCCGGAACCGTCATACCAACGAGACta                         >  1:1504828/1‑149 (MQ=255)
                                                                                                                                  ttaacatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGCCTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGGTTCCCGCTTGACCCCGTACCCGTCCTCCCAACGCCACTCGAcc                     >  1:115024/4‑150 (MQ=255)
                                                                                                                                  aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAAcc                     <  2:657002/150‑1 (MQ=255)
                                                                                                                                    cacaCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCAt                   >  1:814633/1‑150 (MQ=255)
                                                                                                                                     acatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATa                  >  1:1415184/1‑150 (MQ=255)
                                                                                                                                      catcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATAc                 >  2:1312490/4‑150 (MQ=255)
                                                                                                                                      cacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATAc                 >  1:1499407/1‑150 (MQ=255)
                                                                                                                                         cAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGc              >  2:1442250/1‑150 (MQ=255)
                                                                                                                                          aGGCCGGATCACCTCCGGCCTGATGTTAATGATACCCATTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGACCAAGCTCGGTTTCCTGAAAAGCTTCCCGCTAGCCCCCGACCCCGTCATACCAACGCCACTCACGCAGACCGcg             >  1:327549/1‑149 (MQ=255)
                                                                                                                                          aGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGCCTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCt             >  1:1589656/1‑150 (MQ=255)
                                                                                                                                            gCCGGATGTGAGCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTgg           <  1:1246503/150‑1 (MQ=255)
                                                                                                                                              cgtaTGTGATGCGGACTGATGTTAATGAAACCCTTTCTCACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCTGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCATCCATACCGCTGGAg         <  1:374428/147‑1 (MQ=255)
                                                                                                                                               ggATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCACGCTTGACCCCGTAACCGTCATACCAACGACAGTCAACCATACCGCTGGAGt        >  2:282218/1‑150 (MQ=255)
                                                                                                                                                     gATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAg  <  2:183150/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAG  >  NC_000913/2028012‑2028309

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