Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,016,880 G→C 5.1% Q109H (CAG→CAC) ‡ fliI → flagellum‑specific ATP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,016,8800GC 5.1% 265.1 / 5.6 78Q109H (CAG→CAC) ‡fliIflagellum‑specific ATP synthase
Reads supporting (aligned to +/- strand):  ref base G (39/35);  new base C (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

GAAACGCACGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGC  >  NC_000913/2016734‑2017028
                                                                                                                                                  |                                                                                                                                                    
gAACCGCACGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGACGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGtt                                                                                                                                                    >  6:46338/1‑149 (MQ=255)
 aaaCGCACGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTg                                                                                                                                                   <  3:175469/149‑1 (MQ=255)
   aCGCACGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGcc                                                                                                                                                 >  3:146965/1‑149 (MQ=255)
     gCACGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGc                                                                                                                                                  >  5:250028/1‑146 (MQ=255)
     gCACGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGc                                                                                                                                                  <  6:250028/146‑1 (MQ=255)
       aCGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGCCTGCAAAGCGGCAAGCAGTTGCCGCTc                                                                                                                                             >  8:192283/1‑149 (MQ=255)
          aaGTAGAAAGCGAAGGCGTTGGCTTTAACGGTCAACGGCGGGTTTTAATGCCGCTGGAGGAAGTCGAAGGGGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACAGTTCGGCAGAATGGCTGCAAAGCGGCAAGCAGTTGCCGCTCggg                                                                                                                                          >  8:404093/1‑148 (MQ=255)
               gAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCgg                                                                                                                                      <  4:33354/148‑1 (MQ=255)
               gAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGc                                                                                                                                     >  1:210218/1‑149 (MQ=255)
                aaaGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGGTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCg                                                                                                                                    >  4:191758/1‑149 (MQ=255)
                  aGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGtt                                                                                                                                  <  3:51367/149‑1 (MQ=255)
                     gAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGttatt                                                                                                                               <  5:88390/149‑1 (MQ=255)
                               gCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGtt                                                                                                                     >  7:405216/1‑149 (MQ=255)
                               gCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCAGGCGTTATTAGGTCGCGt                                                                                                                      >  7:410815/1‑148 (MQ=255)
                                 tttAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCt                                                                                                                   >  2:161770/1‑149 (MQ=255)
                                 tttAACGGTCAACGGCTGTTTTTAATGCCGATGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGACAAAAACATTTCGGCAGAAGTGCTCCAAAGCGGCAAGCAGTTGCCGCGCGGTCCGGCGTTATTAGGTCGAGGTCt                                                                                                                   >  8:174208/1‑149 (MQ=255)
                                    aaCGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAACTGCTTGCCGCtttgca                                                                                                                                         >  5:25263/1‑119 (MQ=255)
                                    aaCGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAACTGCTTGCCGCtttgca                                                                                                                                         <  6:25263/124‑6 (MQ=255)
                                    aaCGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAACTGCTTGCCGCtttgca                                                                                                                                         <  6:31163/124‑6 (MQ=255)
                                    aaCGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAACTGCTTGCCGCtttgca                                                                                                                                         >  5:31163/1‑119 (MQ=255)
                                     aCGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGTCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGa                                                                                                                <  2:210218/148‑1 (MQ=255)
                                     aCGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGa                                                                                                                <  3:131541/148‑1 (MQ=255)
                                     aCGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGa                                                                                                                <  7:241840/148‑1 (MQ=255)
                                      cGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTTTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACg                                                                                                              <  7:140365/149‑1 (MQ=255)
                                             cGGCTGTTTTTAATGCCGCTGGAGTAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTc                                                                                                                          <  4:36342/130‑1 (MQ=255)
                                             cGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTc                                                                                                                          >  3:36342/1‑130 (MQ=255)
                                                cTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCGCGGGCAGGCGGTATTAGGTCGCGTTCTGGACGGCAGCGGTa                                                                                                     >  6:257901/1‑148 (MQ=255)
                                                   tttttAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGAGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAAcc                                                                                                 >  7:149050/1‑149 (MQ=255)
                                                         aTGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGa                                                                                           <  5:250106/149‑1 (MQ=255)
                                                                  gaggaAGTCGCAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGGCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGcc                                                                                  <  7:404093/149‑1 (MQ=255)
                                                                  gaggaAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGCCGGCAGCGGTAAACCGCTCGa                                                                                           >  8:363797/1‑140 (MQ=255)
                                                                  gaggaAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGa                                                                                           <  7:363797/140‑1 (MQ=255)
                                                                   aggaAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGccc                                                                                 >  1:345366/1‑149 (MQ=255)
                                                                    ggaAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCt                                                                                >  5:111286/1‑149 (MQ=255)
                                                                    ggaAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCt                                                                                >  5:162872/1‑149 (MQ=255)
                                                                        gTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCa                                                                                                           >  7:212902/1‑118 (MQ=255)
                                                                        gTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCa                                                                                                           <  8:212902/118‑1 (MQ=255)
                                                                         tCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGAt                                                                                          <  1:248614/134‑1 (MQ=255)
                                                                         tCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGAt                                                                                          >  2:248614/1‑134 (MQ=255)
                                                                         tCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTccccc                                                                           >  6:222279/1‑149 (MQ=255)
                                                                          cGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTccccc                                                                           >  5:94740/1‑148 (MQ=255)
                                                                           gAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGa                                                                         >  6:359490/1‑149 (MQ=255)
                                                                           gAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCACCCCCGa                                                                         >  5:79191/1‑149 (MQ=255)
                                                                            aaGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACTGCAGCGGTAAACCGCTCGATGGCCTGCCCACCCCCGAt                                                                        <  6:79191/149‑1 (MQ=255)
                                                                              ggTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGGTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATac                                                                      >  2:58905/1‑149 (MQ=255)
                                                                                 gtCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATacgac                                                                   >  7:292162/1‑149 (MQ=255)
                                                                                    cTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACgg                                                                 >  5:173719/1‑148 (MQ=255)
                                                                                    cTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACgg                                                                 >  2:378647/1‑148 (MQ=255)
                                                                                       cccGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAAc                                                             >  6:287648/1‑149 (MQ=255)
                                                                                         cGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCg                                                           <  2:275410/149‑1 (MQ=255)
                                                                                          ggCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCgg                                                          >  7:340571/1‑149 (MQ=255)
                                                                                                   gtTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGATCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGAt                                                 >  7:410287/1‑149 (MQ=255)
                                                                                                        tGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTAccc                                            <  1:378647/149‑1 (MQ=255)
                                                                                                          ccAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTAcccc                                           <  6:94740/148‑1 (MQ=255)
                                                                                                          ccAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCg                                          >  3:117553/1‑149 (MQ=255)
                                                                                                           cAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGc                                         <  5:222279/149‑1 (MQ=255)
                                                                                                           cAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGc                                         <  1:161770/149‑1 (MQ=255)
                                                                                                            aaaaaCATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGc                                         >  1:281313/1‑148 (MQ=255)
                                                                                                            aaaaaCATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGc                                         >  1:240205/1‑148 (MQ=255)
                                                                                                              aaaCATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCAt                                      <  5:233961/149‑1 (MQ=255)
                                                                                                                 cATTTCGGCTGAAGGGCTGCAACGCGGCAAGCAGTTTCCGCTCGGTCCGGAGTGTTTAGGTCGCGTTCTGGACGGAAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGAAGGAAACCGGTGCGCTGATTACCCCGCCATTTa                                   <  8:410815/149‑1 (MQ=255)
                                                                                                                 cATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCTGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCAATTa                                   <  8:405216/149‑1 (MQ=255)
                                                                                                                  aTTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTaa                                  <  8:340571/149‑1 (MQ=255)
                                                                                                                   tttCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAAc                                 <  6:162872/149‑1 (MQ=255)
                                                                                                                        gCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGtt                            <  5:359490/149‑1 (MQ=255)
                                                                                                                        gCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGtt                            <  5:287648/149‑1 (MQ=255)
                                                                                                                        gCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCGCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGtt                            >  6:167748/1‑149 (MQ=255)
                                                                                                                         cAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTg                           <  6:111286/149‑1 (MQ=255)
                                                                                                                          aGAAGGGCTGCAAAGCGGCAAGCAGGTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGc                          <  4:146965/149‑1 (MQ=255)
                                                                                                                            aaGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGTGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCaa                        >  3:240711/1‑149 (MQ=255)
                                                                                                                             aGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAAc                       <  3:191758/149‑1 (MQ=255)
                                                                                                                                  tGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGTGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTaca                  >  2:241962/1‑149 (MQ=255)
                                                                                                                                      aaGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGa              <  2:345366/149‑1 (MQ=255)
                                                                                                                                        gCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGGAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAACCCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGAtt            >  6:47729/1‑149 (MQ=255)
                                                                                                                                        gCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGGCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGAtt            <  5:46338/149‑1 (MQ=255)
                                                                                                                                           gCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGaa         <  4:117553/149‑1 (MQ=255)
                                                                                                                                            cAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCATCCCCCGAGCCGACGGCAACCGGGGCGCTGATGACACCGCCATTTAACCCGTGGCAACGTACACCGATTGaa         >  7:172214/1‑148 (MQ=255)
                                                                                                                                              aGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACAt      <  2:240205/149‑1 (MQ=255)
                                                                                                                                                  gTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTATGGACTGCAGCGGTAAACCGCTCGAAGGCCTGCCGTCCCCAGATACGACGGAAACCGGTGCGCTGATTTCCCCGCCATTTAACCCGTTGCAACGTCCACCGATTGAACATGTGc  >  5:328425/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAAACGCACGAAGTAGAAAGCGAAGTCGTTGGCTTTAACGGTCAACGGCTGTTTTTAATGCCGCTGGAGGAAGTCGAAGGTGTCCTGCCCGGCGCGCGTGTTTATGCCAAAAACATTTCGGCAGAAGGGCTGCAAAGCGGCAAGCAGTTGCCGCTCGGTCCGGCGTTATTAGGTCGCGTTCTGGACGGCAGCGGTAAACCGCTCGATGGCCTGCCCTCCCCCGATACGACGGAAACCGGTGCGCTGATTACCCCGCCATTTAACCCGTTGCAACGTACACCGATTGAACATGTGC  >  NC_000913/2016734‑2017028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: