Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,450,962 T→G 6.7% intergenic (+44/‑96) 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,9620TG 6.7% 316.2 / 7.5 121intergenic (+44/‑96)ytfQ/ytfRgalactofuranose ABC transporter periplasmic binding protein/putative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (54/58);  new base G (3/6);  total (57/64)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

AAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAAT  >  NC_000913/4450813‑4451109
                                                                                                                                                     |                                                                                                                                                   
aaaaGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCt                                                                                                                                                     <  1:1721225/150‑1 (MQ=255)
 aaaGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACACTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGAAAAGGCGCTc                                                                                                                                                    >  1:1407169/1‑150 (MQ=255)
     aCGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgc                                                                                                                                                <  2:1150972/150‑1 (MQ=255)
      cGGCACCATGCCTGAAAAGCTGACGTTAACCAAATACACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                                                                                                                               >  1:193626/1‑150 (MQ=255)
      cGGCACCATGCCTGAAAAGCTGACGTTAACAAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                                                                                                                               <  2:766315/150‑1 (MQ=255)
       ggCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgccg                                                                                                                                              <  1:136600/150‑1 (MQ=255)
            cATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                                         <  1:595204/150‑1 (MQ=255)
              tGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCTCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                                       <  2:519508/150‑1 (MQ=255)
              tGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCgg                                                                                                                                       <  2:828855/150‑1 (MQ=255)
               gCCTGAAAAGCTGACGTTAACCAAAACCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGCGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                                                                                                      <  1:1462789/150‑1 (MQ=255)
                ccTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCg                                                                                                                                     >  1:907124/1‑150 (MQ=255)
                   gAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATg                                                                                                                                  <  1:1417902/150‑1 (MQ=255)
                       aGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTCGAGAAGAAGAAAAAAATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGc                                                                                                                              <  2:715044/150‑1 (MQ=255)
                       aGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGCGGGCTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGc                                                                                                                              <  1:1206245/150‑1 (MQ=255)
                        gCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCa                                                                                                                             >  2:1465145/1‑150 (MQ=255)
                         cTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCAc                                                                                                                            <  2:467893/150‑1 (MQ=255)
                         cTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCAc                                                                                                                            <  2:461538/150‑1 (MQ=255)
                         cTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCAc                                                                                                                            >  1:634792/1‑150 (MQ=255)
                         cTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCAc                                                                                                                            >  1:289430/1‑150 (MQ=255)
                              gTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                                                                                                       >  1:309845/1‑150 (MQ=255)
                              gTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                                                                                                       >  1:1243611/1‑150 (MQ=255)
                              gTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGACGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                                                                                                       >  1:1483386/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)
                                   ccAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                                                                                  >  1:1471788/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)
                                         ccACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGAc                                                                                                            >  2:363744/1‑150 (MQ=255)
                                          cACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGACGCCTGATGCGACg                                                                                                           >  1:316876/1‑150 (MQ=255)
                                           aCCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATTGGGTATTGAGGGTTGCTATGCCTGCTGCCTATTCGTAGGCAGGATAAGGCGCTCGCGCCGCATCCGGAGATGGTGCACTGAAGCCTGATGCTACga                                                                                                          >  1:98172/1‑149 (MQ=255)
                                             ccTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                                                                                        <  2:951406/150‑1 (MQ=255)
                                              cTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                                                                                       <  2:55057/150‑1 (MQ=255)
                                                  aCCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                                                                                                   <  2:1219990/150‑1 (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:1192258/1‑150 (MQ=255)
                                                    ctgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                                                                                 >  2:541061/1‑150 (MQ=255)
                                                     tgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                                                                                >  2:795653/1‑150 (MQ=255)
                                                      gcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                                                                               >  1:1471000/1‑150 (MQ=255)
                                                      gcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                                                                               <  2:634792/150‑1 (MQ=255)
                                                        cttctcacGCAAAAGAAGACTTCGAGACCACGAAAACTCTGGGGTATTGAGGGGTGCTATGCCTGATGCCCAGGCGTAGGCCGGATACGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                                                                             <  1:1431755/143‑1 (MQ=255)
                                                        ctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                                                                             >  1:393304/1‑150 (MQ=255)
                                                         tgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                                                                                            >  2:378887/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)
                                                               cGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                                                                                      <  2:1702534/150‑1 (MQ=255)
                                                               cGCAAAAGAAGAATTAGAGAAGAAGAAAAATAGGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                                                                                      <  1:989808/150‑1 (MQ=255)
                                                                     agaagaATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                                                                                <  2:1218511/150‑1 (MQ=255)
                                                                     agaagaATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                                                                                <  2:1362306/150‑1 (MQ=255)
                                                                      taagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                                                                               <  2:1566248/149‑1 (MQ=255)
                                                                      gaagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTCGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                                                                               >  2:862071/1‑150 (MQ=255)
                                                                       aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                                                                              <  1:633870/150‑1 (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:753158/1‑150 (MQ=255)
                                                                                agaagaagaAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt                                                                     <  1:1332332/150‑1 (MQ=255)
                                                                                               tGGGGTAGTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAAcc                                                      <  1:1528734/150‑1 (MQ=255)
                                                                                                ggggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCa                                                     >  1:392342/1‑150 (MQ=255)
                                                                                                 gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCAt                                                    <  1:49871/150‑1 (MQ=255)
                                                                                                  ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATg                                                   >  1:1103209/1‑150 (MQ=255)
                                                                                                    tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATgac                                                 <  1:1371893/150‑1 (MQ=255)
                                                                                                       tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATgacgac                                              <  1:764508/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)
                                                                                                               gCTATGCATGATGCCGATTCGTATGCCGGATAAGGCGCTCGCTCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAac                                      <  2:1067006/150‑1 (MQ=255)
                                                                                                                cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAacc                                     >  2:279827/1‑149 (MQ=255)
                                                                                                                cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAaca                                     >  1:911065/1‑150 (MQ=255)
                                                                                                                cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAaca                                     <  1:530936/150‑1 (MQ=255)
                                                                                                                  aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACAcc                                   <  2:483809/150‑1 (MQ=255)
                                                                                                                    gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAg                                 <  1:1673510/150‑1 (MQ=255)
                                                                                                                           gCCGATTCGTAGGCCGGAAACGGCGCGCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGAtcc                          <  2:289430/150‑1 (MQ=255)
                                                                                                                                   gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGCTGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGACACCATGACGCCCGACCAACGCCAGGAGATCCTCCGCAcc                  >  1:330689/1‑150 (MQ=255)
                                                                                                                                    atcgCCGGATGCGGCGCGAGCGCCTTATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCCACTGGGAGCACGCTTTACACCTGCGGAAACCATGACGACCGACCAACACAAGGAGATCATCCGCACCg                 >  2:1628929/4‑150 (MQ=255)
                                                                                                                                      cgCCGGATGCGGCGCGAGCGCCTTATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGaa               >  2:281572/2‑150 (MQ=255)
                                                                                                                                        ccGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAgg             <  2:817812/150‑1 (MQ=255)
                                                                                                                                        ccGGATAAGGCGCGCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTGCGCACCGAAgg             <  2:789520/150‑1 (MQ=255)
                                                                                                                                          ggATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACGCCGGGGGATACCAGGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGag           >  1:1099067/1‑149 (MQ=255)
                                                                                                                                          ggATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGAt           <  2:1569219/150‑1 (MQ=255)
                                                                                                                                           gatgcggCGCGAGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGAtt          >  2:1433058/6‑150 (MQ=255)
                                                                                                                                           gATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGAtt          >  2:848382/1‑150 (MQ=255)
                                                                                                                                            aTAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTa         <  1:1157758/150‑1 (MQ=255)
                                                                                                                                            aTAAGGCGCGCGCGCGGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGTATTa         <  1:690360/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)
                                                                                                                                                ggCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGCCGACCGACCAGCTCCAGGAGATCCTCCGCACCGGAGGAGTTAGtt     >  1:293243/1‑149 (MQ=38)
                                                                                                                                                ggCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCCACACCAGGAGATCCTCCGCACCGAAGGATTAAGTa     >  2:631296/1‑150 (MQ=255)
                                                                                                                                                   gcTCGCGCCGAATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGTAGCACGCTTTACACCGGGTGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAAt  <  1:1325225/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGGATTAAGTAAAT  >  NC_000913/4450813‑4451109

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