Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,141,380 A→T 5.2% L281H (CTC→CAC)  frwA ← putative PTS enzyme: Hpr, enzyme I and II components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,141,3800AT 5.2% 376.7 / 6.2 125L281H (CTC→CAC) frwAputative PTS enzyme: Hpr, enzyme I and II components
Reads supporting (aligned to +/- strand):  ref base A (65/53);  new base T (2/5);  total (67/58)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 125 total aligned reads displayed.

ATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAG  >  NC_000913/4141233‑4141526
                                                                                                                                                   |                                                                                                                                                  
aTCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAgg                                                                                                                                                  <  2:1515050/150‑1 (MQ=255)
aTCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAgg                                                                                                                                                  <  1:306894/150‑1 (MQ=255)
   gTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAAGGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAAc                                                                                                                                               >  1:342018/1‑150 (MQ=255)
     ttgttgCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTAGGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAActg                                                                                                                                             <  2:9562/150‑1 (MQ=255)
      tgttgCTGCCACGGAGTAAGGGCATCAATATCCACACCATCCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAActgg                                                                                                                                            <  1:1077706/150‑1 (MQ=255)
       gttgCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCTTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCCTCGAActgga                                                                                                                                           <  2:670390/150‑2 (MQ=255)
        ttgCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCCTCGAActggat                                                                                                                                          <  1:59631/150‑3 (MQ=255)
          gctgcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAActggctgg                                                                                                                                        >  1:765473/1‑150 (MQ=255)
           ctgcCACGGAGTAATGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTggg                                                                                                                                       <  1:282766/150‑1 (MQ=255)
           ctgcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTggg                                                                                                                                       <  2:1279907/150‑1 (MQ=255)
             gcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTggggg                                                                                                                                     <  2:935271/150‑1 (MQ=255)
             gcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGGGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTggggg                                                                                                                                     <  1:1039000/150‑1 (MQ=255)
               cACGGAGTATGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCTAGGATCACCGTATGTGAGGGGGTGTCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTc                                                                                                                                   <  1:1704429/150‑1 (MQ=255)
               cACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTc                                                                                                                                   <  2:830209/150‑1 (MQ=255)
                  ggAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGt                                                                                                                                <  1:1129073/150‑1 (MQ=255)
                  ggAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGt                                                                                                                                >  2:1673703/1‑150 (MQ=255)
                       aaGGGCATCACTATCCACACCACCCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTtcatc                                                                                                                           >  2:1002031/1‑150 (MQ=255)
                       aaGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGTAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTtcatc                                                                                                                           <  2:1314726/150‑1 (MQ=255)
                       aaGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTtcatc                                                                                                                           <  2:1609945/150‑1 (MQ=255)
                          ggCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGc                                                                                                                        >  1:1102318/1‑150 (MQ=255)
                            cATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCa                                                                                                                      >  1:677719/1‑150 (MQ=255)
                             aTCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCGTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCAt                                                                                                                     >  1:693528/1‑150 (MQ=255)
                                    tCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATg                                                                                                              <  1:574064/150‑1 (MQ=255)
                                    tCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGATGGGGGTCAGTTCATCAGCCATACAAATg                                                                                                              >  1:219710/1‑150 (MQ=255)
                                    tCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCAAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATg                                                                                                              <  2:980841/150‑1 (MQ=255)
                                      cacaCCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGc                                                                                                            >  1:539787/1‑150 (MQ=255)
                                       acacCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCg                                                                                                           >  1:1617939/1‑150 (MQ=255)
                                        cacCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCTGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCgg                                                                                                          >  1:264575/1‑150 (MQ=255)
                                        cacCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCgg                                                                                                          >  1:1600783/1‑150 (MQ=255)
                                         accaaccaGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATGTTTATCCAGTTCGAGGAACTGGCGGGGGGTCAGTTCGTCAGCCATACCAATGGCggg                                                                                                         >  1:1708481/1‑150 (MQ=255)
                                          ccaaccaGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATGCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGc                                                                                                        <  2:1201577/150‑1 (MQ=255)
                                          ccaaccaGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGc                                                                                                        <  1:143438/150‑1 (MQ=255)
                                          ccaaccaGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGc                                                                                                        >  2:917970/1‑150 (MQ=255)
                                              ccaGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCAGACAAATGGCGGGCTGCg                                                                                                    >  2:1473507/1‑150 (MQ=255)
                                                 gCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCa                                                                                                 >  1:903304/1‑150 (MQ=255)
                                                    ttGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGCGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGtt                                                                                              >  2:246769/1‑150 (MQ=255)
                                                     tGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGttt                                                                                             >  2:64715/1‑150 (MQ=255)
                                                      ggAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTg                                                                                            >  1:1527575/1‑150 (MQ=255)
                                                      ggAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTg                                                                                            >  2:1454731/1‑150 (MQ=255)
                                                             tGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTg                                                                                     >  2:1016467/1‑150 (MQ=255)
                                                                 cgaacgGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCgg                                                                                 >  1:796188/1‑150 (MQ=255)
                                                                  gaacgGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCggg                                                                                >  2:1592919/1‑150 (MQ=255)
                                                                     cgGGAAAGGATCACCGTATGTGAGGTGGGGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAg                                                                             <  2:1670885/150‑1 (MQ=255)
                                                                         cAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGtt                                                                         <  2:1718596/150‑1 (MQ=255)
                                                                              aTCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCa                                                                    <  1:1673703/150‑1 (MQ=255)
                                                                                    gaaTGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAg                                                              <  1:164804/148‑1 (MQ=255)
                                                                                    gTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAg                                                              >  1:203951/1‑150 (MQ=255)
                                                                                            ggtggtGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTg                                                      <  1:1432105/150‑1 (MQ=255)
                                                                                            ggtggtGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTg                                                      <  2:734100/150‑1 (MQ=255)
                                                                                             gtggtgCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGCACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGGCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTgg                                                     <  1:1002031/150‑1 (MQ=255)
                                                                                             gtggtgCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTgg                                                     <  1:1016467/150‑1 (MQ=255)
                                                                                             gtggtgCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAACCTTTCAACACCGTAGATTTGATgg                                                     >  2:1502070/1‑150 (MQ=255)
                                                                                               ggtgCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAg                                                   >  2:461084/1‑150 (MQ=255)
                                                                                                   ccgccgCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAAc                                               >  2:257446/1‑150 (MQ=255)
                                                                                                    cgccgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAAAGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACt                                              >  2:413195/1‑150 (MQ=255)
                                                                                                     gccgcTTTTGAGCAACAATCCTTGGAGGTGATTTTTATCCAGTTCGAGGAACAGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTg                                             <  2:686063/150‑1 (MQ=255)
                                                                                                       cgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGa                                           <  1:1051172/150‑1 (MQ=255)
                                                                                                         cTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAg                                         <  1:246769/150‑1 (MQ=255)
                                                                                                         cTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAg                                         >  1:1635288/1‑150 (MQ=255)
                                                                                                          ttttGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGc                                        <  2:789634/150‑1 (MQ=255)
                                                                                                          ttttGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGc                                        <  2:796188/150‑1 (MQ=255)
                                                                                                          ttttGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGc                                        <  2:986509/150‑1 (MQ=255)
                                                                                                            ttGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAg                                      <  2:1600783/150‑1 (MQ=255)
                                                                                                               aGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACg                                   >  1:1663410/1‑150 (MQ=255)
                                                                                                                gCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGt                                  <  2:1617939/150‑1 (MQ=255)
                                                                                                                gCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGt                                  >  1:190852/1‑150 (MQ=255)
                                                                                                                 caacaaTCATTGGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTGCATCAGCCAGACAACTGGCGGGCTTGGTCAGTTTGCCCGGTGCCGGGAAGCCATGCTCCGCGTAGAGTGTCGGGAGTACCTGGAAGCAGACGTc                                 >  2:698621/1‑150 (MQ=255)
                                                                                                                  aacaaTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTcg                                <  2:195455/150‑1 (MQ=255)
                                                                                                                    caaTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTcgcg                              >  2:663623/1‑150 (MQ=255)
                                                                                                                    caaTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTcgcg                              >  2:1425040/1‑150 (MQ=255)
                                                                                                                      aTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTa                            >  1:1471230/1‑150 (MQ=255)
                                                                                                                       tCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTAc                           >  2:72236/1‑150 (MQ=255)
                                                                                                                          tttGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTc                        <  2:765473/150‑1 (MQ=255)
                                                                                                                            tGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCa                      >  1:1241309/1‑150 (MQ=255)
                                                                                                                             gAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGTCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAg                     <  1:1454731/150‑1 (MQ=255)
                                                                                                                              cggTGATTTGTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCAGCAGCCATACAAATGGCGGGCTGCGGCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAgg                    <  2:893322/149‑1 (MQ=255)
                                                                                                                              aGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAgg                    <  2:1102318/150‑1 (MQ=255)
                                                                                                                               ggTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAggg                   >  2:904666/1‑150 (MQ=255)
                                                                                                                               ggTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTAAGTTTGACCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTTGAGTAACTGGAAGCAGAAGTCGAGTACGTCCAtgg                   >  2:351625/1‑150 (MQ=255)
                                                                                                                                gTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGCCGTCGCGTACGTCCAGGGc                  >  2:1649199/1‑150 (MQ=255)
                                                                                                                                gTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGc                  <  2:619125/150‑1 (MQ=255)
                                                                                                                                  gATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAAATGGAAGCAGACGTCGCGTACGTCCAGGGCAc                >  2:1239131/1‑150 (MQ=255)
                                                                                                                                   aTTTTTATCCAGTTCGAGGAACTGGCGGGGGGGCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACg               <  1:117733/150‑1 (MQ=255)
                                                                                                                                    tttttATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGt              >  1:630236/1‑150 (MQ=255)
                                                                                                                                    tttttATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGt              >  1:964570/1‑150 (MQ=255)
                                                                                                                                      tttATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTc            >  1:48362/1‑150 (MQ=255)
                                                                                                                                      tttATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTc            >  1:1236405/1‑150 (MQ=255)
                                                                                                                                      tttATCCAGTTCGAGGAACTGGCTGGGCGTCATATCATCAGCCATACAAATGGAGTGCTGCGTCAGTTTGCCCGGTGCAGTTAAGCGTTGCTCACCGTAGATTTCCTGGAGTAACAGCAAGCAGAAGTCGCGTACTTCCAGGGCACGTTc            >  1:159780/1‑150 (MQ=255)
                                                                                                                                       ttATCCTGTTCGAGGAACTGGCTGGGTGTTTGTTCATCTTGCATACAATTGGCGGGCTGCGTCAGTTTGCGCGGTGCCTGGAATCGTTGCTCACCGTTGTTTTGCTGGTGTAACTGGAAGCAGACGTCGCGTACGTGCAGGGCACGTTCt           <  1:355859/150‑1 (MQ=255)
                                                                                                                                       ttATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCt           <  2:1527575/150‑1 (MQ=255)
                                                                                                                                         agccAGTTCCTCGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTTCGTCCAGGGCACGTTCTTg         >  2:1031149/3‑150 (MQ=255)
                                                                                                                                         aTCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTg         <  2:903304/150‑1 (MQ=255)
                                                                                                                                         aTCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTg         >  2:1376663/1‑150 (MQ=255)
                                                                                                                                         aTCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTAGCTCACCGTAGCATTGCTGGAGTAACTGGAAGCAGACCTCGCGTCCGTCCAGCGCACGTTCTTg         >  1:1035894/1‑150 (MQ=255)
                                                                                                                                          tCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGc        >  2:366049/1‑150 (MQ=255)
                                                                                                                                             aGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCa                 >  1:144416/1‑138 (MQ=255)
                                                                                                                                             aGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCa                 <  2:144416/138‑1 (MQ=255)
                                                                                                                                             aGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAgg     <  2:190852/150‑1 (MQ=255)
                                                                                                                                                tCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAg  >  1:715088/1‑150 (MQ=255)
                                                                                                                                                tCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAg  >  1:867288/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAG  >  NC_000913/4141233‑4141526

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