Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,196,230 C→T 14.3% intergenic (‑131/+102) thiC ← / ← rsd phosphomethylpyrimidine synthase/regulator of sigma D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,196,2300CT14.3% 148.9 / 14.8 70intergenic (‑131/+102)thiC/rsdphosphomethylpyrimidine synthase/regulator of sigma D
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base T (6/4);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

TGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGC  >  NC_000913/4196088‑4196371
                                                                                                                                              |                                                                                                                                             
tgttgCGGACATAGCTCATTCGAGAAAGTTATGGACGGGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCcgc                                                                                                                                               <  2:863356/143‑1 (MQ=255)
  ttgCAGACATAGCTCATTCCAAAAAGGTAAGGACGTGGCTGGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGACTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGGTCCGCgg                                                                                                                                             >  2:318195/1‑143 (MQ=255)
   tgCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGa                                                                                                                                            >  1:986932/1‑143 (MQ=255)
    gCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGAt                                                                                                                                           >  1:289240/1‑143 (MQ=255)
      aGACATAGCTGATTCCAAAAAGTTAAGGCCGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATcc                                                                                                                                         <  1:879800/143‑1 (MQ=255)
      aGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGGTCAGGTTCCGCGGATcc                                                                                                                                         >  1:95904/1‑143 (MQ=255)
      aGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATcc                                                                                                                                         >  1:1053971/1‑143 (MQ=255)
         cATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGa                                                                                                                                      <  1:989674/143‑1 (MQ=255)
              cTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACACTTGTTGTGTTCGCAGGTATTAGCCTGATCAGGTTCCGTGGATCCCGAATaaa                                                                                                                                 >  1:668594/1‑143 (MQ=255)
               tCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGTGGATCCCGAATAAAc                                                                                                                                >  2:980626/1‑143 (MQ=255)
                   gccAAAACGTTAAGGCCGTGGCTTGTCAGACGACGGATGACGCAAGAGCCGATCGCCGCAGGGGCGATGCGATAGCAGACTAACTCTTGTTCCCTTCGCGGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtc                                                                                                                            <  2:968275/142‑1 (MQ=255)
                    ccAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtct                                                                                                                           <  1:94946/143‑1 (MQ=255)
                     cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGCATAAACGGgctc                                                                                                                          >  2:805874/1‑143 (MQ=255)
                     cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtctc                                                                                                                          >  1:862865/1‑143 (MQ=255)
                      aaaaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGACTAAACGGTCTCa                                                                                                                         >  1:651612/1‑143 (MQ=255)
                      aaaaaGTGAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGGTGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCa                                                                                                                         <  2:208877/143‑1 (MQ=255)
                         gagctAAGGACCTGGCTTGTCAGACGACGGATGAAGCAAGAGTCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGcc                                                                                                                      <  2:812188/139‑1 (MQ=255)
                               atgACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTaa                                                                                                                <  2:621095/141‑1 (MQ=255)
                                    gtgtCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTATTTTTTGTTTTTTTTTCAGGTATTAGCCTGATCAGGTGCCGCGGAGCCCGAATAAACGGTCTCAGCCAGTTAAGGCAc                                                                                                           <  1:48870/143‑1 (MQ=255)
                                    gTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAACCGGTCTCCGCCAGTTCACGCAc                                                                                                           >  1:209134/1‑143 (MQ=255)
                                    gTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGACGCGATAGCAGATTAACACTTGTTCCCTTCGCAGGGATTAGCCTGATCACGTTCCGCGGATCCCGCATACACAGTCTCAGCCAGTTACGGCAc                                                                                                           >  2:199713/1‑143 (MQ=255)
                                       gCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTcc                                                                                                        >  1:571411/1‑143 (MQ=255)
                                           gTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGTGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACa                                                                                                    <  1:166533/143‑1 (MQ=255)
                                               gacgacGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTACGACAATCCGACAAGaa                                                                                                >  1:860455/1‑143 (MQ=255)
                                                acgacgGATGAAGCAAGAGCCGTTCGCCGCAGGGGCGATGCGATAGCAGATTACCTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGaaa                                                                                               <  2:683522/143‑1 (MQ=255)
                                                 cgacgGATGAAGCAAGGGAGGATCGCGGCAGGGGGGATGCGATAGCAGATTAACGCTTGGTCCCTTCGCAGGGATTAGCCTGATCAGGTGCCGCGGATCCCGAAGAAACGGTCGCAGCCAGTTAAGGCACTCCGACAAGAAAt                                                                                              <  1:810344/143‑1 (MQ=255)
                                                  gacgGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTGCGCGGATCCCGAATAAACGGTCTCAGCCGGGTAAGGCACTCCGGCAAGAAAtt                                                                                             >  1:208139/1‑143 (MQ=255)
                                                  gacgGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGAATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGaaact                                                                                             >  2:423068/1‑141 (MQ=255)
                                                  gacgGATGAAGCAAGAGACGATCGCCGCAGGGGAGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCACGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTCAGGCACCCCGACAAGaaaat                                                                                             >  2:145373/1‑141 (MQ=255)
                                                      gATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGGATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAAGCACTCCGACAAGAAATTAGcc                                                                                         >  2:191428/1‑143 (MQ=255)
                                                       aTGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGTGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGccc                                                                                        >  2:421940/1‑143 (MQ=255)
                                                             cAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCcaa                                                                                  >  1:188892/1‑143 (MQ=255)
                                                             cAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCaaa                                                                                  >  2:150790/1‑143 (MQ=255)
                                                             cAACAGACGATCGCCGCAGGGGCGAGGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCAAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCaaa                                                                                  <  2:49643/143‑1 (MQ=255)
                                                              aaGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAAt                                                                                 <  2:606425/143‑1 (MQ=255)
                                                               agagACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACGCCGACAAGAAATTAGACCCGCAAATg                                                                                >  2:312267/1‑143 (MQ=255)
                                                                      aTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAtt                                                                         >  2:666645/1‑143 (MQ=255)
                                                                      aTCGCCGCAGGGGCGAGGCGATAGCACATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCATGTTCCGCGGATCCCGAATAACCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAtt                                                                         <  1:430006/143‑1 (MQ=255)
                                                                           cgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATg                                                                    >  1:617131/1‑143 (MQ=255)
                                                                              aGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGTGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTaa                                                                 <  1:980626/143‑1 (MQ=255)
                                                                                   cgatgcgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAc                                                            <  2:542660/143‑1 (MQ=255)
                                                                                    gatgcgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGAACCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAAATACCCCCGCAAAGAGGGCATTGAACGTTAATTAcg                                                           >  2:309381/1‑143 (MQ=255)
                                                                                     atgcgagAGCAGATTATCTCTTGTTCCCTTCGCAGGTAATAGCCTGATCAGGTTCCGTGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAcgc                                                          <  1:884619/143‑1 (MQ=255)
                                                                                        cgatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAACCGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAGAAGGGGCATTGAATGTACATTCCGCGtt                                                       >  1:839635/1‑143 (MQ=255)
                                                                                         gatAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGTGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTa                                                      <  1:122901/143‑1 (MQ=255)
                                                                                              cAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGCATAAACGGTCTCAGCCAGTTAAGGCACTCCGCCAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAACTTACGCGTCAACAgc                                                 >  2:822198/1‑143 (MQ=255)
                                                                                               aGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAgcg                                                <  1:151749/143‑1 (MQ=255)
                                                                                                gATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAgcgc                                               >  2:733649/1‑143 (MQ=255)
                                                                                                gATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGAGGCATAGAAGCTAAATTACGCGTTAACAacgc                                               >  2:258732/1‑143 (MQ=255)
                                                                                                 aTTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGACAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAAAGAAAGTAAATTACGCGTGAACAGCGCa                                              >  1:392752/1‑143 (MQ=255)
                                                                                                     aCTCTTGTTCCCTTCGGAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACGCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAAc                                          <  1:860811/143‑1 (MQ=255)
                                                                                                       tctTGTTCCCTTCGCAGGTATTATCCTGATAAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTc                                        <  1:655280/143‑1 (MQ=255)
                                                                                                        ctTGTGACCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCa                                       <  2:95904/143‑1 (MQ=255)
                                                                                                         ttGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGTGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTACATTACGCGTAAACAGCGCAGAGCTCaa                                      >  2:545876/1‑143 (MQ=255)
                                                                                                         ttGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCaa                                      >  1:273326/1‑143 (MQ=255)
                                                                                                         ttGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCaa                                      >  2:401118/1‑143 (MQ=255)
                                                                                                         ttGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGGTAACAGCGCAGAACTCaa                                      >  1:317677/1‑143 (MQ=255)
                                                                                                          tGTTCCCTTCGCAGGTATTAGCCCGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAg                                     <  2:655904/143‑1 (MQ=255)
                                                                                                            ttCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGCTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCa                                   <  1:538254/143‑1 (MQ=255)
                                                                                                               ccTTCGCAGGTATTAGCCTGATCAGGTTCCGTGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTGAACAGCGCAGAACTCAAGCAGGa                                >  1:428303/1‑143 (MQ=255)
                                                                                                                   cGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGtt                            <  1:583109/143‑1 (MQ=255)
                                                                                                                    gCAGGTATTAGCCTGATCATGTTCCGCGGTTCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGttt                           <  2:597208/143‑1 (MQ=255)
                                                                                                                      aGGTATTAGCCTGATCAGGTACCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGa                         <  1:225107/143‑1 (MQ=255)
                                                                                                                         gatgaGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCGGAACTCAAGCAGGATGTTTGAcgc                      <  1:865082/139‑1 (MQ=255)
                                                                                                                               ccgcaTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCACGTAAGGCACTCCGACAAGAAATTAGCCCCGCAAACGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcgggcgg                <  2:208139/139‑1 (MQ=255)
                                                                                                                               ccTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACgcgggcgg                <  1:423068/143‑1 (MQ=255)
                                                                                                                                 tggtCAGGGTCCGCGGATCCCGAGTAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATGAGCCCCGCAAATGGGGCATTGAATGTAAAGTACGCGTGAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGgcggcg              <  2:659713/140‑1 (MQ=255)
                                                                                                                                  gcccAGGTCCCGCGCATCCAGAATAAACGGTCTCAGCCAGTCAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGt             <  2:571411/140‑1 (MQ=255)
                                                                                                                                         ttCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAAATACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCAc      >  1:727101/1‑143 (MQ=255)
                                                                                                                                           ccGTGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACca    >  2:787895/1‑143 (MQ=255)
                                                                                                                                           ccGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACca    >  1:914550/1‑143 (MQ=255)
                                                                                                                                            cgcgGATCCCGAATACACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGCTTGACGCCGCCCGCGTCaaacatcct   <  2:40038/143‑8 (MQ=255)
                                                                                                                                            cgcgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTCACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcag   >  2:331321/1‑143 (MQ=255)
                                                                                                                                            cgcgGATCACGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGACCTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcag   >  1:383163/1‑143 (MQ=255)
                                                                                                                                             gcgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGCACTCCAGCAGGATGTTTGACCCGGGCTGCTTAAATCACcaga  >  1:132173/1‑142 (MQ=255)
                                                                                                                                             gcgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagc  >  1:872819/1‑143 (MQ=255)
                                                                                                                                             gcgGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGATAACAGCGCAGCACTCACGCAGGATGTTTGACGCGGGCGGCGTCAAGCACcagc  >  1:215663/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                             
TGTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGCGGGCGGCGTCAAGCACCAGC  >  NC_000913/4196088‑4196371

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