Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 245,789 A→C 5.7% L6* (TTA→TGA) ‡ yafK ← L,D‑transpeptidase domain‑containing protein YafK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913245,7890AC 5.7% 419.7 / 15.1 141L6* (TTA→TGA) ‡yafKL,D‑transpeptidase domain‑containing protein YafK
Reads supporting (aligned to +/- strand):  ref base A (74/59);  new base C (5/3);  total (79/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 143 total aligned reads displayed.

CGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGCCTGTCGGCGCAGGGTTT  >  NC_000913/245640‑245937
                                                                                                                                                     |                                                                                                                                                    
ggAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTa                                                                                                                                                      <  1:365644/149‑1 (MQ=255)
cGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTa                                                                                                                                                      >  1:314120/1‑150 (MQ=255)
  aGCGTACGTTCTTCCTTGAAGATCTGTATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAAt                                                                                                                                                    <  2:1889639/150‑1 (MQ=255)
     gTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCg                                                                                                                                                 >  1:5260/1‑150 (MQ=255)
        cGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGAtt                                                                                                                                              >  2:2305282/1‑150 (MQ=255)
        cGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGAtt                                                                                                                                              >  1:700073/1‑150 (MQ=255)
        cGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGAtt                                                                                                                                              <  1:996692/150‑1 (MQ=255)
           tcttcCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCATTAATTCTTGCGAtgttg                                                                                                                                           <  2:865385/150‑2 (MQ=255)
              tcCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGc                                                                                                                                        >  2:2063076/1‑150 (MQ=255)
                cTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCAt                                                                                                                                      >  1:2101648/1‑150 (MQ=255)
                 ttGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCATTAATTCTTGCGAtgttgttga                                                                                                                                       >  1:1080092/1‑143 (MQ=255)
                 ttGAAGATCTGGATGTAGACAGGGGATCCAATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATa                                                                                                                                     >  2:792966/1‑150 (MQ=255)
                 ttGAAGATCTGGATGTAGACAGGGGATCACATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCATTAATTCTTGCGAtgttgttga                                                                                                                                       <  2:1080092/148‑6 (MQ=255)
                   gAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACt                                                                                                                                   <  1:11143/150‑1 (MQ=255)
                    aaGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTg                                                                                                                                  <  2:1283934/150‑1 (MQ=255)
                      gATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCt                                                                                                                                <  2:934385/150‑1 (MQ=255)
                      gATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCt                                                                                                                                <  1:2360723/150‑1 (MQ=255)
                      gATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCt                                                                                                                                <  1:908811/150‑1 (MQ=255)
                        tCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTa                                                                                                                              >  1:551177/1‑150 (MQ=255)
                         cTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTAt                                                                                                                             <  2:2035/150‑1 (MQ=255)
                          tGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCCGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTAtt                                                                                                                            >  2:2363212/1‑150 (MQ=255)
                          tGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTAtt                                                                                                                            >  2:764000/1‑150 (MQ=255)
                           ggATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTc                                                                                                                           >  1:1443190/1‑150 (MQ=255)
                             aTGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCt                                                                                                                         >  1:1392709/1‑150 (MQ=255)
                               gTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTg                                                                                                                       <  2:19846/150‑1 (MQ=255)
                                tAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGc                                                                                                                                        >  1:669265/1‑132 (MQ=255)
                                  gACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATa                                                                                                                    >  1:1722069/1‑150 (MQ=255)
                                  gACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATa                                                                                                                    >  1:1074028/1‑150 (MQ=255)
                                  gACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATa                                                                                                                    >  2:955592/1‑150 (MQ=255)
                                      ggggATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAAc                                                                                                                >  2:1247809/1‑150 (MQ=255)
                                      ggggATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAAc                                                                                                                <  2:721658/150‑1 (MQ=255)
                                      ggggATCCCATCAACTGCCGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAAc                                                                                                                <  1:1488915/150‑1 (MQ=255)
                                       gggATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAACCGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACg                                                                                                               >  1:448300/1‑150 (MQ=255)
                                       gggATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACg                                                                                                               >  1:722080/1‑150 (MQ=255)
                                       gggATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACg                                                                                                               >  2:2063197/1‑150 (MQ=255)
                                        ggATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACgg                                                                                                              >  1:2010309/1‑150 (MQ=255)
                                         gATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGt                                                                                                             <  2:1789501/150‑1 (MQ=255)
                                         gATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGCGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGt                                                                                                             >  1:761768/1‑150 (MQ=255)
                                               aTCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTacaca                                                                                                       <  2:1392709/150‑1 (MQ=255)
                                                tCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTacacac                                                                                                      >  1:70869/1‑150 (MQ=255)
                                                 cAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACg                                                                                                     <  2:2155253/150‑1 (MQ=255)
                                                    ctgctgCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCa                                                                                                  <  1:1501570/150‑1 (MQ=255)
                                                     tgctgcTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAg                                                                                                 >  2:685360/1‑150 (MQ=255)
                                                     tgctgcTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAg                                                                                                 <  1:347548/150‑1 (MQ=255)
                                                      gctgctTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAgg                                                                                                >  2:1341348/1‑150 (MQ=255)
                                                            tataCTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCa                                                                                          >  2:11152/1‑150 (MQ=255)
                                                               aCTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCaaaa                                                                                       >  1:1653189/1‑150 (MQ=255)
                                                                 tctTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATg                                                                                     <  1:1802342/150‑1 (MQ=255)
                                                                  ctTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATgg                                                                                    >  1:2337937/1‑150 (MQ=255)
                                                                     tGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATc                                                                                 <  2:93405/150‑1 (MQ=255)
                                                                     tGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATc                                                                                 >  1:1118323/1‑150 (MQ=255)
                                                                       cTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATccc                                                                               <  2:407620/150‑1 (MQ=255)
                                                                       cTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATccc                                                                               >  1:198885/1‑150 (MQ=255)
                                                                       cTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATccc                                                                               >  1:1004666/1‑150 (MQ=255)
                                                                        tCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCa                                                                              <  1:1145271/150‑1 (MQ=255)
                                                                         cacaGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCaa                                                                             >  2:594583/1‑150 (MQ=255)
                                                                             ggCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCg                                                                         >  2:967916/1‑150 (MQ=255)
                                                                                gTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGaa                                                                      <  1:2186179/150‑1 (MQ=255)
                                                                                  ggTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATa                                                                    >  2:595093/1‑150 (MQ=255)
                                                                                     agagCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATGAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTc                                                                 <  1:970077/150‑1 (MQ=255)
                                                                                        gCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTgg                                                              <  2:5260/150‑1 (MQ=255)
                                                                                          tACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCtt                                                                                                                        <  2:2221153/90‑1 (MQ=255)
                                                                                          tACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCtt                                                                                                                        >  1:2221153/1‑90 (MQ=255)
                                                                                          tACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGAt                                                            <  1:1742569/150‑1 (MQ=255)
                                                                                              gctgcCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTgg                                                        <  2:2235811/150‑1 (MQ=255)
                                                                                              gctgcCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTgg                                                        <  1:1742/150‑1 (MQ=255)
                                                                                               ctgcCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGa                                                       <  2:1647141/150‑1 (MQ=255)
                                                                                               ctgcCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGa                                                       <  1:1162066/150‑1 (MQ=255)
                                                                                                tgcCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGaa                                                      <  2:336753/150‑1 (MQ=255)
                                                                                                 gcCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAg                                                     >  1:1318862/1‑150 (MQ=255)
                                                                                                   ccAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGc                                                   <  1:771721/150‑1 (MQ=255)
                                                                                                      gCAAACCGTCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTca                                                <  1:1228156/150‑1 (MQ=255)
                                                                                                           ccGGCAAATGAAAAGCACGGGATCAACAACATCGAAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGACAGGACGGCAAAATGTATAAAAAAACGGAATAGACTGGATTTGGAAGGCTAACATTa                                           >  2:1033332/1‑150 (MQ=255)
                                                                                                           caggCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCAGAATAGTCTGGATTTGGAAGGCTAACATTa                                           >  2:500807/3‑150 (MQ=255)
                                                                                                             ggCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATc                                         >  1:2319066/1‑150 (MQ=255)
                                                                                                             ggCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATc                                         >  1:1852038/1‑150 (MQ=255)
                                                                                                              gCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCa                                        <  1:902556/150‑1 (MQ=255)
                                                                                                                  aTGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCaa                                    >  2:2327939/1‑150 (MQ=255)
                                                                                                                     aaaCGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAag                                 <  2:1852038/150‑1 (MQ=255)
                                                                                                                     aaaCGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAag                                 <  2:1849441/150‑1 (MQ=255)
                                                                                                                        cGCACGGGATCAACAACATCGCAAGAATTAATGCTATTTTACGCATACTGCTTATTCCTTGATAAAAGGGTTACACACGCCAGGACTTCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGt                              <  1:603601/150‑1 (MQ=255)
                                                                                                                        cGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGt                              >  1:1019461/1‑150 (MQ=255)
                                                                                                                          cACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGttt                            >  2:1308302/1‑150 (MQ=255)
                                                                                                                            cGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTa                          >  1:1696345/1‑150 (MQ=255)
                                                                                                                             gggATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTAc                         <  1:1767211/150‑1 (MQ=255)
                                                                                                                               gATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACTGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGc                       <  1:16095/150‑1 (MQ=255)
                                                                                                                                  gcgtaaaaTCGCATTAATTCTTGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAg                    >  1:1420060/8‑150 (MQ=255)
                                                                                                                                  caacaaCATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAg                    >  2:673646/1‑150 (MQ=255)
                                                                                                                                   aacaacATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGc                   >  2:2451585/1‑150 (MQ=255)
                                                                                                                                    gtaaaaTCGCATTAATTCTTGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGcc                  >  1:2437523/6‑150 (MQ=255)
                                                                                                                                    acaacaTCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGcc                  >  2:885587/1‑150 (MQ=255)
                                                                                                                                         aTCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGTAGGCTCACATTATAACCAAAAGAGTTTTACGCAATCCTGTCg             <  2:198885/150‑1 (MQ=255)
                                                                                                                                         aTCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTGACATTATCACCAAAAGAGTTTTACGCAAGCCTGTCg             >  1:1038954/1‑150 (MQ=255)
                                                                                                                                          tCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGCCTGTCgg            <  2:2101648/150‑1 (MQ=255)
                                                                                                                                           cGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGCCTGTCGgc           >  2:2469244/1‑150 (MQ=255)
                                                                                                                                                gAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGCCTGTCGGCGCAgg      <  2:551177/150‑1 (MQ=255)
                                                                                                                                                 aaTTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGCCTGTCGGCGCAggg     <  1:1322092/150‑1 (MQ=255)
                                                                                                                                                   ttAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTGTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGCCTGTCGGCGCAGGGtt   <  2:1722069/150‑1 (MQ=255)
                                                                                                                                                    tAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGCCTGTCGGCGCAGGGttt  >  1:492601/1‑150 (MQ=255)
                                                                                                                                                    tAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGCCTGTCGGCGCAGGGttt  >  2:1903120/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAGCGTACGTTCTTCCTTGAAGATCTGGATGTAGACAGGGGATCCCATCAACTGCTGCTTATACTCTTTGCTCACAGGCGTGGTAGAGCTACTGCTGCCCAGCAAACCGGCAAATGAAACGCACGGGATCAACAACATCGCAAGAATTAATGCGATTTTACGCATACTGCTTATTCCTTGATAAAACGGTTACACACGCCAGGACGGCAAAATGGATCCCAAATCGGAATAGTCTGGATTTGGAAGGCTCACATTATCACCAAAAGAGTTTTACGCAAGCCTGTCGGCGCAGGGTTT  >  NC_000913/245640‑245937

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