Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,707,906 G→A 5.2% T283T (ACC→ACT rseB ← anti‑sigma factor stabilizing protein RseB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,707,9060GA 5.2% 243.9 / 2.1 76T283T (ACC→ACTrseBanti‑sigma factor stabilizing protein RseB
Reads supporting (aligned to +/- strand):  ref base G (33/39);  new base A (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

CTGCCCGTTTTGCCAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGA  >  NC_000913/2707757‑2708052
                                                                                                                                                     |                                                                                                                                                  
cTGCCCGTTTTGCCAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACg                                                                                                                                                    <  1:306257/150‑1 (MQ=255)
   cccGTTTTGCCAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGtt                                                                                                                                                 <  1:744080/150‑1 (MQ=255)
      gTTTTGCCAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTg                                                                                                                                              <  1:44225/150‑1 (MQ=255)
           gCCAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTcc                                                                                                                                         <  2:834746/150‑1 (MQ=255)
           gCCAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTcc                                                                                                                                         <  2:670584/150‑1 (MQ=255)
           gCCAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACAGTTCTGCGTcc                                                                                                                                         >  1:1271833/1‑150 (MQ=255)
           gCCAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACAGTTCTGCGTcc                                                                                                                                         >  1:1276998/1‑150 (MQ=255)
             cAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCgg                                                                                                                                       <  2:96731/150‑1 (MQ=255)
             cAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCgg                                                                                                                                       >  2:595732/1‑150 (MQ=255)
             cAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCgg                                                                                                                                       >  1:869173/1‑150 (MQ=255)
                gagaCGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTgc                                                                                                                                    <  1:196112/150‑1 (MQ=255)
                  gacgacGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTgcgc                                                                                                                                  >  2:775737/1‑150 (MQ=255)
                   acgacgGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCa                                                                                                                                 <  2:1005029/150‑1 (MQ=255)
                   acgacgGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCa                                                                                                                                 <  2:474738/150‑1 (MQ=255)
                        ggTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATc                                                                                                                            >  1:1322384/1‑150 (MQ=255)
                        ggTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTCCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATc                                                                                                                            <  2:901960/150‑1 (MQ=255)
                             cccACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCTCAACATCTGATc                                                                                                                       <  2:1150429/150‑1 (MQ=255)
                              ccACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCg                                                                                                                      >  1:352858/1‑150 (MQ=255)
                               cACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCgg                                                                                                                     <  1:337212/150‑1 (MQ=255)
                                 ctctTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTg                                                                                                                   <  1:1059788/150‑1 (MQ=255)
                                 cgctTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTg                                                                                                                   <  2:1279432/148‑1 (MQ=255)
                                     ttGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCg                                                                                                               >  1:80577/1‑150 (MQ=255)
                                       gATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGAt                                                                                                             >  1:188819/1‑150 (MQ=255)
                                       gATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGAt                                                                                                             <  2:1133057/150‑1 (MQ=255)
                                         tCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATgg                                                                                                           <  2:891181/150‑1 (MQ=255)
                                         tCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATgg                                                                                                           >  2:1353163/1‑150 (MQ=255)
                                          cATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGc                                                                                                          >  2:829507/1‑150 (MQ=255)
                                             tGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTa                                                                                                       >  1:811099/1‑150 (MQ=255)
                                              gcgcTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAg                                                                                                      >  1:8746/1‑150 (MQ=255)
                                                   gCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCgg                                                                                                 <  1:1094709/150‑1 (MQ=255)
                                                           cTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGtt                                                                                         <  2:677631/150‑1 (MQ=255)
                                                             tAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGtt                                                                                                                                                 >  1:762408/1‑92 (MQ=255)
                                                             tAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGtt                                                                                                                                                 <  2:762408/92‑1 (MQ=255)
                                                             tAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTa                                                                                       >  1:1039000/1‑150 (MQ=255)
                                                             tAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGCTAACGTTTa                                                                                       <  1:603362/150‑1 (MQ=255)
                                                               atatTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTAcc                                                                                     <  2:1252696/150‑1 (MQ=255)
                                                                 atTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTAc                                                                                      >  2:246602/1‑147 (MQ=255)
                                                                 atTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTAc                                                                                      <  1:246602/147‑1 (MQ=255)
                                                                  ttCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCgag                                                                                  <  2:572310/150‑1 (MQ=255)
                                                                  ttCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCgag                                                                                  <  2:1202496/150‑1 (MQ=255)
                                                                        gCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTg                                                                            >  2:720582/1‑150 (MQ=255)
                                                                        gCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTg                                                                            >  1:773548/1‑150 (MQ=255)
                                                                          aaTGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGaa                                                                          <  1:829507/150‑1 (MQ=255)
                                                                           aTGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGaat                                                                         <  1:595732/150‑1 (MQ=255)
                                                                            tGCGTTTCGCCGTTTCCGCCGGCAGTTCACCGACAATGGTGATTTCGCCGTTATCACCTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACAACTGATCGGTGCTCGATGCCGTACCGCGGTTAACGTTTACCGAGAAGCTGaata                                                                        >  2:792197/1‑150 (MQ=255)
                                                                               gTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATc                                                                     <  1:456095/150‑1 (MQ=255)
                                                                                tttCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATcc                                                                    <  2:157112/150‑1 (MQ=255)
                                                                                     ccGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACAGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCgg                                                               <  2:1271833/150‑1 (MQ=255)
                                                                                     ccGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACAGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCgg                                                               <  2:1276998/150‑1 (MQ=255)
                                                                                       gTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGaa                                                             <  2:1322384/150‑1 (MQ=255)
                                                                                          tGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATag                                                          <  2:811099/150‑1 (MQ=255)
                                                                                           gcggcgGCAGTTGACCGAGATTGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCTCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATaga                                                         <  1:1354760/150‑1 (MQ=255)
                                                                                             ggcggcAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAgaga                                                       >  2:243795/1‑150 (MQ=255)
                                                                                                  cAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGa                                                  <  2:188819/150‑1 (MQ=255)
                                                                                                  cAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGa                                                  >  1:281284/1‑150 (MQ=255)
                                                                                                     ttCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTc                                               >  1:1322440/1‑150 (MQ=255)
                                                                                                     ttCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCATCCGGTGCGCAACATCTGATCGGTGCTCCATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTc                                               >  2:870591/1‑150 (MQ=255)
                                                                                                      tCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCg                                              >  1:932797/1‑150 (MQ=255)
                                                                                                       cACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGa                                             >  1:842098/1‑150 (MQ=255)
                                                                                                               aTGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCAtg                                     >  1:279330/1‑150 (MQ=255)
                                                                                                                     aTTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTcc                               <  2:352858/150‑1 (MQ=255)
                                                                                                                     aTTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTcc                               >  1:804325/1‑150 (MQ=255)
                                                                                                                      tttCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTTTCCa                              >  1:1305392/1‑150 (MQ=255)
                                                                                                                       ttCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGAt                                            >  2:579154/1‑135 (MQ=255)
                                                                                                                       ttCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGAt                                            <  1:579154/135‑1 (MQ=255)
                                                                                                                          ggCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCAtcgt                          >  1:1327446/1‑150 (MQ=255)
                                                                                                                            cGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCAtcgtcg                        <  1:1353163/150‑1 (MQ=255)
                                                                                                                                   aCGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCgg                 >  2:897375/1‑150 (MQ=255)
                                                                                                                                   aCGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCgg                 <  1:243795/150‑1 (MQ=255)
                                                                                                                                         gCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCg           <  2:473250/150‑1 (MQ=255)
                                                                                                                                          cTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGa          >  1:1351897/1‑150 (MQ=255)
                                                                                                                                          cTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGa          >  2:961031/1‑150 (MQ=255)
                                                                                                                                           ttGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGac         <  2:869173/150‑1 (MQ=255)
                                                                                                                                               aCTGACGGTTCTGCGTCCGTTGCGCAAGTTGTGATCGGTGCTCGATGGCGTAGCGCGGTTAACTTTTACGGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGTCATGTTGTCCATCGTCGGTAGCGGACTTCGactact     <  2:804325/150‑1 (MQ=255)
                                                                                                                                               aCTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGactact     >  2:512625/1‑150 (MQ=255)
                                                                                                                                                cTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTg    <  2:8746/150‑1 (MQ=255)
                                                                                                                                                  gACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGa  >  1:1262801/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTGCCCGTTTTGCCAGGAGACGACGGTAGCCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGA  >  NC_000913/2707757‑2708052

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: