Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,059,948 T→A 5.1% intergenic (+22/+16) asnU → / ← cbl tRNA‑Asn/ssuEADCB/tauABCD operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,059,9480TA 5.1% 283.9 / 4.6 79intergenic (+22/+16)asnU/cbltRNA‑Asn/ssuEADCB/tauABCD operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (45/30);  new base A (2/2);  total (47/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GGTTTGACAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCT  >  NC_000913/2059804‑2060091
                                                                                                                                                |                                                                                                                                               
ggTTTGACAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCGGGCttttcagg                                                                                                                                             <  4:573963/149‑5 (MQ=255)
   ttGACAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCgg                                                                                                                                          <  3:625043/149‑1 (MQ=255)
   ttGACAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCgg                                                                                                                                          <  3:29694/149‑1 (MQ=255)
     gACAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGa                                                                                                                                        >  3:450096/1‑149 (MQ=255)
      aCAAAAGATTTTTCACCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGAt                                                                                                                                       >  3:55860/1‑149 (MQ=255)
       cAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAATCCCGCTATAAAAGCGGGCtt                                                                                                                                      <  3:329498/149‑1 (MQ=21)
       cAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGAtt                                                                                                                                      >  4:135128/1‑149 (MQ=255)
          aaGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCATAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGAttttt                                                                                                                                   >  4:526164/1‑149 (MQ=255)
            gATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGc                                                                                                                                 >  2:317158/1‑149 (MQ=255)
                 ttCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTatat                                                                                                                            <  4:625196/149‑1 (MQ=255)
                      cGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGAt                                                                                                                       >  4:86059/1‑149 (MQ=255)
                      cGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGAt                                                                                                                       <  4:55860/149‑1 (MQ=255)
                       gTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATa                                                                                                                      >  1:211971/1‑149 (MQ=255)
                        ttAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGa                                                                                                                        <  1:147824/146‑1 (MQ=255)
                        ttAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGa                                                                                                                        >  2:147824/1‑146 (MQ=255)
                           aGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTtata                                                                                                                                              >  1:180428/1‑121 (MQ=255)
                           aGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTtata                                                                                                                                              <  2:180428/121‑1 (MQ=255)
                                 gCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATaa                                                                                                                     <  1:329632/140‑1 (MQ=255)
                                 gCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATaa                                                                                                                     >  2:329632/1‑140 (MQ=255)
                                   cTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCctc                                                                                                          >  3:282758/1‑149 (MQ=255)
                                   cTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCctc                                                                                                          >  1:494173/1‑149 (MQ=255)
                                      aacaacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTc                                                                                                       >  2:284974/1‑149 (MQ=255)
                                       acaacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTtata                                                                                                                                              >  2:631465/1‑109 (MQ=255)
                                       acaacGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTtata                                                                                                                                              <  1:631465/109‑1 (MQ=255)
                                               tCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTc                                                                                              <  3:422935/149‑1 (MQ=255)
                                               tCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTc                                                                                              >  3:206267/1‑149 (MQ=255)
                                                  tctGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCAt                                                                                           >  3:274630/1‑149 (MQ=255)
                                                    tGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCGGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCa                                                                                         >  1:798472/1‑149 (MQ=255)
                                                    tGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTtata                                                                                                                                              <  1:656571/96‑1 (MQ=255)
                                                    tGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTtata                                                                                                                                              >  2:656571/1‑96 (MQ=255)
                                                      tAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATaa                                                                                                                     <  3:353368/119‑1 (MQ=255)
                                                      tAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATaa                                                                                                                     >  4:353368/1‑119 (MQ=255)
                                                              tCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGGCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCtt                                                                               <  2:494173/149‑1 (MQ=255)
                                                                 gTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGat                                                                            <  3:135128/149‑1 (MQ=255)
                                                                    gAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATc                                                                         >  1:607056/1‑149 (MQ=255)
                                                                    gAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATc                                                                         <  1:124512/149‑1 (MQ=255)
                                                                       cggcggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTc                                                                      <  3:245869/149‑1 (MQ=255)
                                                                       cggcggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTc                                                                      >  4:245869/1‑149 (MQ=255)
                                                                          cggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAAc                                                                   <  3:248324/149‑1 (MQ=255)
                                                                           ggACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACt                                                                  >  4:564417/1‑149 (MQ=255)
                                                                            gACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTg                                                                 <  2:211971/149‑1 (MQ=255)
                                                                             aCTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGa                                                                <  3:86059/149‑1 (MQ=255)
                                                                                 ttAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGt                                                            >  4:38432/1‑149 (MQ=255)
                                                                                   aaTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTcc                                                          >  2:382119/1‑149 (MQ=255)
                                                                                    aTCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCt                                                         >  3:101181/1‑149 (MQ=255)
                                                                                     tCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGat                                                                            <  3:291185/129‑1 (MQ=255)
                                                                                     tCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGat                                                                            >  4:291185/1‑129 (MQ=255)
                                                                                     tCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTg                                                        >  4:202878/1‑149 (MQ=255)
                                                                                              cACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACaaa                                               <  3:697959/149‑1 (MQ=255)
                                                                                                cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGAttttt                                                                                                                                   <  3:192136/63‑1 (MQ=255)
                                                                                                cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGAttttt                                                                                                                                   >  4:192136/1‑63 (MQ=255)
                                                                                                cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGt                                             >  1:760875/1‑149 (MQ=255)
                                                                                                cTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGt                                             >  3:67067/1‑149 (MQ=255)
                                                                                                 tGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGtt                                            >  3:705859/1‑149 (MQ=255)
                                                                                                   gTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTcc                                          >  4:804332/1‑149 (MQ=255)
                                                                                                    ttCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCa                                         >  1:400288/1‑149 (MQ=255)
                                                                                                         gTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAg                                    >  4:505097/1‑149 (MQ=255)
                                                                                                          tCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAgc                                   >  2:691943/1‑149 (MQ=255)
                                                                                                               tCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAg                              >  1:780150/1‑149 (MQ=255)
                                                                                                                 agagGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTcc                                                                                             >  2:371858/1‑84 (MQ=255)
                                                                                                                 agagGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTcc                                                                                             <  1:371858/84‑1 (MQ=255)
                                                                                                                    ggagCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATa                         <  3:38432/149‑1 (MQ=255)
                                                                                                                      agCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGc                                                       <  4:601250/117‑1 (MQ=255)
                                                                                                                      agCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGc                                                       >  3:601250/1‑117 (MQ=255)
                                                                                                                       gCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAAc                                                                   >  4:522252/1‑104 (MQ=255)
                                                                                                                       gCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAAc                                                                   <  3:522252/104‑1 (MQ=255)
                                                                                                                          aaaTTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACg                   >  3:412436/1‑149 (MQ=255)
                                                                                                                          aaaTTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACg                   >  4:44235/1‑149 (MQ=255)
                                                                                                                           aaTTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGc                  <  2:798472/149‑1 (MQ=255)
                                                                                                                           aaTTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGc                  <  4:705859/149‑1 (MQ=255)
                                                                                                                            aTTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCt                 <  4:450096/149‑1 (MQ=255)
                                                                                                                               ccTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGaa              <  4:274630/149‑1 (MQ=255)
                                                                                                                                gcaaaaaTCCCGCTATAAAAGCGGGCTTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGaa              <  1:163930/145‑1 (MQ=255)
                                                                                                                                gcaaaaaTCCCGCTATAAAAGCGGGCTTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGaa              >  2:163930/4‑148 (MQ=255)
                                                                                                                                cTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATa                         <  4:217489/137‑1 (MQ=255)
                                                                                                                                cTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATa                         >  3:217489/1‑137 (MQ=255)
                                                                                                                                 tGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                           <  4:732663/102‑1 (MQ=255)
                                                                                                                                 tGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTc                                                           >  3:732663/1‑102 (MQ=255)
                                                                                                                                 tGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGt                       <  2:654474/138‑1 (MQ=255)
                                                                                                                                 tGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGt                       >  1:654474/1‑138 (MQ=255)
                                                                                                                                   aaaaGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTg          >  2:44141/1‑149 (MQ=255)
                                                                                                                                     aaGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTc        <  4:206267/149‑1 (MQ=255)
                                                                                                                                      atcccGCTATAAAAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTcc       >  1:409522/3‑149 (MQ=255)
                                                                                                                                         ccGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCg    >  1:371623/1‑149 (MQ=255)
                                                                                                                                           gCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCt  >  4:244334/1‑149 (MQ=255)
                                                                                                                                           gCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCt  >  2:466724/1‑149 (MQ=255)
                                                                                                                                               ttATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTg                                                                              >  3:770567/1‑69 (MQ=255)
                                                                                                                                               ttATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTg                                                                              <  4:770567/69‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GGTTTGACAAAAGATTTTTCGCCGTTAAGATGTGCCTCAACAACGATTCCTCTGTAGTTCAGTCGGTAGAACGGCGGACTGTTAATCCGTATGTCACTGGTTCGAGTCCAGTCAGAGGAGCCAAATTCCTGAAAAGCCCGCTTTTATAGCGGGATTTTTGCTATATCTGATAATCAATTTCCTCTTCACTGCTTTCCATCACCTGCCGCTTGATATCCTCAACTGACAGTCCTGCATTACAAAGTTCCAGAAAGCGCCAGACATAGTTACGCTGAAGTTGTCCTCGCT  >  NC_000913/2059804‑2060091

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