Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,199,607 A→G 6.8% intergenic (+2309/‑672) yehF → / → yehI putative protein YehF/DUF4132 domain‑containing protein YehI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,199,6070AG 6.8% 395.7 / 13.7 133intergenic (+2309/‑672)yehF/yehIputative protein YehF/DUF4132 domain‑containing protein YehI
Reads supporting (aligned to +/- strand):  ref base A (51/73);  new base G (4/5);  total (55/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 139 total aligned reads displayed.

CCAAACCGAACAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCTTTTC  >  NC_000913/2199458‑2199756
                                                                                                                                                     |                                                                                                                                                     
ccAAACCGAACAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGaaa                                                                                                                                                       <  1:1510457/150‑1 (MQ=255)
  aaaCCGAACAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAAc                                                                                                                                                     <  1:1929557/150‑1 (MQ=255)
  aaaCCGAACAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAAc                                                                                                                                                     <  2:1540347/150‑1 (MQ=255)
      cGAACAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTcc                                                                                                                                                 <  1:1263769/150‑1 (MQ=255)
         aCAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGc                                                                                                                                              <  2:1267384/150‑1 (MQ=255)
          cAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGct                                                                                                                                             >  1:1732344/1‑150 (MQ=255)
           aaaaGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAGTTTTCCAGcga                                                                                                                                            <  2:337511/150‑3 (MQ=255)
            aaaGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAGTTTTCCAGcgac                                                                                                                                           <  1:2201363/150‑4 (MQ=255)
            aaaGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGctct                                                                                                                                           >  1:71295/1‑150 (MQ=255)
            aaaGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGctct                                                                                                                                           >  1:1148962/1‑150 (MQ=255)
            aaaGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGctct                                                                                                                                           >  1:1249458/1‑150 (MQ=255)
             aaGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCtt                                                                                                                                          <  2:1787645/150‑1 (MQ=255)
                cGATGATGAAACTCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCt                                                                                                                                       <  2:1857803/150‑1 (MQ=255)
                  atgatgAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTgg                                                                                                                                     >  2:62206/1‑150 (MQ=255)
                   tgatgaAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGt                                                                                                                                    <  1:2341012/150‑1 (MQ=255)
                   tgatgaAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGt                                                                                                                                    >  1:9808/1‑150 (MQ=255)
                    gatgaAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGtt                                                                                                                                   <  1:1837427/150‑1 (MQ=255)
                     atgaAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTg                                                                                                                                  <  2:837027/150‑1 (MQ=255)
                       gaAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAg                                                                                                                                <  2:1178813/150‑1 (MQ=255)
                        aaaCGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGc                                                                                                                               >  1:796649/1‑150 (MQ=255)
                        aaaCGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGc                                                                                                                               >  1:1794507/1‑150 (MQ=255)
                           cgcgGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATc                                                                                                                            >  2:783587/1‑150 (MQ=255)
                                 cGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTgg                                                                                                                      <  2:190902/150‑1 (MQ=255)
                                 cGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTgg                                                                                                                      <  2:2274030/150‑1 (MQ=255)
                                 cGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTgg                                                                                                                      >  1:124304/1‑150 (MQ=255)
                                  gAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGt                                                                                                                     >  1:2130382/1‑150 (MQ=255)
                                       aTGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTc                                                                                                                <  2:1999120/150‑1 (MQ=255)
                                       aTGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTc                                                                                                                <  1:512916/150‑1 (MQ=255)
                                       aTGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTc                                                                                                                <  1:336146/150‑1 (MQ=255)
                                         gCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCAc                                                                                                              <  1:1830365/150‑1 (MQ=255)
                                         gCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCAc                                                                                                              <  1:575969/150‑1 (MQ=255)
                                           ggTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACtt                                                                                                            >  2:1688659/1‑150 (MQ=255)
                                           ggTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCCGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACtt                                                                                                            <  2:444787/150‑1 (MQ=255)
                                                 ccGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACccg                                                                                                      <  1:1850375/150‑1 (MQ=255)
                                                   gCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAAcc                                                                                                        >  1:1105799/1‑146 (MQ=255)
                                                      aCAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTc                                                                                                 >  2:2451678/1‑150 (MQ=255)
                                                      aCAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTc                                                                                                 <  2:2130382/150‑1 (MQ=255)
                                                       cAAATTGCTGAAAAAAGAGGCGCGTACCATCGCCGCCCAGCAGGCAGCAAGGCGGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCt                                                                                                <  2:1790485/150‑1 (MQ=255)
                                                         aaTTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGa                                                                                              >  1:1169177/1‑150 (MQ=255)
                                                              cTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTgg                                                                                         <  1:1483220/150‑1 (MQ=255)
                                                               tGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTggg                                                                                        <  1:1175483/150‑1 (MQ=255)
                                                               tGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTggg                                                                                        >  1:1373166/1‑150 (MQ=255)
                                                                  aaaaaGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCg                                                                                     >  2:483069/1‑150 (MQ=255)
                                                                       gATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTtat                                                                                >  1:2269179/1‑150 (MQ=255)
                                                                        aTGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTtata                                                                               <  1:2405012/150‑1 (MQ=255)
                                                                        aTGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTtata                                                                               <  2:358865/150‑1 (MQ=255)
                                                                        aTGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTtata                                                                               <  2:1291090/150‑1 (MQ=255)
                                                                        aTGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTtata                                                                               <  1:956787/150‑1 (MQ=255)
                                                                           cgcgTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAgcg                                                                            >  1:1758466/1‑150 (MQ=255)
                                                                           cgcgTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAAGCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAgcg                                                                            <  1:2153311/150‑1 (MQ=255)
                                                                            gcgTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAgcgc                                                                           <  1:1197367/150‑1 (MQ=255)
                                                                            gcgTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAgcgc                                                                           <  1:2039673/150‑1 (MQ=255)
                                                                               tACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGTAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGa                                                                        <  1:2483683/150‑1 (MQ=255)
                                                                               tACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGa                                                                        >  1:1048412/1‑150 (MQ=255)
                                                                               tACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGa                                                                        >  1:2485045/1‑150 (MQ=255)
                                                                               tACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGa                                                                        >  1:2038886/1‑150 (MQ=255)
                                                                                  cATCGCCGCCAAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCAGCCGTCTGATTTGGGGCGTTTATAGCGCCGaaaa                                                                     >  1:541103/1‑150 (MQ=255)
                                                                                    tCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAAcc                                                                   >  1:1643095/1‑150 (MQ=255)
                                                                                     cgccgcCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCa                                                                  <  2:1148962/150‑1 (MQ=255)
                                                                                     cgccgcCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCa                                                                  <  2:1014536/150‑1 (MQ=255)
                                                                                     cgccgcCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCa                                                                  <  1:748836/150‑1 (MQ=255)
                                                                                       ccgccCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGc                                                                <  1:675255/150‑1 (MQ=255)
                                                                                       ccgccCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGc                                                                <  2:711030/150‑1 (MQ=255)
                                                                                         gccCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTa                                                              >  1:870086/1‑150 (MQ=255)
                                                                                          cccAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTAc                                                             >  2:2092375/1‑150 (MQ=255)
                                                                                             agcagGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTgg                                                          >  2:2479650/1‑150 (MQ=255)
                                                                                                 ggTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGgcttg                                                      >  2:2286185/1‑150 (MQ=255)
                                                                                                 ggTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGgcttg                                                      <  2:2167006/150‑1 (MQ=255)
                                                                                                 ggTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGgcttg                                                      <  2:44781/150‑1 (MQ=255)
                                                                                                   tAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGgcttgct                                                    <  2:836043/150‑1 (MQ=255)
                                                                                                       aaGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTcg                                                <  1:1688659/150‑1 (MQ=255)
                                                                                                       aaGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTcg                                                <  2:71295/150‑1 (MQ=255)
                                                                                                       aaGGCTGGAATCCGCCATGTGCCTGCGCCGCCACTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTcg                                                <  2:2098426/150‑1 (MQ=255)
                                                                                                         ggCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTcgcg                                              <  2:2361038/150‑1 (MQ=255)
                                                                                                         ggCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTcgcg                                              <  2:796649/150‑1 (MQ=255)
                                                                                                            tGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAg                                           <  2:1859124/150‑1 (MQ=255)
                                                                                                              gAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCa                                         <  1:2223581/150‑1 (MQ=255)
                                                                                                                 accGGCATGTTCCTGCGGCGGCGGTGGGTGGTGGAAAACTTCCATGTCTTCCTTGTTGTTCAATCGGTTGTTGCGCACTTAAACCGGCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCagaa                                      <  2:1668781/149‑1 (MQ=255)
                                                                                                                       aTGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAAc                                <  1:119053/150‑1 (MQ=255)
                                                                                                                       aTGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAAc                                <  2:2272670/150‑1 (MQ=255)
                                                                                                                        tgtgCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGATGATAACa                               <  1:827332/150‑1 (MQ=255)
                                                                                                                        tgtgCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACa                               >  2:1841586/1‑150 (MQ=255)
                                                                                                                           gCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCt                            >  1:1544273/1‑150 (MQ=255)
                                                                                                                            ccTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTc                           >  2:322942/1‑150 (MQ=255)
                                                                                                                                   cgccgcTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCAcc                    >  2:797411/1‑150 (MQ=255)
                                                                                                                                    gccgcTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCg                   >  1:2248550/1‑150 (MQ=255)
                                                                                                                                      cgcTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCt                 >  1:886375/1‑150 (MQ=255)
                                                                                                                                         gaagagCTGGAAGTTTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGAc              >  2:163814/6‑150 (MQ=255)
                                                                                                                                         gaagagCTGGAAGTTTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGAc              >  2:273717/6‑150 (MQ=255)
                                                                                                                                          aagagCTGGAAGTTTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACg             >  1:140145/5‑150 (MQ=255)
                                                                                                                                            tCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACgat           <  2:1169177/150‑1 (MQ=255)
                                                                                                                                              gCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACcgatgacgatga         >  2:2073146/1‑150 (MQ=255)
                                                                                                                                               cTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACgatgat        <  1:641408/150‑1 (MQ=255)
                                                                                                                                                tGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCTCGTAGCAGAAGATAACAGCTCGATCACCGCTGACGATGATc       <  2:541103/150‑1 (MQ=255)
                                                                                                                                                tGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATc       <  1:62206/150‑1 (MQ=255)
                                                                                                                                                tGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATc       <  2:1249458/150‑1 (MQ=255)
                                                                                                                                                 ggAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCt      <  1:1990671/150‑1 (MQ=255)
                                                                                                                                                   aaaaCTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCttt    <  1:1132141/150‑1 (MQ=255)
                                                                                                                                                     aaCTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCTTTTc  <  2:1373166/150‑1 (MQ=255)
                                                                                                                                                     aaCTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCTTTTc  >  2:2439541/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAAACCGAACAAAAGCGATGATGAAACGCGGGCGAACGATGCGGTTAACCGCTACAAATTGCTGAAAAAAGATGCGCGTACCATCGCCGCCCAGCAGGTAGCAAGGCTGGAATCCGCCATGTGCCTGCGCCGCCGCTGGTCGCTGGAAAACTTCCAGCTCTTCCTGGTTGAGCATCCGCTGGTTCGTCACTTAACCCGCCGTCTGATTTGGGGCGTTTATAGCGCCGAAAACCAGCTACTGGCTTGCTTTCGCGTAGCAGAAGATAACAGCTCCAGCACCGCTGACGATGATCTTTTC  >  NC_000913/2199458‑2199756

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