Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,559,335 A→C 5.8% F64L (TTT→TTG)  iadA ← isoaspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,559,3350AC 5.8% 2319.8 / 37.0 672F64L (TTT→TTG) iadAisoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (271/351);  new base C (48/2);  total (319/353)
Fisher's exact test for biased strand distribution p-value = 1.77e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 691 total aligned reads displayed.

TCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGC  >  NC_000913/4559186‑4559484
                                                                                                                                                     |                                                                                                                                                     
tCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATa                                                                                                                                                       >  2:1399122/1‑150 (MQ=255)
 ccGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATaa                                                                                                                                                      <  2:5298172/150‑1 (MQ=255)
   gggTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATaaaa                                                                                                                                                    >  2:1751837/1‑150 (MQ=255)
   gggTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATaaaa                                                                                                                                                    >  1:4677277/1‑150 (MQ=255)
   gggTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATaaaa                                                                                                                                                    >  1:3381658/1‑150 (MQ=255)
   gggTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATaaaa                                                                                                                                                    >  2:159366/1‑150 (MQ=255)
   gggTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATaaaa                                                                                                                                                    <  1:3181985/150‑1 (MQ=255)
      tGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCt                                                                                                                                                 <  2:220618/150‑1 (MQ=255)
      tGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCt                                                                                                                                                 <  1:2667259/150‑1 (MQ=255)
       ggCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGAACTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTg                                                                                                                                                <  1:6329741/150‑1 (MQ=255)
         cGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATCAAACCTggg                                                                                                                                              >  1:6186721/1‑150 (MQ=255)
          gagagaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGc                                                                                                                                             >  2:3356631/1‑150 (MQ=255)
          gagagaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGc                                                                                                                                             <  1:4481971/150‑1 (MQ=255)
           agagaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATCACACCTGGGCa                                                                                                                                            >  1:6963677/1‑150 (MQ=255)
             agaTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCaga                                                                                                                                          <  2:4589868/150‑1 (MQ=255)
               aTAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAgg                                                                                                                                        <  1:2485271/150‑1 (MQ=255)
                tAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGa                                                                                                                                       <  1:7062155/150‑1 (MQ=255)
                     tCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGc                                                                                                                                  <  1:2387297/150‑1 (MQ=255)
                     tCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGc                                                                                                                                  >  1:6751050/1‑150 (MQ=255)
                     tCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGc                                                                                                                                  <  2:1210711/150‑1 (MQ=255)
                       ggTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGccc                                                                                                                                <  2:890624/150‑1 (MQ=255)
                        gTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCa                                                                                                                               <  1:7647713/150‑1 (MQ=255)
                         tGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCag                                                                                                                              <  1:3088808/150‑2 (MQ=255)
                          gCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACt                                                                                                                             <  1:3957772/150‑1 (MQ=255)
                          gCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACt                                                                                                                             <  1:6846517/150‑1 (MQ=255)
                          gCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACt                                                                                                                             <  2:5896257/150‑1 (MQ=255)
                              agcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTgaga                                                                                                                         >  2:6882508/1‑150 (MQ=255)
                              agcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTgaga                                                                                                                         <  2:3872824/150‑1 (MQ=255)
                               gcagACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGAt                                                                                                                        <  2:8010437/150‑1 (MQ=255)
                                  gACCAACCACTGACGTGCCGCCCGCGTCCGTCCGGCGACTTAGCGCCACGTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAACTGGCCGTGTTGATCAATCAAACCTGGGCAGAGGATCTGCCCACTGAGATCGa                                                                                                                     <  2:7530138/150‑1 (MQ=255)
                                  gACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGa                                                                                                                     <  2:2932904/150‑1 (MQ=255)
                                   accaaccaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCCATACAACCTGGGCAGAGGATCTGCCCACTGAGATCGAc                                                                                                                    >  1:1007506/1‑150 (MQ=255)
                                   accaaccaCTGACGTGACCCCCGCTACCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGAGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTCAGATCGAc                                                                                                                    >  2:8584793/1‑150 (MQ=255)
                                      aaccaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAAc                                                                                                                 >  1:8162925/1‑150 (MQ=255)
                                      aaccaCTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAAc                                                                                                                 >  1:8935890/1‑150 (MQ=255)
                                                 tgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATCCCACCTGGGCAGAGGATCTGCCCACGGAGATCGACAACCGTGCAGTTCg                                                                                                      >  1:4088628/1‑150 (MQ=255)
                                                 tgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCg                                                                                                      <  1:3285031/150‑1 (MQ=255)
                                                 tgacgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCg                                                                                                      >  2:4300504/1‑150 (MQ=255)
                                                    cgCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTa                                                                                                   <  1:2878586/150‑1 (MQ=255)
                                                     gCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTGCGGTAc                                                                                                  >  1:9180164/1‑150 (MQ=255)
                                                       ccGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACaa                                                                                                >  2:6187228/1‑150 (MQ=255)
                                                       ccGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACaa                                                                                                <  2:1294027/150‑1 (MQ=255)
                                                         gCTTCCGTCAGGCGACGTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAACTGGACGTGTTGATCAATCAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATg                                                                                              <  2:4896540/150‑1 (MQ=255)
                                                               gTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGaa                                                                                        >  2:3497397/1‑150 (MQ=255)
                                                               gTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGaa                                                                                        <  1:7430242/150‑1 (MQ=255)
                                                               gTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGaa                                                                                        <  1:6541412/150‑1 (MQ=255)
                                                               gTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGaa                                                                                        >  1:5270302/1‑150 (MQ=255)
                                                                   ggCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGa                                                                                    >  1:8791851/1‑150 (MQ=255)
                                                                    gCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGat                                                                                   <  2:2722381/150‑1 (MQ=255)
                                                                         ttAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGc                                                                              >  2:78430/1‑150 (MQ=255)
                                                                         ttAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGc                                                                              >  2:3061609/1‑150 (MQ=255)
                                                                         ttAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGc                                                                              >  2:7685045/1‑150 (MQ=255)
                                                                         ttAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGc                                                                              >  2:4437100/1‑150 (MQ=255)
                                                                            gcgcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTgg                                                                           >  2:5517611/1‑150 (MQ=255)
                                                                            gcgcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTgg                                                                           >  2:8517860/1‑150 (MQ=255)
                                                                              gcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCa                                                                         >  2:7427377/1‑150 (MQ=255)
                                                                              gcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCa                                                                         >  2:5781632/1‑150 (MQ=255)
                                                                                cACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAAc                                                                       >  2:7089893/1‑150 (MQ=255)
                                                                                        gCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATaa                                                               >  2:1483775/1‑150 (MQ=255)
                                                                                          gTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAAtt                                                             <  1:8388876/150‑1 (MQ=255)
                                                                                           tGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAAttt                                                            <  2:3674642/150‑1 (MQ=255)
                                                                                               cgtcgtGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCGGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATATTTTTGcc                                                        <  2:8752607/150‑1 (MQ=255)
                                                                                                 tcgtGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGt                                                      >  2:6747430/1‑150 (MQ=255)
                                                                                                   gtGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTa                                                    <  2:5023600/150‑1 (MQ=255)
                                                                                                      ggACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcg                                                 >  1:8057094/1‑150 (MQ=255)
                                                                                                      ggACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcg                                                 >  2:3995004/1‑150 (MQ=255)
                                                                                                       gACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGGTCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcga                                                >  2:6272043/1‑150 (MQ=255)
                                                                                                        aCCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcgac                                               <  2:9115887/150‑1 (MQ=255)
                                                                                                         ccTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATACAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGcgacg                                              >  1:4783680/1‑150 (MQ=255)
                                                                                                           tGCTTCGCCGCCACCGCCAATCAAATGGACTTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCTGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAg                                            <  1:5585026/150‑1 (MQ=255)
                                                                                                             cTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATACAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGa                                          >  1:6076548/1‑150 (MQ=255)
                                                                                                                   cgccaccgccaATCAAATGGACGTGTTGATCAATCCAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGc                                    >  1:1952779/1‑150 (MQ=255)
                                                                                                                   cgccaccgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGc                                    >  2:2298761/1‑150 (MQ=255)
                                                                                                                    gccaccgccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGAGATTGCTGGCAACGGCGATCATTTTGCCGTTAGCGACGAGGACCTCGCa                                   >  2:8290922/1‑150 (MQ=255)
                                                                                                                      caccgccaATCAACTGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCaaa                                 <  1:7160238/150‑1 (MQ=255)
                                                                                                                         cgccaATCAAAGGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTc                              <  1:979509/150‑1 (MQ=255)
                                                                                                                         cgcaaATCAACTGGACGTGTTGATCAATAAAACCTGGGCCGAGGATCTGCCCTCTGAGATCCAAAACCGTGCTGTTCGGTACAATTTCAGAAGGTATATTGCTGGCAACGCCGATAATTTTGCCGTTAGCGACGAGGTCATCGCAAATTc                              <  1:8517860/150‑1 (MQ=255)
                                                                                                                           ccaATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTccc                            <  2:2372373/150‑1 (MQ=255)
                                                                                                                            caATCAAATGGACGTGATGATCAATAACACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACGGTGCAGTTCGGTAGAATGTCTGAAGGGATAATGCTGGCAACGGCGATACTTTTGCCGTTAGCGTCGAGGACATCGCAAATTcccc                           <  1:8584793/150‑1 (MQ=255)
                                                                                                                             aaTCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCg                          >  1:1218851/1‑150 (MQ=255)
                                                                                                                              aTCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGa                         <  1:159595/150‑1 (MQ=255)
                                                                                                                               tCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGAt                        >  2:6585715/1‑150 (MQ=255)
                                                                                                                               tCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGAt                        >  1:6411430/1‑150 (MQ=255)
                                                                                                                                     ggACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCg                  <  2:7201454/150‑1 (MQ=255)
                                                                                                                                        cGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGgcg               <  1:576197/150‑1 (MQ=255)
                                                                                                                                        cGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGgca               <  1:4891674/150‑2 (MQ=255)
                                                                                                                                          tgtTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCa             >  1:1158014/1‑150 (MQ=255)
                                                                                                                                          tgtTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCa             <  1:7821019/150‑1 (MQ=255)
                                                                                                                                              gATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACa         >  1:4290882/1‑150 (MQ=255)
                                                                                                                                               aTCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACaa        >  2:1133796/1‑150 (MQ=255)
                                                                                                                                               aTCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACaa        >  2:4170636/1‑150 (MQ=255)
                                                                                                                                                 cAATACAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAAt      >  1:36155/1‑150 (MQ=255)
                                                                                                                                                 cAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACACCCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAAt      <  2:8057129/150‑1 (MQ=255)
                                                                                                                                                 cAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAAt      <  1:8704940/150‑1 (MQ=255)
                                                                                                                                                 cAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAAt      <  1:5423504/150‑1 (MQ=255)
                                                                                                                                                   aTAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATgc    >  1:21301/1‑150 (MQ=255)
                                                                                                                                                   aTAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATgc    >  2:2032511/1‑150 (MQ=255)
                                                                                                                                                   aTAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTAGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATgc    >  1:5049751/1‑150 (MQ=255)
                                                                                                                                                    tAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATgcg   >  2:4792658/1‑150 (MQ=255)
                                                                                                                                                     aaaaCCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATgcgc  <  1:1399122/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCAAATTCCCCGATCTTCCGGCGCATACAAATGCGC  >  NC_000913/4559186‑4559484

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