Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,450,956 A→C 5.6% intergenic (+38/‑102) ytfQ → / → ytfR galactofuranose ABC transporter periplasmic binding protein/putative sugar ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,450,9560AC 5.6% 342.0 / 6.3 123intergenic (+38/‑102)ytfQ/ytfRgalactofuranose ABC transporter periplasmic binding protein/putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (55/60);  new base C (2/6);  total (57/66)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 128 total aligned reads displayed.

ATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGT  >  NC_000913/4450807‑4451105
                                                                                                                                                     |                                                                                                                                                     
aTACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATaa                                                                                                                                                       >  1:1702534/1‑150 (MQ=255)
aTACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGGAACCCAAATCCACCCTCTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGACGAAAAATATGGGGTATGGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATaa                                                                                                                                                       <  2:718291/150‑1 (MQ=255)
  aCAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAgg                                                                                                                                                     <  1:1289503/150‑1 (MQ=255)
  aCAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAgg                                                                                                                                                     <  1:168184/150‑1 (MQ=255)
   cAAAAAAGACGGCACCATGCCTGAATAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGATATGCCTGATGCCGATTCGTCGGCCGGATAAGgc                                                                                                                                                    >  1:314649/1‑150 (MQ=255)
   cAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgc                                                                                                                                                    >  1:817812/1‑150 (MQ=255)
   cAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgc                                                                                                                                                    <  1:629312/150‑1 (MQ=255)
   cAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgc                                                                                                                                                    >  2:1199069/1‑150 (MQ=255)
     aaaaaGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgcgc                                                                                                                                                  >  2:2125/1‑150 (MQ=255)
       aaaGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACACTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGAAAAGGCGCTc                                                                                                                                                >  1:1407169/1‑150 (MQ=255)
        aaGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGgcgcgag                                                                                                                                               <  2:499984/150‑4 (MQ=255)
           aCGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgc                                                                                                                                            <  2:1150972/150‑1 (MQ=255)
            cGGCACCATGCCTGAAAAGCTGACGTTAACCAAATACACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                                                                                                                           >  1:193626/1‑150 (MQ=255)
            cGGCACCATGCCTGAAAAGCTGACGTTAACAAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                                                                                                                           <  2:766315/150‑1 (MQ=255)
                  cATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                                     <  1:595204/150‑1 (MQ=255)
                    tGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCTCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                                   <  2:519508/150‑1 (MQ=255)
                     gCCTGAAAAGCTGACGTTAACCAAAACCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGCGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                                                                                                  <  1:1462789/150‑1 (MQ=255)
                      ccTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCg                                                                                                                                 >  1:907124/1‑150 (MQ=255)
                             aGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGCGGGCTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGc                                                                                                                          <  1:1206245/150‑1 (MQ=255)
                              gCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCa                                                                                                                         >  2:1465145/1‑150 (MQ=255)
                               cTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCAc                                                                                                                        <  2:461538/150‑1 (MQ=255)
                               cTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCAc                                                                                                                        >  1:634792/1‑150 (MQ=255)
                                    gTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                                                                                                   >  1:309845/1‑150 (MQ=255)
                                    gTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGACGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                                                                                                   >  1:1483386/1‑150 (MQ=255)
                                      tAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                                                                                 >  2:1145919/1‑150 (MQ=255)
                                      tAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGcc                                                                                                                 <  1:1713729/150‑1 (MQ=255)
                                       aaCCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                                                                                                >  2:236407/1‑150 (MQ=255)
                                        aCCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                                                                               >  2:49871/1‑150 (MQ=255)
                                        aCCAAATCCACCCTTTACCTGCCTGAAACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATACGTAGGCCGGATAAGGCGCTCGCGCAGCATCCGGCGATGGTGCACTGAAGCCTg                                                                                                               >  1:221163/1‑150 (MQ=255)
                                             aTCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCg                                                                                                          >  2:19866/1‑150 (MQ=255)
                                             aTCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCg                                                                                                          <  2:709630/150‑1 (MQ=255)
                                              tCCAGCCTTTACCTGCCTGATACCGCAAAAGTAGAATTAGAGTAGAAGAAAAATATGGGGATTTGAGGGGTGCTATGACTGATGCCGAGTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATTGTGCACTGAAGCCTGATGCGa                                                                                                         <  2:955809/150‑1 (MQ=38)
                                                cACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGACGCCTGATGCGACg                                                                                                       >  1:316876/1‑150 (MQ=255)
                                                 aCCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATTGGGTATTGAGGGTTGCTATGCCTGCTGCCTATTCGTAGGCAGGATAAGGCGCTCGCGCCGCATCCGGAGATGGTGCACTGAAGCCTGATGCTACga                                                                                                      >  1:98172/1‑149 (MQ=255)
                                                    cTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                                                                                   <  2:55057/150‑1 (MQ=255)
                                                        aCCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                                                                                               <  2:1219990/150‑1 (MQ=255)
                                                         cctgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                                                                              >  1:321439/1‑150 (MQ=255)
                                                         cctgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                                                                              <  1:742969/150‑1 (MQ=255)
                                                         cctgcctgATACCGCAAAACTATAATTATACAAGATTAAAAATATGGGGTTTTGAGGGTTGCTATGCCTGATGCCGTTTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                                                                              <  1:672393/150‑1 (MQ=38)
                                                          ctgcctgATACCGCAAACTAACAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                                                                             <  1:1333726/150‑1 (MQ=255)
                                                          ctgcctgATACCGCAAAAGAAGAATTAGATAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                                                                             <  2:1529894/150‑1 (MQ=255)
                                                          ctgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                                                                             >  2:541061/1‑150 (MQ=255)
                                                          ctgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                                                                             >  2:1192258/1‑150 (MQ=255)
                                                           tgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                                                                            >  2:802389/1‑150 (MQ=255)
                                                           tgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                                                                            >  2:795653/1‑150 (MQ=255)
                                                           tgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                                                                            >  2:1332332/1‑150 (MQ=255)
                                                            gcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                                                                           <  2:634792/150‑1 (MQ=255)
                                                            gcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                                                                           >  1:1471000/1‑150 (MQ=255)
                                                              cttctcacGCAAAAGAAGACTTCGAGACCACGAAAACTCTGGGGTATTGAGGGGTGCTATGCCTGATGCCCAGGCGTAGGCCGGATACGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                                                                         <  1:1431755/143‑1 (MQ=255)
                                                              ctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                                                                         >  1:393304/1‑150 (MQ=255)
                                                               tgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                                                                                        >  2:378887/1‑150 (MQ=255)
                                                               tgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                                                                                        >  2:879487/1‑150 (MQ=255)
                                                                gATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATTCGACGCTTACCGCGTCTTAta                                                                                       >  2:545196/1‑149 (MQ=255)
                                                                    ccGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                                                                                   >  2:499767/1‑150 (MQ=255)
                                                                     cGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGCAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                                                                                  >  2:1322430/1‑150 (MQ=255)
                                                                     cGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                                                                                  <  2:103298/150‑1 (MQ=255)
                                                                     cGCAAAAGAAGAATTAGAGAAGAAGAAAAATAGGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                                                                                  <  1:989808/150‑1 (MQ=255)
                                                                            gaagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTCGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                                                                           >  2:862071/1‑150 (MQ=255)
                                                                              agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                                                                         <  1:668852/150‑1 (MQ=255)
                                                                              agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCAGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                                                                         <  2:323270/150‑1 (MQ=255)
                                                                               gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc                                                                        >  1:1565945/1‑150 (MQ=255)
                                                                                    agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt                                                                   >  2:287778/1‑150 (MQ=255)
                                                                                    agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt                                                                   >  2:753158/1‑150 (MQ=255)
                                                                                      agaagaagaAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt                                                                 <  1:1332332/150‑1 (MQ=255)
                                                                                              aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCggg                                                         <  1:2125/150‑1 (MQ=255)
                                                                                                     tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACAGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAAcc                                                  >  1:394421/1‑150 (MQ=255)
                                                                                                     tGGGGTAGTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAAcc                                                  <  1:1528734/150‑1 (MQ=255)
                                                                                                      ggggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCa                                                 >  1:392342/1‑150 (MQ=255)
                                                                                                        ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATg                                               >  1:1103209/1‑150 (MQ=255)
                                                                                                          tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATgac                                             <  1:1371893/150‑1 (MQ=255)
                                                                                                          tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATgac                                             <  1:1385574/150‑1 (MQ=255)
                                                                                                           aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATgacg                                            <  1:1465145/150‑1 (MQ=255)
                                                                                                              gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGAcgccc                                         >  2:966492/1‑150 (MQ=255)
                                                                                                                gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGAcgaccga                                       >  1:914265/1‑150 (MQ=255)
                                                                                                                  gtggCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACgaccgacc                                     <  2:1360115/147‑1 (MQ=255)
                                                                                                                  gctGCTATCTCGGATGCCGATACTTTGGCCGGAAAAGGCGCTCGCGCCGAATCCTGCGATGGTGCACTGAAGCCTGTTCCGACGATTACAGCGTCTTATCATGGCTACTGGTAGAACGCTTTACACCGGGGGAAACCATGACgaccgacc                                     <  1:340940/148‑1 (MQ=255)
                                                                                                                   ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCa                                    <  2:991775/150‑1 (MQ=255)
                                                                                                                     gCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAac                                  <  1:654120/150‑1 (MQ=255)
                                                                                                                      cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAaca                                 >  1:911065/1‑150 (MQ=255)
                                                                                                                      cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAaca                                 <  1:530936/150‑1 (MQ=255)
                                                                                                                        aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACAcc                               <  2:393304/150‑1 (MQ=255)
                                                                                                                        aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACAcc                               <  2:483809/150‑1 (MQ=255)
                                                                                                                          gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAg                             <  1:1673510/150‑1 (MQ=255)
                                                                                                                            cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGga                           <  2:1699444/150‑1 (MQ=255)
                                                                                                                                 gCCGATTCGTAGGCCGGAAACGGCGCGCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGAtcc                      <  2:289430/150‑1 (MQ=255)
                                                                                                                                  ccGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGAtcct                     >  2:7256/1‑150 (MQ=255)
                                                                                                                                          atcgCCGGATGCGGCGCGAGCGCCTTATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCCACTGGGAGCACGCTTTACACCTGCGGAAACCATGACGACCGACCAACACAAGGAGATCATCCGCACCg             >  2:1628929/4‑150 (MQ=255)
                                                                                                                                           aGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTGCCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGCCGACCGACCTACGCCAGGAGCTCCTCCGCACCgc            >  1:750403/1‑149 (MQ=38)
                                                                                                                                            cgCCGGATGCGGCGCGAGCGCCTTATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGaa           >  2:281572/2‑150 (MQ=255)
                                                                                                                                              ccGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAgg         <  2:817812/150‑1 (MQ=255)
                                                                                                                                               cGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGa        <  1:1283049/150‑1 (MQ=255)
                                                                                                                                                ggATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACGCCGGGGGATACCAGGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGag       >  1:1099067/1‑149 (MQ=255)
                                                                                                                                                ggATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGAt       <  2:1569219/150‑1 (MQ=255)
                                                                                                                                                 gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGAtt      >  2:848382/1‑150 (MQ=255)
                                                                                                                                                  aTAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTa     <  1:1157758/150‑1 (MQ=255)
                                                                                                                                                   tAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTaa    <  2:1565945/150‑1 (MQ=255)
                                                                                                                                                   tAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTaa    >  2:834754/1‑150 (MQ=255)
                                                                                                                                                    aaGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAg   <  1:753158/150‑1 (MQ=255)
                                                                                                                                                    aaGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAg   <  2:1664531/150‑1 (MQ=255)
                                                                                                                                                     aGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGt  >  2:663593/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGT  >  NC_000913/4450807‑4451105

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: