Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,900,672 G→A 6.5% T305M (ACG→ATG) ‡ yoaE ← putative membrane protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,900,6720GA 6.5% 232.4 / 4.0 77T305M (ACG→ATG) ‡yoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base G (32/40);  new base A (3/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCG  >  NC_000913/1900528‑1900815
                                                                                                                                                |                                                                                                                                               
cgTTCTTTAGCACGTACACTACCGATGATTTCATCCAGTTCACGGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGGTGCTCGCGGATTTCATCGACCCCGTGATTAGCGGCTACCCAGCTTATTTCACCGCGCGGCGTCATg                                                                                                                                            <  1:127573/149‑1 (MQ=255)
agTTCTTTAGCACGTACATTACCGATGATTTCATCCAGTTCACCGCTACATACCGGGAAGAGACTGTGCGGTGAAGAGAGCAGTTGCTCGAGGATTTCTTCGTCCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATg                                                                                                                                            <  2:751995/150‑1 (MQ=255)
agTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCATTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATg                                                                                                                                            <  2:1131362/150‑1 (MQ=255)
 gTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGGTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGa                                                                                                                                           <  2:1152299/150‑1 (MQ=255)
 gTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGa                                                                                                                                           >  2:40277/1‑150 (MQ=255)
     tttAGCACGTACAATACCGATGAGTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCATCTTATTTCACCGCGCGGCGTCATGATccc                                                                                                                                       <  1:730684/150‑1 (MQ=255)
         gCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGAGTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCgcgc                                                                                                                                   <  1:1100927/150‑1 (MQ=255)
           aCGTACAATACCGATGATATCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTAGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCTCGCGGGATCATGACGCCgcgcgg                                                                                                                                 <  1:753062/150‑3 (MQ=255)
             gTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCagag                                                                                                                               <  2:1336123/150‑1 (MQ=255)
              tACAATACCGATGATTTCATCCAGTTCACCTCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCCAGATTAGAGTAAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCTCGCAgaga                                                                                                                              >  1:1055973/1‑150 (MQ=255)
                cAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAAc                                                                                                                            <  2:371363/150‑1 (MQ=255)
                   tACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAAcgcg                                                                                                                         >  2:1437481/1‑150 (MQ=255)
                     ccGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGAc                                                                                                                       <  1:1726154/150‑1 (MQ=255)
                     ccGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCAGGATCCCGCGCAGAGAACGCGAc                                                                                                                       >  2:1276727/1‑150 (MQ=255)
                      cGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACg                                                                                                                      >  1:1265932/1‑150 (MQ=255)
                       gatgatTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGc                                                                                                                     <  2:1371367/150‑1 (MQ=255)
                        atgatTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGcc                                                                                                                    <  2:581925/150‑1 (MQ=255)
                          gatTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGcccg                                                                                                                  >  1:1603000/1‑148 (MQ=255)
                          gatTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAg                                                                                                                  <  2:776883/150‑1 (MQ=255)
                             ttCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGCCGCCAGCGt                                                                                                               >  1:230434/1‑150 (MQ=255)
                               cATCCAGTTCACCGCGACATACCGGGAACAGCCTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGGATCATGACGCCGCGCGGTGaaa                                                                                                                            <  1:1222555/135‑2 (MQ=255)
                               cATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGGATCATGACGCCGCGCGGTGaaa                                                                                                                            >  2:1222555/1‑134 (MQ=255)
                                aTCCAGTTCACCGCGACATACCGGGAACAGACTGTTCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGGGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAg                                                                                                            <  2:1013480/150‑1 (MQ=255)
                                aTCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTCGCGTCACCCCAGCTTATTTCACAGCGCGGCGTCATGCTCCCGCGCAGAGCACGCGACGCCAGCGTCAg                                                                                                            >  1:1171732/1‑150 (MQ=255)
                                   cAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTATCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGAAGAGAACGCGACGCCAGCGTCAGTAc                                                                                                         >  2:822741/1‑150 (MQ=255)
                                          ccGCGACATACCGGGATCAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGGCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTa                                                                                                  >  2:1120234/1‑150 (MQ=255)
                                              gACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCa                                                                                              <  1:388394/150‑1 (MQ=255)
                                              gACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACACCAGCGTCAGTACGCCGTTAATCa                                                                                              >  1:1596678/1‑150 (MQ=255)
                                                    ccGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAAc                                                                                        >  1:389804/1‑150 (MQ=255)
                                                        gAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttc                                                                                    >  1:1191598/1‑150 (MQ=255)
                                                        gAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGAGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGCCGCCAGCGTCAGTACGCCGTTAATCATGTAACGgtc                                                                                    >  2:1395768/1‑150 (MQ=255)
                                                         aaCAGACTGTTCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttct                                                                                   <  2:415874/150‑1 (MQ=255)
                                                           cAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttcttc                                                                                 <  2:279876/150‑1 (MQ=255)
                                                           cAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttcttc                                                                                 <  1:1147732/150‑1 (MQ=255)
                                                           cAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttcttc                                                                                 <  1:1048837/150‑1 (MQ=255)
                                                             gACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttcttctt                                                                               >  1:1656651/1‑150 (MQ=255)
                                                               cTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTCCTTCTTCg                                                                             >  1:1663166/1‑150 (MQ=255)
                                                                 gtgCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGc                                                                           >  1:1399155/1‑150 (MQ=255)
                                                                     ggTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAAt                                                                       >  1:1733780/1‑150 (MQ=255)
                                                                        gAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCa                                                                    >  1:980523/1‑150 (MQ=255)
                                                                          agagagCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCAcc                                                                  <  1:729664/150‑1 (MQ=255)
                                                                                       tCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggc                                                     <  1:1578913/150‑1 (MQ=255)
                                                                                        cgcgGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggca                                                    <  2:1214609/150‑1 (MQ=255)
                                                                                         gcgGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGACATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTACATCTTCGGCAAAAGCACCTTCCGGGatcggcat                                                   >  1:202161/1‑150 (MQ=255)
                                                                                          cgGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggcatc                                                  >  1:736755/1‑150 (MQ=255)
                                                                                            gATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggcatcgg                                                >  2:207899/1‑150 (MQ=255)
                                                                                            gATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGatcggcatcgg                                                >  1:1694463/1‑150 (MQ=255)
                                                                                                 cATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGgcatcggcatcg                                           >  2:51673/1‑150 (MQ=255)
                                                                                                   tCGACGCCGAGATTAGCGTCAACCCAGCTTATTTCATCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCggg                                         <  1:1676291/150‑1 (MQ=255)
                                                                                                   tCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCggg                                         <  1:844197/150‑1 (MQ=255)
                                                                                                   tCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCggg                                         <  1:396530/150‑1 (MQ=255)
                                                                                                    cGACCCCGAGATTAGCGTCCACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGt                                        >  2:1294522/1‑150 (MQ=255)
                                                                                                    cGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGt                                        <  2:1033594/150‑1 (MQ=255)
                                                                                                       ccccGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTAt                                     >  2:633640/1‑150 (MQ=255)
                                                                                                        cccGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATc                                    <  2:1265932/150‑1 (MQ=255)
                                                                                                        cccGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATc                                    >  2:786657/1‑150 (MQ=255)
                                                                                                         ccGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCg                                   <  2:980523/150‑1 (MQ=255)
                                                                                                          cGAGATTAGCGTCAACCCAGCTTATTTCACCGGGCTGCGTCATGATCCCGGGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCgg                                  <  1:673365/150‑1 (MQ=255)
                                                                                                           gagaTTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGc                                 <  2:1198265/150‑1 (MQ=255)
                                                                                                            agaTTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCa                                >  2:588764/1‑150 (MQ=255)
                                                                                                                  gCGTCAAGCCAGCTTATTTCACCGCGCGGCGGAATGATCCCGCGCAGAGAACGCGTCGCCAGCTTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGc                          <  2:1603000/150‑1 (MQ=255)
                                                                                                                   cGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCt                         >  1:1346175/1‑150 (MQ=255)
                                                                                                                    ggcAACGCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCTTCAGTACGCCGTTAATCATGTAACGTTCTTCATCGGCAAATGCACCTTGCTGGATCGGCATCGGCATCGGGTAATCGGCATCGTGCTg                        <  2:312027/148‑1 (MQ=255)
                                                                                                                    gTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGAAGAGAACGAGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTg                        <  2:1072190/150‑1 (MQ=255)
                                                                                                                      cAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCAACTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGaa                      >  2:869487/1‑150 (MQ=255)
                                                                                                                        aCCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTACGGCAAATGCACCTTCCGGGAGCGGCATCGGCAGCGGGGTATCGGCTTCGTGCTGAACa                    >  2:279484/1‑150 (MQ=255)
                                                                                                                         cccAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACAt                   >  2:640236/1‑150 (MQ=255)
                                                                                                                            aGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTgg                <  1:1314669/150‑1 (MQ=255)
                                                                                                                              cTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGcc              <  2:389804/150‑1 (MQ=255)
                                                                                                                               ttATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCt             <  2:1360988/150‑1 (MQ=255)
                                                                                                                                   ttCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGGTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACg         <  1:855730/150‑1 (MQ=255)
                                                                                                                                      tctgcgcgGGATCATGACGCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGttt      >  1:1701444/4‑150 (MQ=255)
                                                                                                                                       ctgcgcgGGATCATGACGCCGCGCAGAGAACGCGACGCAAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCGTCGTCAACTGCACCTTCCGGGATCGGCATCGGCATCGGGTGATCGGCAGCGTGCTGAACATTGGCCGGCCGTGTc     >  1:847816/3‑150 (MQ=255)
                                                                                                                                       ccGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTc     <  1:771253/150‑1 (MQ=255)
                                                                                                                                        cgcgcgGCGTCATGATCCCGCGCAGAGAACGCGACTCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCCTCGTGCTGAACATTGGCCTGACGTTTcc    <  2:1504105/150‑1 (MQ=255)
                                                                                                                                          cgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCg  <  2:1656651/150‑1 (MQ=255)
                                                                                                                                          cgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCg  <  2:1346175/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AGTTCTTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCG  >  NC_000913/1900528‑1900815

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