Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 397,277 C→G 5.7% I213M (ATC→ATG)  sbmA → peptide antibiotic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913397,2770CG 5.7% 316.8 / 9.1 88I213M (ATC→ATG) sbmApeptide antibiotic transporter
Reads supporting (aligned to +/- strand):  ref base C (45/38);  new base G (3/2);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

GTACAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTAGGG  >  NC_000913/397138‑397424
                                                                                                                                           |                                                                                                                                                   
gTACAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTg                                                                                                                                            >  1:204474/1‑149 (MQ=255)
   cAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATc                                                                                                                                         >  4:160914/1‑149 (MQ=255)
    aGCGATGAACGAATATTACATCGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCg                                                                                                                                        <  4:584874/149‑1 (MQ=255)
     gCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGc                                                                                                                                       <  2:204474/149‑1 (MQ=255)
       gATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCt                                                                                                                                     >  3:202832/1‑149 (MQ=255)
          gAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTcc                                                                                                                                  >  3:52398/1‑149 (MQ=255)
            aCGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTg                                                                                                                                >  2:11603/1‑149 (MQ=255)
              gAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGACGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGcc                                                                                                                              <  3:71466/149‑1 (MQ=255)
                   ttACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGc                                                                                                                         <  3:144925/149‑1 (MQ=255)
                   ttACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGc                                                                                                                         <  2:538562/149‑1 (MQ=255)
                      cATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTgg                                                                                                                      <  1:11603/149‑1 (MQ=255)
                         ggCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGAc                                                                                                                                                <  3:439141/120‑1 (MQ=255)
                         ggCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGAc                                                                                                                                                >  4:439141/1‑120 (MQ=255)
                         ggCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTaa                                                                                                                   <  1:104087/149‑1 (MQ=255)
                          gCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAAc                                                                                                                  <  3:343581/149‑1 (MQ=255)
                                tGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCtctc                                                                                                            >  2:203612/1‑149 (MQ=255)
                                       aacTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCAt                                                                                                     >  1:448076/1‑149 (MQ=255)
                                       aacTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCAt                                                                                                     >  2:13437/1‑149 (MQ=255)
                                         cTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTcc                                                                                                           >  1:326676/1‑141 (MQ=255)
                                         cTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTcc                                                                                                           <  2:326676/141‑1 (MQ=255)
                                          tGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCAtgtg                                                                                                  >  1:67891/1‑149 (MQ=255)
                                          tGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCAtgtg                                                                                                  >  4:610382/1‑149 (MQ=255)
                                           gCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGc                                                                                                 <  3:175864/149‑1 (MQ=255)
                                           gCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGc                                                                                                 <  3:568967/149‑1 (MQ=255)
                                             gtgatatCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCg                                                                                               >  4:102301/4‑149 (MQ=255)
                                              tCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCgg                                                                                              >  4:289007/1‑149 (MQ=255)
                                               cATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGa                                                                                             >  1:15519/1‑149 (MQ=255)
                                                    cGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGc                                                                                        <  2:448347/149‑1 (MQ=255)
                                                      aaGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCg                                                                                      >  3:45812/1‑149 (MQ=255)
                                                      aaGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCg                                                                                      >  1:124050/1‑149 (MQ=255)
                                                           gCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTAt                                                                                 >  2:570406/1‑149 (MQ=255)
                                                            ccGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATc                                                                                >  1:545311/1‑149 (MQ=255)
                                                               cacaGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCggg                                                                             <  1:203612/149‑1 (MQ=255)
                                                                          cAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTa                                                                  <  1:1011/149‑1 (MQ=255)
                                                                            ggAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATg                                                                <  1:457496/149‑1 (MQ=255)
                                                                             gAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATgg                                                               >  2:494075/1‑149 (MQ=255)
                                                                              aaGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGt                                                              >  2:199718/1‑149 (MQ=255)
                                                                                  cacCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTgg                                                          <  3:102301/149‑1 (MQ=255)
                                                                                     cATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGa                                                       <  1:570406/149‑1 (MQ=255)
                                                                                      aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGAt                                                      >  3:394364/1‑149 (MQ=255)
                                                                                      aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGAt                                                      >  1:475542/1‑149 (MQ=255)
                                                                                      aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGAt                                                      <  4:45812/149‑1 (MQ=255)
                                                                                      aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGAt                                                      >  2:148083/1‑149 (MQ=255)
                                                                                      aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGAt                                                      <  2:124050/149‑1 (MQ=255)
                                                                                          gTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGcc                                                  <  3:334042/149‑1 (MQ=255)
                                                                                            tttGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGc                                                >  4:470954/1‑149 (MQ=255)
                                                                                            tttGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGc                                                >  4:514550/1‑149 (MQ=255)
                                                                                              tGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCaa                                              <  3:242390/149‑1 (MQ=255)
                                                                                               gCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAAt                                             <  4:394364/149‑1 (MQ=255)
                                                                                               gCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAAt                                             <  4:52398/149‑1 (MQ=255)
                                                                                                 ttCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTa                                                                  <  1:243960/126‑1 (MQ=255)
                                                                                                 ttCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTa                                                                  >  2:243960/1‑126 (MQ=255)
                                                                                                 ttCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCg                                           >  3:55694/1‑149 (MQ=255)
                                                                                                 ttCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCg                                           >  3:488929/1‑149 (MQ=255)
                                                                                                 ttCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTACGGTCTGGTGATTGCCGCAATCg                                           >  3:351947/1‑149 (MQ=255)
                                                                                                  tCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGt                                          >  4:181609/1‑149 (MQ=255)
                                                                                                   cAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTc                                         <  3:610382/149‑1 (MQ=255)
                                                                                                   cAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTc                                         <  3:181609/149‑1 (MQ=255)
                                                                                                    aaCGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCt                                        <  3:289007/149‑1 (MQ=255)
                                                                                                       gCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCCGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGt                                     >  4:438711/1‑149 (MQ=255)
                                                                                                          ggAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGc                                  <  4:249712/149‑1 (MQ=255)
                                                                                                          ggAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGc                                  >  3:269600/1‑149 (MQ=255)
                                                                                                                  tGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATgggg                          <  1:494075/149‑1 (MQ=255)
                                                                                                                   gggCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGa                         >  3:312685/1‑149 (MQ=255)
                                                                                                                     gCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGAcc                       >  2:559079/1‑149 (MQ=255)
                                                                                                                     gCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGAcc                       >  4:200353/1‑149 (MQ=255)
                                                                                                                         cAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCCGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGAt                   <  3:438711/149‑1 (MQ=255)
                                                                                                                          aGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCTATCGTCTGGTCGCTGATGGGGACCGGAtt                  >  3:591040/1‑149 (MQ=255)
                                                                                                                           gTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATc                                            <  3:577461/122‑1 (MQ=255)
                                                                                                                           gTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATc                                            >  4:577461/1‑122 (MQ=255)
                                                                                                                           gTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGCGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCCGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTg                 >  2:361301/1‑149 (MQ=255)
                                                                                                                               gaTCAACGTCATGATGGCGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTa                                                                  <  2:346477/95‑1 (MQ=255)
                                                                                                                               gaTCAACGTCATGATGGCGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTa                                                                  >  1:346477/2‑96 (MQ=255)
                                                                                                                               gaTCAACGTCATGATGGCGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTgg             >  2:61525/2‑149 (MQ=255)
                                                                                                                                aTCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGc            >  3:168227/1‑149 (MQ=255)
                                                                                                                                 tCACCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCTATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCa           <  4:591040/149‑1 (MQ=255)
                                                                                                                                 tCAACGTCATGATGGCGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATggg                           <  2:538756/133‑1 (MQ=255)
                                                                                                                                 tCAACGTCATGATGGCGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATggg                           >  1:538756/1‑133 (MQ=255)
                                                                                                                                 tCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCa           <  3:160914/149‑1 (MQ=255)
                                                                                                                                   aaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGa                               <  3:617218/127‑1 (MQ=255)
                                                                                                                                   aaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGa                               >  4:617218/1‑127 (MQ=255)
                                                                                                                                   aaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGt         >  1:32851/1‑149 (MQ=255)
                                                                                                                                   aaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGt         >  3:331148/1‑149 (MQ=255)
                                                                                                                                    aCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTg        <  4:353589/149‑1 (MQ=255)
                                                                                                                                    aCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTg        >  3:84173/1‑149 (MQ=255)
                                                                                                                                      gCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGt      <  2:545311/149‑1 (MQ=255)
                                                                                                                                       ccATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTa     >  4:229309/1‑149 (MQ=255)
                                                                                                                                       ccATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTa     <  3:522849/149‑1 (MQ=255)
                                                                                                                                        catcatGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTAg    <  1:13437/149‑1 (MQ=255)
                                                                                                                                         atcatGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTAgg   <  2:67891/149‑1 (MQ=255)
                                                                                                                                          tcatGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTAggg  >  4:486973/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GTACAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTAGGG  >  NC_000913/397138‑397424

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