Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,017,114 T→C 6.3% G187G (GGT→GGC fliI → flagellum‑specific ATP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,017,1140TC 6.3% 197.6 / 2.2 63G187G (GGT→GGCfliIflagellum‑specific ATP synthase
Reads supporting (aligned to +/- strand):  ref base T (21/38);  new base C (2/2);  total (23/40)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGATTGC  >  NC_000913/2016972‑2017260
                                                                                                                                              |                                                                                                                                                  
gTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAgt                                                                                                                                              <  2:365461/149‑1 (MQ=255)
gTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAgt                                                                                                                                              <  2:359244/149‑1 (MQ=255)
gTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAgt                                                                                                                                              <  4:674484/149‑1 (MQ=255)
          tACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTgg                                                                                                                                    <  4:153056/149‑1 (MQ=255)
          tACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTgg                                                                                                                                    <  2:326031/149‑1 (MQ=255)
            ccccGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTgggg                                                                                                                                  <  1:391233/149‑1 (MQ=255)
             cccGCCATTTAACTCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGCAAAAGTGTGCTGCTGGGga                                                                                                                                 >  3:624551/1‑149 (MQ=255)
             cccGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGga                                                                                                                                 <  1:496324/149‑1 (MQ=255)
               cGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGgatg                                                                                                                               >  2:591428/1‑149 (MQ=255)
                    tttAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGc                                                                                                                          >  4:429001/1‑149 (MQ=255)
                      tAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCAc                                                                                                                        <  4:672893/149‑1 (MQ=255)
                          ccGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTa                                                                                                                    >  4:512295/1‑149 (MQ=255)
                           cGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTac                                                                                                                   <  3:429001/149‑1 (MQ=255)
                                       cacCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGAt                                                                                                       <  3:263435/149‑1 (MQ=255)
                                        acCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATg                                                                                                      >  1:69595/1‑149 (MQ=255)
                                        acCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATg                                                                                                      >  1:670097/1‑149 (MQ=255)
                                        acCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATg                                                                                                      >  1:300728/1‑149 (MQ=255)
                                        acCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATg                                                                                                      >  3:283057/1‑149 (MQ=255)
                                                 aCATGCGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCg                                                                                             >  4:492309/1‑149 (MQ=255)
                                                  cATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGt                                                                                            <  2:143645/149‑1 (MQ=255)
                                                   aTGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTg                                                                                           <  2:152726/149‑1 (MQ=255)
                                                             cacCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTg                                                                                 >  2:421827/1‑149 (MQ=255)
                                                              acCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTgg                                                                                <  4:464563/149‑1 (MQ=255)
                                                               ccGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGt                                                                               <  1:702973/149‑1 (MQ=255)
                                                                cGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTg                                                                              >  4:139575/1‑149 (MQ=255)
                                                                 ggCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGa                                                                             >  2:420350/1‑149 (MQ=255)
                                                                     tGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGt                                                                         >  3:68257/1‑149 (MQ=255)
                                                                          ccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTggg                                                                                                                                   >  1:236856/1‑86 (MQ=255)
                                                                          ccAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTggg                                                                                                                                   <  2:236856/86‑1 (MQ=255)
                                                                             atcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGa                                                                 >  4:185608/1‑149 (MQ=255)
                                                                              tcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGaa                                                                <  4:204975/149‑1 (MQ=255)
                                                                              tcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGaa                                                                <  2:87090/149‑1 (MQ=255)
                                                                              tcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGaa                                                                <  1:560123/149‑1 (MQ=255)
                                                                              tcaatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGCAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGaa                                                                <  4:624551/149‑1 (MQ=255)
                                                                              tcaatGCCCTGCGTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGaa                                                                <  2:463098/149‑1 (MQ=255)
                                                                                aatGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGt                                                              <  3:286964/149‑1 (MQ=255)
                                                                                   gCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaaa                                                           >  3:284017/1‑149 (MQ=255)
                                                                                    cccTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTaaaa                                                          <  2:543433/149‑1 (MQ=255)
                                                                                        gCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGAtt                                                      <  1:591428/149‑1 (MQ=255)
                                                                                         cTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGAttt                                                     <  3:492309/149‑1 (MQ=255)
                                                                                          ttACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGACGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGAtttt                                                    >  2:197904/1‑149 (MQ=255)
                                                                                           tACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTa                                                   <  2:96160/149‑1 (MQ=255)
                                                                                            aCCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTAt                                                  <  3:527020/149‑1 (MQ=255)
                                                                                             ccGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTAtt                                                 <  1:553280/149‑1 (MQ=255)
                                                                                             ccGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTAtt                                                 <  2:444523/149‑1 (MQ=255)
                                                                                                ttGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTgag                                              <  2:670097/149‑1 (MQ=255)
                                                                                                ttGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGCAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTgag                                              >  3:398289/1‑149 (MQ=255)
                                                                                                  gggCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGaa                                            <  2:197185/149‑1 (MQ=255)
                                                                                                        gggCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCt                                      >  3:301447/1‑149 (MQ=255)
                                                                                                         ggCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTc                                     <  2:130977/149‑1 (MQ=255)
                                                                                                              cGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGc                                <  3:478553/149‑1 (MQ=255)
                                                                                                              cGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGc                                >  1:583720/1‑149 (MQ=255)
                                                                                                                 aTGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGa                             >  2:271581/1‑149 (MQ=255)
                                                                                                                    gggCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAgg                          <  1:271581/149‑1 (MQ=255)
                                                                                                                    gggCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAgg                          <  3:543150/149‑1 (MQ=255)
                                                                                                                      gCTGTTTGGCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGGTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGc                        >  2:387953/1‑149 (MQ=255)
                                                                                                                             gCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACg                 >  3:524821/1‑149 (MQ=255)
                                                                                                                                   tCCGGCGTTGGCAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAgt           <  4:398289/149‑1 (MQ=255)
                                                                                                                                    ccGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAgtg          <  3:179376/149‑1 (MQ=255)
                                                                                                                                    ccGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAgtg          <  2:300728/149‑1 (MQ=255)
                                                                                                                                     cGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAgtgg         <  3:122656/149‑1 (MQ=255)
                                                                                                                                            ggTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGATTGc  <  4:524821/149‑1 (MQ=255)
                                                                                                                                            ggTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGATTGc  <  2:69595/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGCTGGACACCGGCGTGCGCCCAATCAATGCCCTGCTTACCGTTGGGCGTGGGCAGCGTATGGGGCTGTTTGCCGGGTCCGGCGTTGGTAAAAGTGTGCTGCTGGGGATGATGGCACGTTACACCCGCGCCGATGTCATTGTCGTGGGTTTGATTGGTGAACGTGGGCGCGAAGTAAAAGATTTTATTGAGAACATCCTCGGTGCCGAAGGGCGTGCACGCTCAGTGGTGATTGC  >  NC_000913/2016972‑2017260

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