Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,944,233 G→T 8.9% intergenic (+34/+113) znuB → / ← ruvB zinc ABC transporter permease/ATP‑dependent DNA helicase, component of RuvABC resolvasome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,944,2330GT 8.9% 302.7 / 13.0 89intergenic (+34/+113)znuB/ruvBzinc ABC transporter permease/ATP‑dependent DNA helicase, component of RuvABC resolvasome
Reads supporting (aligned to +/- strand):  ref base G (43/38);  new base T (2/6);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGGTTTAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATT  >  NC_000913/1944091‑1944381
                                                                                                                                              |                                                                                                                                                    
tGGCGGTTTAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGcc                                                                                                                                                >  3:521802/1‑149 (MQ=255)
 ggCGGTTTAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGcc                                                                                                                                                >  2:438382/1‑148 (MQ=255)
    ggTTTAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTAt                                                                                                                                            <  4:762/149‑1 (MQ=255)
     gTTTAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTAt                                                                                                                                            <  1:810812/148‑1 (MQ=255)
     gTTTAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGCGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATc                                                                                                                                           <  4:878989/149‑1 (MQ=255)
      tttAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATAAGGCGCTCGCGccgcatc                                                                                                                                           <  1:632477/148‑6 (MQ=255)
        tAACCTTTTTCGCGGTTTACGATACGCCGGCGGGTGCGGTGGTGGTGCTCTGTACGGCGCTGTTATTTATTCTCAGTCTGATAAAAAAGTATTCCCGCTCTTCTATCGCTGAACACATTTCTCGGATGCGGCGATATCGCCTTGTCCGa                                                                                                                                        >  4:527853/1‑149 (MQ=255)
        tAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGa                                                                                                                                        >  4:37233/1‑149 (MQ=255)
        tAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATAAGGCGCTCGCGCCGCATCCg                                                                                                                                         <  1:495881/148‑1 (MQ=255)
           ccTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTCTTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATAAGGCGCTCGCGCCGCATCCGAcaa                                                                                                                                     <  4:539133/149‑3 (MQ=255)
           ccTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGAcc                                                                                                                                      >  2:64424/1‑148 (MQ=255)
             ttttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTa                                                                                                                                    >  1:468420/1‑148 (MQ=255)
             ttttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTAc                                                                                                                                   >  3:302221/1‑149 (MQ=255)
             ttttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATAAGGCGCTCGCGCCGCATCCGAcaaat                                                                                                                                   <  3:811901/149‑5 (MQ=255)
             ttttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATAAGGCGCTCGCGCCGCATCCGAcaaa                                                                                                                                    <  1:529492/148‑4 (MQ=255)
              tttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTAc                                                                                                                                   >  2:304249/1‑148 (MQ=255)
               ttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACg                                                                                                                                  <  1:71247/148‑1 (MQ=255)
               ttCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACg                                                                                                                                  <  2:236971/148‑1 (MQ=255)
                   gcgGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTc                                                                                                                              <  1:442510/148‑1 (MQ=255)
                       tttACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTAt                                                                                                                         >  3:674165/1‑149 (MQ=255)
                          aCGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTAtctc                                                                                                                      <  3:777244/149‑1 (MQ=255)
                            gATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGTATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTg                                                                                                                    >  4:837475/1‑149 (MQ=255)
                            gATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATctct                                                                                                                     <  2:413249/148‑1 (MQ=255)
                            gATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTg                                                                                                                    <  3:209362/149‑1 (MQ=255)
                             aTACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTg                                                                                                                    >  1:604948/1‑148 (MQ=255)
                                  ccGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGt                                                                                                               <  2:524431/148‑1 (MQ=255)
                                    ggcggGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTcc                                                                                                             >  2:490063/1‑148 (MQ=255)
                                    ggcggGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTcc                                                                                                             <  1:685227/148‑1 (MQ=255)
                                       gggTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTgg                                                                                                          <  1:392504/148‑1 (MQ=255)
                                           cggTCGGTGGTGCTATGTGCGGCAGTGTTATTTATTCTCAGTATGATTAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCAGATGCGGCGCGAGCGTCGTATCCGACCTACGGTTCGGTATCTCTGGGAGTCTTGGTAAGa                                                                                                     <  4:191392/147‑1 (MQ=255)
                                                 gtggtgCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGa                                                                                                <  2:260849/148‑1 (MQ=255)
                                                     tgCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAg                                                                                            <  2:179337/148‑1 (MQ=255)
                                                      gCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCg                                                                                          >  3:847412/1‑149 (MQ=255)
                                                        tATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGt                                                                                         <  2:211509/148‑1 (MQ=255)
                                                        tATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGt                                                                                         >  2:370105/1‑148 (MQ=255)
                                                        tATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTc                                                                                        <  4:553612/149‑1 (MQ=255)
                                                             gCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCAt                                                                                    <  2:875304/148‑1 (MQ=255)
                                                              cGGCAGGGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATGTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATc                                                                                   <  2:617039/148‑1 (MQ=255)
                                                               ggCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCa                                                                                  <  1:64424/148‑1 (MQ=255)
                                                                     gTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCata                                                                            >  1:666379/1‑148 (MQ=255)
                                                                       tatttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCCGGCactgg                                                                          >  1:301152/1‑144 (MQ=255)
                                                                        atttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTg                                                                         <  1:919896/148‑1 (MQ=255)
                                                                          ttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGcc                                                                       <  1:801956/148‑1 (MQ=255)
                                                                          ttattCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCa                                                                      <  3:448601/149‑1 (MQ=255)
                                                                               ctcAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGcc                                                                 >  3:223096/1‑149 (MQ=255)
                                                                                tcAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCg                                                                >  3:216166/1‑149 (MQ=255)
                                                                                 cAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTc                                                                                                                              >  2:346180/1‑86 (MQ=255)
                                                                                 cAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTc                                                                                                                              <  1:346180/86‑1 (MQ=255)
                                                                                 cAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCg                                                                >  1:416073/1‑148 (MQ=255)
                                                                                  aGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCgg                                                               <  2:705317/148‑1 (MQ=255)
                                                                                   gTATGATGAAAAATCAGGCCAGCTAATCTATCGCTCAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGAt                                                             >  4:876420/1‑149 (MQ=255)
                                                                                    tATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGAt                                                             >  1:643917/1‑148 (MQ=255)
                                                                                      tgatgaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATgcg                                                          >  4:610425/1‑149 (MQ=255)
                                                                                      tgatgaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATgcg                                                          >  3:741005/1‑149 (MQ=255)
                                                                                       gatgaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACAGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATgcgg                                                         >  4:507007/1‑149 (MQ=255)
                                                                                          gaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTaa                                                                                                       >  2:40661/1‑100 (MQ=255)
                                                                                          gaAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTaa                                                                                                       <  1:40661/100‑1 (MQ=255)
                                                                                           aaaaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGgcgc                                                      >  2:14687/1‑148 (MQ=255)
                                                                                           aaaaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGgcgc                                                      >  2:147052/1‑148 (MQ=255)
                                                                                            aaaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCgcga                                                    <  3:828301/149‑1 (MQ=255)
                                                                                             aaaGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCgcgag                                                   <  3:876420/149‑1 (MQ=255)
                                                                                                gCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGAATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCgcgagcg                                                 <  1:700428/148‑1 (MQ=255)
                                                                                                 cAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCgcgagcga                                                >  2:462706/1‑148 (MQ=255)
                                                                                                  aGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGAcc                                              >  3:564188/1‑149 (MQ=255)
                                                                                                      cAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAAt                                           <  2:416073/148‑1 (MQ=255)
                                                                                                       aGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGAGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATcc                                         <  4:674165/149‑1 (MQ=255)
                                                                                                        gCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCg                                        >  4:276859/1‑149 (MQ=255)
                                                                                                           aaTCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACt                                     <  4:345678/149‑1 (MQ=255)
                                                                                                            aTCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGc                                                                               <  2:870035/106‑1 (MQ=255)
                                                                                                            aTCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGc                                                                               >  1:870035/1‑106 (MQ=255)
                                                                                                             tCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACtt                                    <  2:162618/148‑1 (MQ=255)
                                                                                                               tATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGa                                       >  4:223607/1‑143 (MQ=255)
                                                                                                               tATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGa                                       <  3:223607/143‑1 (MQ=255)
                                                                                                                 tCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACgg                                >  1:896549/1‑148 (MQ=255)
                                                                                                                     tGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATc                                                                                   >  4:700627/1‑93 (MQ=255)
                                                                                                                     tGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATc                                                                                   <  3:700627/93‑1 (MQ=255)
                                                                                                                        acacATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCt                         >  2:432085/1‑148 (MQ=255)
                                                                                                                        acacATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCt                         >  2:715101/1‑148 (MQ=255)
                                                                                                                         cacaTTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTgg                       >  3:15344/1‑149 (MQ=255)
                                                                                                                          acaTTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTggc                      >  4:460005/1‑149 (MQ=255)
                                                                                                                            aTTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTggcgg                    >  4:177716/1‑149 (MQ=255)
                                                                                                                            aTTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTggcgg                    >  4:461576/1‑149 (MQ=255)
                                                                                                                              ttGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTggcggc                   >  1:661688/1‑148 (MQ=255)
                                                                                                                              ttGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTggcggc                   >  2:535787/1‑148 (MQ=255)
                                                                                                                               tGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGgcggcg                  <  2:661688/148‑1 (MQ=255)
                                                                                                                                    gATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATATCTGGCGGCGGTACGc            >  4:806818/1‑149 (MQ=255)
                                                                                                                                     ataaggcgCTCGCGCCGCATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATgcgg                                                         <  3:964/99‑1 (MQ=255)
                                                                                                                                      taaggcgCTCGCGCCGCATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGcc           >  1:525956/4‑148 (MQ=255)
                                                                                                                                        cggcgCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCaa         <  2:213924/148‑1 (MQ=255)
                                                                                                                                              gAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGGCGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCACTTCTGGCGGCGGTATGCCAAAGTGAtt  >  4:824891/1‑149 (MQ=255)
                                                                                                                                              gAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGAtt  <  3:603578/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TGGCGGTTTAACCTTTTCCGCGGTTTACGATACGCCGGCGGGTCCGTCGGTGGTGCTATGTGCGGCACTGTTATTTATTCTCAGTATGATGAAAAAGCAGGCCAGCTAATCTATCGCTGAACACATTTGTCGGATGCGGCGCGAGCGCCTTATCCGACCTACGGTTCGGTATCTCTGGTAGTCCTGGTAAGACGCGAACAGCGTCGCATCAGGCATATTGCCAGTGCCGGATGCGGCGCGAGCGACCAATCCGACTTACGGCATTTCTGGCGGCGTTATGCCAAAGTGATT  >  NC_000913/1944091‑1944381

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: