Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,584,864 T→A 6.1% L541F (TTA→TTT ydeP ← putative oxidoreductase YdeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,584,8640TA 6.1% 247.0 / 5.6 82L541F (TTA→TTTydePputative oxidoreductase YdeP
Reads supporting (aligned to +/- strand):  ref base T (39/38);  new base A (3/2);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

ACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTT  >  NC_000913/1584716‑1585013
                                                                                                                                                    |                                                                                                                                                     
aCTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTa                                                                                                                                                      <  2:1241374/150‑1 (MQ=255)
aCTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTa                                                                                                                                                      >  1:1127826/1‑150 (MQ=255)
  tctGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAac                                                                                                                                                    <  2:116598/150‑1 (MQ=255)
     ggCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAAcacg                                                                                                                                                 <  1:290080/150‑1 (MQ=255)
     ggCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGAGTTCAGCATTACACCGGCGGGTTTTAAcacg                                                                                                                                                 >  2:1281263/1‑150 (MQ=255)
          cacaGCTTCAATGTCATTGCGAATTCGATCATAATCTTCCACCAGATACTACCAGGCTACCACGCTCTGGGGTAGTGATGCCTCCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATAACACCGGCGGGTTTTAAcaagccacg                                                                                                                                            >  2:1436311/1‑150 (MQ=255)
          cacaGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTAAAAcccgccggt                                                                                                                                            <  1:935903/150‑4 (MQ=255)
          cacaGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTAAAAcccgccggt                                                                                                                                            <  1:407355/150‑4 (MQ=255)
            caGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCAcgcg                                                                                                                                          <  2:267611/150‑1 (MQ=255)
            caGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCAcgcg                                                                                                                                          >  2:408210/1‑150 (MQ=255)
              gCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAg                                                                                                                                        <  2:90742/150‑1 (MQ=255)
              gCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAg                                                                                                                                        >  1:548957/1‑150 (MQ=255)
              gCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAg                                                                                                                                        <  1:199873/150‑1 (MQ=255)
                  cAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCAt                                                                                                                                    <  2:213236/150‑1 (MQ=255)
                    aTGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGa                                                                                                                                  >  1:771630/1‑150 (MQ=255)
                      gTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGTGTTTTAACACGCCACGCGAGGCATGAat                                                                                                                                <  2:1577221/150‑1 (MQ=255)
                         aTTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAatcat                                                                                                                             >  1:1419449/1‑150 (MQ=255)
                            gCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGa                                                                                                                          >  2:78805/1‑150 (MQ=255)
                                 tGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                                                                                                     <  1:95306/150‑1 (MQ=255)
                                 tGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTCCGCGATTCCCCCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTg                                                                                                                     >  2:221269/1‑150 (MQ=255)
                                   cGATCATATTCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGaa                                                                                                                   >  2:1613516/1‑150 (MQ=255)
                                   cGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGaa                                                                                                                   >  2:1057693/1‑150 (MQ=255)
                                      tCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATcc                                                                                                                <  2:288147/150‑1 (MQ=255)
                                                ccaccaGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTa                                                                                                      <  1:1302591/150‑1 (MQ=255)
                                                     aGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGcc                                                                                                 <  1:951857/150‑1 (MQ=255)
                                                        tACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTgc                                                                                              >  1:1071427/1‑150 (MQ=255)
                                                         aCTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTgcg                                                                                             <  1:548863/150‑1 (MQ=255)
                                                          cTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTgcgc                                                                                            <  2:1419449/150‑1 (MQ=255)
                                                            cccAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCCCGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCAc                                                                                          >  2:590882/1‑150 (MQ=255)
                                                            cccAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCAc                                                                                          <  2:125926/150‑1 (MQ=255)
                                                                 gCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGttt                                                                                     <  1:446459/150‑1 (MQ=255)
                                                                  cTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGtttt                                                                                    <  1:711298/150‑1 (MQ=255)
                                                                  cTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGATACCGCCTGCGCACCGtttt                                                                                    >  2:847490/1‑150 (MQ=255)
                                                                    aCCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGtttttt                                                                                  <  1:499670/150‑1 (MQ=255)
                                                                       aCGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGc                                                                               >  2:672586/1‑150 (MQ=255)
                                                                       aCGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGc                                                                               >  2:60325/1‑150 (MQ=255)
                                                                          ctctGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACAATGAATCCTCAACGGTTACCGCCTGCGCACAGTTTTTTTGCATg                                                                            <  2:644258/150‑1 (MQ=255)
                                                                                ggTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAAtc                                                                      >  1:106240/1‑150 (MQ=255)
                                                                                   aGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCTCACCGTTTTTTTGCATGTCAATCTCg                                                                   <  2:771630/150‑1 (MQ=255)
                                                                                      gctgcCTGCGCGATTCCCGCGACCACTTCACACTCTGATTTCAGCATTACACCGGCGTGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTa                                                                <  1:844177/150‑1 (MQ=255)
                                                                                      gctgcCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTa                                                                >  1:1525379/1‑150 (MQ=255)
                                                                                        tgcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACg                                                              >  2:227816/1‑150 (MQ=255)
                                                                                        tgcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACg                                                              >  1:471222/1‑150 (MQ=255)
                                                                                            tgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTccc                                                          >  1:774349/1‑150 (MQ=255)
                                                                                               gccaTTCCCTCGACCACTGCACACTATGATTTAAGCATTACACAGGCGGGTTTTAACACGCAACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAgg                                                       >  2:594238/4‑150 (MQ=255)
                                                                                                  aTTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAcc                                                    <  2:106240/150‑1 (MQ=255)
                                                                                                  aTTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAcc                                                    >  1:579706/1‑150 (MQ=255)
                                                                                                  aTTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAcc                                                    >  2:1304759/1‑150 (MQ=255)
                                                                                                  aTTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAcc                                                    >  1:1087865/1‑150 (MQ=255)
                                                                                                    tCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCgg                                                  >  1:1649977/1‑150 (MQ=255)
                                                                                                     cccGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGc                                                 <  1:507316/150‑1 (MQ=255)
                                                                                                      ccGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCa                                                >  2:1566771/1‑150 (MQ=255)
                                                                                                      ccGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCa                                                >  2:1572589/1‑150 (MQ=255)
                                                                                                         cgACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGaa                                             >  2:411973/1‑150 (MQ=255)
                                                                                                             cACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAAtata                                         <  1:1057693/150‑1 (MQ=255)
                                                                                                              aCTGCACACAGTGATTTAAGCATTACACCGGCGGGTTTTAACACGCCCCGCGAGGCATGAATCATCGGCATTGCATCCTCAACGGTTTCCGCCTGCGCACCGTTTTTTTTCATGTGAATCACGCTACGTCCCAGGCCCGGCAGAATATAg                                        >  2:1027800/1‑150 (MQ=255)
                                                                                                               cTGCACACTCTGATTTCAGCATTACACCGGCTGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGc                                       <  1:1456546/150‑1 (MQ=255)
                                                                                                               cTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGc                                       >  2:1493932/1‑150 (MQ=255)
                                                                                                               cTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGc                                       >  1:222341/1‑150 (MQ=255)
                                                                                                                 gCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAg                                                        <  1:55173/131‑1 (MQ=255)
                                                                                                                 gCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAg                                                        >  2:55173/1‑131 (MQ=255)
                                                                                                                 gCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTa                                     <  2:343516/150‑1 (MQ=255)
                                                                                                                   acacTGTGATTTCAGCATTACACCGGCGGGTTATAACACGCCACGCGAGGCATGAATCATCGTCATTGAATCCTCATCTGTTAACGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCGAGGAGCGGCAGAATATAGCTAtt                                   <  2:1577258/150‑2 (MQ=255)
                                                                                                                   acacTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATg                                   <  1:408210/150‑1 (MQ=255)
                                                                                                                        ctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCaa                                                                        >  2:421939/1‑108 (MQ=255)
                                                                                                                        ctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCaa                                                                        <  1:421939/108‑1 (MQ=255)
                                                                                                                           aTTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCgg                           <  1:221269/150‑1 (MQ=255)
                                                                                                                             ttCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTc                         >  2:1175049/1‑150 (MQ=255)
                                                                                                                                aGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAAc                      <  1:227816/150‑1 (MQ=255)
                                                                                                                                 gCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACa                     >  1:1316812/1‑150 (MQ=255)
                                                                                                                                  cATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAg                    <  1:636520/150‑1 (MQ=255)
                                                                                                                                   aTTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGa                   >  1:581204/1‑150 (MQ=255)
                                                                                                                                    ttACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGAt                  <  1:1199105/150‑1 (MQ=255)
                                                                                                                                     tACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATg                 <  1:672586/150‑1 (MQ=255)
                                                                                                                                       gcgtggcgtGTTAAAACCCGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATgag               >  1:1299214/5‑150 (MQ=255)
                                                                                                                                       gcgtggcgtGTTAAAACCCGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATgag               >  2:1513479/5‑150 (MQ=255)
                                                                                                                                       cacCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATgag               <  1:219132/150‑1 (MQ=255)
                                                                                                                                        cgtggcgtGTTAAAACCCGCCACGCGAGGCATGAATCATCGACATTGAATCATCAACGGTTACCGCCTGCGCACCGTTTTTTTGCAAGTCAATCTCGCTACGACCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGAAgaga              >  2:911987/4‑150 (MQ=255)
                                                                                                                                         ccGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATg                 <  2:899964/146‑1 (MQ=255)
                                                                                                                                         ccGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATg                 >  1:899964/1‑146 (MQ=255)
                                                                                                                                         ccGGCGGGTTTTAACACGCCACGCGAGGCAAGAATCATCGACATTGAATCCTCAACGGTTACCGACTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATCCCGTGCGGTAAACAGATGagag             >  2:1173414/1‑150 (MQ=255)
                                                                                                                                               ggTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTa       <  1:1436311/150‑1 (MQ=255)
                                                                                                                                                 ttttAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAg     <  2:1649977/150‑1 (MQ=255)
                                                                                                                                                 ttttAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAg     <  2:1279421/150‑1 (MQ=255)
                                                                                                                                                    tAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCtt  >  1:179652/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTT  >  NC_000913/1584716‑1585013

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