Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,900,679 T→C 6.4% I303V (ATC→GTC)  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,6790TC 6.4% 226.7 / 5.4 78I303V (ATC→GTC) yoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base T (35/38);  new base C (3/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

TTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCA  >  NC_000913/1900533‑1900828
                                                                                                                                                  |                                                                                                                                                     
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                                                                                              <  1:1048837/150‑1 (MQ=255)
                                                      cAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttcttc                                                                                              <  1:1147732/150‑1 (MQ=255)
                                                      cAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGttcttc                                                                                              <  2:279876/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:786657/1‑150 (MQ=255)
                                                                                                   cccGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATc                                                 <  2:1265932/150‑1 (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:1346175/150‑1 (MQ=255)
                                                                                                                                     cgcgGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCg               <  2:1656651/150‑1 (MQ=255)
                                                                                                                                            tCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCa        >  1:1310149/1‑150 (MQ=255)
                                                                                                                                            tCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCa        <  2:230434/150‑1 (MQ=255)
                                                                                                                                               tGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAAc     <  1:1563686/150‑1 (MQ=255)
                                                                                                                                               tGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAAc     >  2:580104/1‑150 (MQ=255)
                                                                                                                                               tGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAAc     >  2:573094/1‑150 (MQ=255)
                                                                                                                                                ggtCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCATCTTCGGCAAATGCTCCTTCCGGGATGGGCATCGGCATCGGGTTATCGGCATCGTGCTGTACATTGGCCTGACGTTTCCCGCCCAGAAAACg    >  2:163410/3‑150 (MQ=255)
                                                                                                                                                gATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCAGCGGGGGATCGGCCTCGTGCTGAACATTGGCCTGCCGTTTCCCGCCCATCAAACg    >  1:1603889/1‑150 (MQ=255)
                                                                                                                                                  tCCCGCGCAGAGAACGCGACGCCAGCGTCCGTACGACGTGAATCATGTCTCGTTCTTCTTCGGCAAATGCACCTGCCGGGATCGGCATCGGCATCGGGTTATCGGCAGCGTGCTGAACAGTGGCCTGCCGTTTCCCGCCCAGCACACGCa  >  1:250564/1‑150 (MQ=255)
                                                                                                                                                  tCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGcc  >  2:707900/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTTAGCACGTACAATACCGATGATTTCATCCAGTTCACCGCGACATACCGGGAACAGACTGTGCGGTGAAGAGAGCAGTTGCTCGCGGATTTCATCGACCCCGAGATTAGCGTCAACCCAGCTTATTTCACCGCGCGGCGTCATGATCCCGCGCAGAGAACGCGACGCCAGCGTCAGTACGCCGTTAATCATGTAACGTTCTTCTTCGGCAAATGCACCTTCCGGGATCGGCATCGGCATCGGGTTATCGGCATCGTGCTGAACATTGGCCTGACGTTTCCCGCCCATCAAACGCA  >  NC_000913/1900533‑1900828

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: