Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,707,923 G→A 5.8% R278C (CGC→TGC)  rseB ← anti‑sigma E factor, binds RseA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,707,9230GA 5.8% 277.2 / 3.2 85R278C (CGC→TGC) rseBanti‑sigma E factor, binds RseA
Reads supporting (aligned to +/- strand):  ref base G (42/38);  new base A (2/3);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGT  >  NC_000913/2707786‑2708071
                                                                                                                                         |                                                                                                                                                    
cccACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGAt                                                                                                                                           >  8:34133/1‑149 (MQ=255)
 ccACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATc                                                                                                                                          >  3:137899/1‑149 (MQ=255)
 ccACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATc                                                                                                                                          <  4:344271/149‑1 (MQ=255)
 ccACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATc                                                                                                                                          >  6:62385/1‑149 (MQ=255)
  cACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCg                                                                                                                                         <  8:336179/149‑1 (MQ=255)
    ctctTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGt                                                                                                                                       <  7:271630/149‑1 (MQ=255)
    ctctTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGt                                                                                                                                       >  1:278514/1‑149 (MQ=255)
     tctTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGa                                                                                                                                            >  8:312814/1‑143 (MQ=255)
     tctTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGa                                                                                                                                            <  7:312814/143‑1 (MQ=255)
     tctTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTg                                                                                                                                      >  8:68940/1‑149 (MQ=255)
       tttGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCt                                                                                                                                    >  1:387934/1‑149 (MQ=255)
       tttGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCt                                                                                                                                    >  4:144865/1‑149 (MQ=255)
         tGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTc                                                                                                                                   <  2:387934/148‑1 (MQ=255)
             cATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATgg                                                                                                                              >  5:366515/1‑149 (MQ=255)
               ttGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCg                                                                                                                            >  3:372114/1‑149 (MQ=255)
                       ccccGAACTTAATATTCTCGGCAAGGCGTTTCGCCGTTTGCGGCTGCAGTTCCCCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACGGACGGTTCTGCGTCCGGGGCGCAACATCTGATCGGTGCTCGATGGCGTAGcgcg                                                                                                                     <  2:140222/148‑1 (MQ=255)
                       ccccGAACGAAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGGTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCgg                                                                                                                    <  3:242347/149‑1 (MQ=255)
                        cccGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGt                                                                                                                   >  1:316475/1‑149 (MQ=255)
                             aCTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACg                                                                                                              <  7:68940/149‑1 (MQ=255)
                             aCTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCACAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACg                                                                                                              >  4:339158/1‑149 (MQ=255)
                             aCTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCACAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACg                                                                                                              >  4:344869/1‑149 (MQ=255)
                                   tatTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTAcc                                                                                                        <  6:369347/149‑1 (MQ=255)
                                        tcGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCCCGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACTTCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCgag                                                                                                     >  1:32447/1‑147 (MQ=255)
                                        tcGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGCTGGCGTAGCGCGGTTAACGTTTACCgag                                                                                                     >  1:394888/1‑147 (MQ=255)
                                        tcGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCgag                                                                                                     >  7:321535/1‑147 (MQ=255)
                                        tcGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCgag                                                                                                     >  7:403050/1‑147 (MQ=255)
                                          ggCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGaa                                                                                                   <  7:243994/147‑1 (MQ=255)
                                          ggCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAg                                                                                                  <  1:106252/148‑1 (MQ=255)
                                              aTGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGaa                                                                                             >  4:245478/1‑149 (MQ=255)
                                                 cGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATTGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGaata                                                                                           <  6:270138/148‑1 (MQ=255)
                                                  gTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGaataat                                                                                         >  4:248541/1‑149 (MQ=255)
                                                   tttCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATc                                                                                        >  4:96843/1‑149 (MQ=255)
                                                      cGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGt                                                                                     >  5:50644/1‑149 (MQ=255)
                                                        ccGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCg                                                                                   >  2:152352/1‑149 (MQ=255)
                                                                ggcggcAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATagag                                                                           <  1:211514/149‑1 (MQ=255)
                                                                ggcggcAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATagag                                                                           <  5:134400/149‑1 (MQ=255)
                                                                      aGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGa                                                                     >  3:128341/1‑149 (MQ=255)
                                                                       gTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGAt                                                                    >  7:49077/1‑149 (MQ=255)
                                                                           aCCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGa                                                                <  2:316475/149‑1 (MQ=255)
                                                                           aCCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGa                                                                <  5:267422/149‑1 (MQ=255)
                                                                            ccGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGAt                                                               <  3:248541/149‑1 (MQ=255)
                                                                            ccGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGAt                                                               <  1:152352/149‑1 (MQ=255)
                                                                            ccGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGAGTCCGGTGCGCAACAGCTGATCGGTGCTCGATGGCGTAGCTCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGAAGTGATTCGa                                                                >  5:92808/1‑148 (MQ=255)
                                                                            ccGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTGCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGAt                                                               <  2:211092/149‑1 (MQ=255)
                                                                             cGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATa                                                                              >  4:244440/1‑133 (MQ=255)
                                                                             cGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATa                                                                              <  3:244440/133‑1 (MQ=255)
                                                                             cGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGa                                                                >  1:19597/1‑147 (MQ=255)
                                                                             cGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATa                                                              <  3:144865/149‑1 (MQ=255)
                                                                             cGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGGTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATa                                                              <  5:248371/149‑1 (MQ=255)
                                                                              gACAATGGTGATTTCGGCGTTATCACGTACGATTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTAGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAg                                                             <  6:366515/149‑1 (MQ=255)
                                                                               aCAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAgg                                                            >  5:173453/1‑149 (MQ=255)
                                                                                   tGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCAtg                                                        <  8:403050/149‑1 (MQ=255)
                                                                                   tGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCAtg                                                        >  8:64074/1‑149 (MQ=255)
                                                                                     gTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCAtgtt                                                      <  4:372114/149‑1 (MQ=255)
                                                                                       gATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCAGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGAGATTCGATAGGCatgatga                                                    >  2:2605/1‑148 (MQ=255)
                                                                                        aTTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCAtgttgt                                                    <  6:93354/148‑1 (MQ=255)
                                                                                        aTTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCACAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCAtgttgt                                                    <  3:344869/148‑1 (MQ=255)
                                                                                        aTTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCACAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCAtgttgt                                                    <  3:339158/148‑1 (MQ=255)
                                                                                           tCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCa                                                 <  2:278514/148‑1 (MQ=255)
                                                                                                 ttATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCgg                                          <  3:245478/149‑1 (MQ=255)
                                                                                                 ttATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCgg                                          >  6:218170/1‑149 (MQ=255)
                                                                                                 ttATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCgg                                          >  8:175205/1‑149 (MQ=255)
                                                                                                 ttATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCgg                                          >  1:167479/1‑149 (MQ=255)
                                                                                                    tCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAg                                       >  1:47153/1‑149 (MQ=255)
                                                                                                    tCACGTACGCTTGTACTGACGGTTCTGCGGCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTa                                        <  7:164190/148‑1 (MQ=255)
                                                                                                        gTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGa                                   <  1:160249/149‑1 (MQ=255)
                                                                                                            gCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTc                               >  2:32176/1‑149 (MQ=255)
                                                                                                                ggaCTGACGGGTCTGCGTCCGGTGCGCAACATCTGATCGGTGCGCGATGGCGTAGCGCGGTTAACGTGTCCCGAGAAGCTGTATACTCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGa                             <  6:374072/145‑1 (MQ=255)
                                                                                                                gTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGGCGGTAGCGGACGTCGa                             <  6:371593/147‑1 (MQ=255)
                                                                                                                   cTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGactact                        >  1:24995/1‑149 (MQ=255)
                                                                                                                    tGACGGTTCTGCGTCCGGTGCGCAAAATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATACGTCGGAATAGAGACTAGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGCCTACTg                       >  6:46392/1‑149 (MQ=255)
                                                                                                                    tGACGGTTCTGCGTCAGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTAGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTg                       >  2:347128/1‑149 (MQ=255)
                                                                                                                      aCGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTg                       <  1:32176/147‑1 (MQ=255)
                                                                                                                        ggTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGAt                                                               <  5:150936/105‑1 (MQ=255)
                                                                                                                        ggTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGAt                                                               >  6:150936/1‑105 (MQ=255)
                                                                                                                        ggTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGaaa                   <  5:62385/149‑1 (MQ=255)
                                                                                                                         gTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAAc                  >  4:308025/1‑149 (MQ=255)
                                                                                                                              gCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCg              >  2:137218/1‑148 (MQ=255)
                                                                                                                              gCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCg              <  1:137218/148‑1 (MQ=255)
                                                                                                                               cGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGc             <  4:148074/148‑1 (MQ=255)
                                                                                                                               cGTCCGGTGCACAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCg              <  7:216290/147‑1 (MQ=255)
                                                                                                                                gTCCGGTGCGCAACATCTAATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCt            >  8:192362/1‑148 (MQ=255)
                                                                                                                                   cGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTaaaa        >  4:4897/1‑149 (MQ=255)
                                                                                                                                   cGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTaaaa        <  2:47153/149‑1 (MQ=255)
                                                                                                                                   cGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAAGAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTaaaa        >  3:338714/1‑149 (MQ=255)
                                                                                                                                        cgcAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCTTCGGTAGCGGACGTCGAATACTGGAAACTTCGCTAAACCCCtg   >  6:281520/1‑149 (MQ=255)
                                                                                                                                        agcAACATCTGATCGGTGCTCGATGGCGTAGCGCGTTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACTCGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCt    <  6:231531/147‑1 (MQ=255)
                                                                                                                                         gcAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCtgt  >  6:244394/1‑149 (MQ=255)
                                                                                                                                         gcAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCtgt  <  7:64074/149‑1 (MQ=255)
                                                                                                                                         gcAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCAATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCtgt  >  6:135121/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
CCCACTCTTTGATCATTGCGCTGCCCCGAACTTAATATTCTCGGCAATGCGTTTCGCCGTTTGCGGCGGCAGTTCACCGACAATGGTGATTTCGGCGTTATCACGTACGCTTGTACTGACGGTTCTGCGTCCGGTGCGCAACATCTGATCGGTGCTCGATGGCGTAGCGCGGTTAACGTTTACCGAGAAGCTGAATAATCCGTCGGAATAGAGACGTGATTCGATAGGCATGTTGTCCATCGTCGGTAGCGGACGTCGACTACTGGAAACTTCGCTAAAACCCTGT  >  NC_000913/2707786‑2708071

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