Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,158 G→C 10.7% intergenic (‑141/+30) 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,1580GC10.7% 361.8 / 29.4 121intergenic (‑141/+30)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (51/56);  new base C (6/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 = 9.23e-01

Only 100 of 122 total aligned reads displayed.

GTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGT  >  NC_000913/2028009‑2028303
                                                                                                                                                      |                                                                                                                                                 
gTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtg                                                                                                                                                   >  1:872221/1‑150 (MQ=255)
   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)
             acgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCTTACCGCAAAACAGCGATAATTTTATAATGGTGGGATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                                                                      <  1:953531/150‑2 (MQ=255)
             acgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                                                                      >  1:1453365/1‑150 (MQ=255)
             acgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                                                                      <  1:1579207/150‑1 (MQ=255)
             acgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                                                                      <  2:1408279/150‑2 (MQ=255)
             acgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                                                                      <  2:459248/150‑2 (MQ=255)
              cgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                                                                     <  1:535298/150‑1 (MQ=255)
               gacTCACAACCTGCTTTTTTATGATTCGGGCATAACCTGGCGATAGCGTACCGCAATACAGCGATAATTTTATAATGGTGGCATGAAAAGACGCGACTATCAATCATTGCACGCATTAGCACACCAGGCCGGATCACATCCGGCCTGGTg                                                                                                                                    <  1:311015/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)
                     cAACCTTCTTTTTTATTATTCTGGCATAACCTGGCGATAGCGTAACGCAATACAGCTATTTTTTTATATTGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCGGGCCTGATGTTAATg                                                                                                                              <  2:1264819/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)
                                                   ccTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATATGACGCGCCTATCATTCATTGCAGGCATTAGCACCCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTACACACTTCTTCCGCAGAACg                                                                                                 >  1:597724/1‑150 (MQ=255)
                                                     tGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAc                                                                                               <  2:813777/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)
                                                                                     ttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGa                                                               <  1:912376/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)
                                                                                                     aTAAGACTCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCGCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGa                                               <  1:1021921/150‑1 (MQ=255)
                                                                                                        aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAccc                                            >  1:723809/1‑150 (MQ=255)
                                                                                                        aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAccc                                            >  2:1648600/1‑150 (MQ=255)
                                                                                                        aGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAccc                                            <  1:1444911/150‑1 (MQ=255)
                                                                                                        aGACGCGACTATCAGTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGcccc                                            >  1:1354180/1‑150 (MQ=255)
                                                                                                          aCGTGACTATCATTCTTTGCAGGCATTAGCACACCAGGCCGGATGTGCTCCGGCCTGATGTTAAAGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCACCGGAACATGCTAGGTTTCATGAATAGCTTCCCGCTTGACCCCg                                          <  2:1053596/150‑1 (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                     <  1:645646/150‑1 (MQ=255)
                                                                                                                               ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGacac                     <  1:943231/150‑1 (MQ=255)
                                                                                                                               ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGacac                     >  2:1588181/1‑150 (MQ=255)
                                                                                                                               ggCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAGGTTTATGATACCCTTTCACACTTCTTCCGCAGACCGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAATGCTGCCCGCTTGAGCGAGTAACCGTCATCCCAGCGacac                     >  2:1430890/1‑150 (MQ=255)
                                                                                                                                 cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc                   <  1:465255/150‑1 (MQ=255)
                                                                                                                                 cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTc                   >  2:1433260/1‑150 (MQ=255)
                                                                                                                                 cATTAGCACACCAGGCCGGATGTGATCCGGCCGGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTGCATGAAAAGCTTCCCGCTTGACCCCGGAACCGTCATACCAACGAGACta                   >  1:1504828/1‑149 (MQ=255)
                                                                                                                                     ttaacatCAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGCCTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGGTTCCCGCTTGACCCCGTACCCGTCCTCCCAACGCCACTCGAcc               >  1:115024/4‑150 (MQ=255)
                                                                                                                                     aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAAcc               <  2:1574002/150‑1 (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)
                                                                                                                                               gCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACTAGCTCGGTTTCATGAAAAGCTTCCCGCTTGCCCCCGTAACCGACTTACCAACGACACTCAACCATACCGCTgg     >  2:150361/1‑150 (MQ=255)
                                                                                                                                               gCCGGATGTGAGCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTgg     <  1:1246503/150‑1 (MQ=255)
                                                                                                                                                 cgtaTGTGATGCGGACTGATGTTAATGAAACCCTTTCTCACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCTGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCATCCATACCGCTGGAg   <  1:374428/147‑1 (MQ=255)
                                                                                                                                                 cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGAATCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTACCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAg   <  2:1694038/150‑1 (MQ=255)
                                                                                                                                                  ggATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCACGCTTGACCCCGTAACCGTCATACCAACGACAGTCAACCATACCGCTGGAGt  >  2:282218/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGT  >  NC_000913/2028009‑2028303

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: