Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 397,273 C→T 5.8% A212V (GCC→GTC)  sbmA → peptide antibiotic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913397,2730CT 5.8% 290.5 / 5.6 86A212V (GCC→GTC) sbmApeptide antibiotic transporter
Reads supporting (aligned to +/- strand):  ref base C (44/37);  new base T (3/2);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

GTACAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTA  >  NC_000913/397138‑397421
                                                                                                                                       |                                                                                                                                                    
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                                                                                                                                             >  4:439141/1‑120 (MQ=255)
                         ggCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGAc                                                                                                                                             <  3:439141/120‑1 (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:568967/149‑1 (MQ=255)
                                           gCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGc                                                                                              <  3:175864/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                                                                                   >  1:124050/1‑149 (MQ=255)
                                                      aaGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCg                                                                                   >  3:45812/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                                                   <  4:45812/149‑1 (MQ=255)
                                                                                      aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGAt                                                   >  2:148083/1‑149 (MQ=255)
                                                                                      aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGAt                                                   <  2:124050/149‑1 (MQ=255)
                                                                                      aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGAt                                                   >  1:475542/1‑149 (MQ=255)
                                                                                      aTGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGAt                                                   >  3:394364/1‑149 (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:488929/1‑149 (MQ=255)
                                                                                                 ttCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCg                                        >  3:55694/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                               >  3:269600/1‑149 (MQ=255)
                                                                                                          ggAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGc                               <  4:249712/149‑1 (MQ=255)
                                                                                                                  tGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATgggg                       <  1:494075/149‑1 (MQ=255)
                                                                                                                   gggCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGa                      >  3:312685/1‑149 (MQ=255)
                                                                                                                     gCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGAcc                    >  4:200353/1‑149 (MQ=255)
                                                                                                                     gCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGAcc                    >  2:559079/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                        >  1:538756/1‑133 (MQ=255)
                                                                                                                                 tCAACGTCATGATGGCGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATggg                        <  2:538756/133‑1 (MQ=255)
                                                                                                                                 tCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCa        <  3:160914/149‑1 (MQ=255)
                                                                                                                                   aaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGa                            <  3:617218/127‑1 (MQ=255)
                                                                                                                                   aaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGa                            >  4:617218/1‑127 (MQ=255)
                                                                                                                                   aaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGt      >  3:331148/1‑149 (MQ=255)
                                                                                                                                   aaCGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGt      >  1:32851/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)
                                                                                                                                       |                                                                                                                                                    
GTACAGCGATGAACGAATATTACATGGCGAACTGGCAACAACTGCGTCATATCGAAGGGGCCGCACAGCGTGTGCAGGAAGACACCATGCGTTTTGCTTCAACGCTGGAGAATATGGGCGTCAGTTTTATCAACGCCATCATGACGTTGATCGCCTTCCTGCCGGTGCTGGTAACGCTCTCCGCGCATGTGCCGGAGCTGCCGATTATCGGGCACATTCCGTATGGTCTGGTGATTGCCGCAATCGTCTGGTCGCTGATGGGGACCGGATTGCTGGCAGTGGTA  >  NC_000913/397138‑397421

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: ‑