Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,087,225 G→A 10.1% intergenic (+38/+880) PP_1865 → / ← aroH transposase/3‑deoxy‑7‑phosphoheptulonate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,087,2250GA10.1% 178.2 / 8.8 69intergenic (+38/+880)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base G (34/28);  new base A (5/2);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

CCACGAGGCCCCAAGGTGCTCAAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAT  >  NC_002947/2087077‑2087373
                                                                                                                                                       |                                                                                                                                                     
ccACGAGGCCCCAAGGTGCTCAAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTgg                                                                                                                                                       >  3:421522/1‑149 (MQ=40)
ccACGAGGCCCCAAGGTGCTCAAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTgg                                                                                                                                                       <  8:2956/149‑1 (MQ=40)
         cccAAGGTGCTCAAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAgg                                                                                                                                               >  3:385074/1‑148 (MQ=255)
         cccAAGGTGCTCAAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGa                                                                                                                                              <  5:440777/149‑1 (MQ=255)
            aaGGTTCTCAAACCGAACGCCTGGGTGCAAGGGACGGCTGTACATTGTCATTTTTCGACGGTTCGGGTCATCAAACCCGCCAAAAACAAAAGACCTTGAAAGGGATGGCCCTTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAAc                                                                                                                                            <  1:291582/149‑1 (MQ=17)
            aaGGTGCTCAAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAAc                                                                                                                                            >  8:410860/1‑148 (MQ=255)
                   tCAAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGCCACGACAAAAAggga                                                                                                                                             <  3:255471/140‑3 (MQ=14)
                   tCAAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGCCACGACAAAAAggga                                                                                                                                             >  4:255471/1‑138 (MQ=14)
                    cAAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCg                                                                                                                                   >  6:360123/1‑149 (MQ=255)
                      aaGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGTATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTc                                                                                                                                 <  8:302117/149‑1 (MQ=255)
                       aGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCg                                                                                                                                <  1:317344/149‑1 (MQ=255)
                         cGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGCCACGACAAAAAgg                                                                                                                                               <  6:65911/132‑1 (MQ=255)
                         cGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGCCACGACAAAAAgg                                                                                                                                               >  5:65911/1‑132 (MQ=255)
                               ccTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGAGGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCACAACCCGTCGGCTCCcat                                                                                                                        >  6:95732/1‑149 (MQ=255)
                                cTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCcat                                                                                                                        >  8:290962/1‑148 (MQ=255)
                                   ggTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGTCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTACCatcatc                                                                                                                    >  4:111586/1‑149 (MQ=255)
                                    gTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATcc                                                                                                                   <  7:410860/149‑1 (MQ=255)
                                    gTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATcc                                                                                                                   >  8:192550/1‑149 (MQ=255)
                                     tGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATAATCCg                                                                                                                  >  8:193307/1‑149 (MQ=255)
                                       cAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGt                                                                                                                                  <  6:7033/131‑1 (MQ=255)
                                       cAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGt                                                                                                                                  >  5:7033/1‑131 (MQ=255)
                                       cAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGc                                                                                                                <  8:207633/149‑1 (MQ=255)
                                        aaGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGcc                                                                                                               >  5:341485/1‑149 (MQ=255)
                                         cgggACGGCGGTCCATCCTCATGTGTCTACGGATCGGGTCATCAAAGCCTCCAAAACCAAAAGACCTTGAAAGGCACGGCCCTTACGGCGGGTCCCTTTTTGGCGTGGTAAAAACGAACCAAAAACCGTCGGCTCCCATCATCCGGccc                                                                                                              <  4:45792/148‑1 (MQ=18)
                                               ggCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGc                                                                                                        >  5:96465/1‑149 (MQ=255)
                                                    tACATGCTCATGTTTCGACGGATCGGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTgcg                                                                                                    >  1:5384/1‑149 (MQ=255)
                                                     aCATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCt                                                                                                  >  3:157968/1‑149 (MQ=255)
                                                       aTGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTcc                                                                                                >  5:369718/1‑149 (MQ=255)
                                                          cTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCggg                                                                                             <  6:341485/149‑1 (MQ=255)
                                                           tCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCCTCCGCTCCCATCAAACGGCCCCAACACTGCGCTcacggg                                                                                            >  2:216699/1‑144 (MQ=18)
                                                                tttCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTcccc                                                                                       >  3:185383/1‑149 (MQ=255)
                                                                tttCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTccc                                                                                        <  3:413210/148‑1 (MQ=255)
                                                                tttCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTccc                                                                                        <  3:419761/148‑1 (MQ=255)
                                                                     aCGGATCGGGTCATCAAAGCCGCCAAAGACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGc                                                                                   >  1:240072/1‑149 (MQ=255)
                                                                       ggATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTcc                                                                                >  3:127530/1‑149 (MQ=255)
                                                                              ggTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCtt                                                                          <  4:157968/149‑1 (MQ=255)
                                                                              ggTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCtt                                                                          >  3:208610/1‑149 (MQ=255)
                                                                                        gccgccAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGa                                                                <  7:290962/149‑1 (MQ=255)
                                                                                          cgacaaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGc                                                              <  5:360123/146‑1 (MQ=255)
                                                                                           gccAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCa                                                             >  2:384885/1‑149 (MQ=255)
                                                                                           gccAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCa                                                             >  5:106597/1‑149 (MQ=255)
                                                                                           gccAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCa                                                             >  5:112765/1‑149 (MQ=255)
                                                                                              aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCACGCACCCGCCAAGCGCCCCGGCGCAccc                                                          >  5:351813/1‑147 (MQ=17)
                                                                                               aaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCAcgcg                                                         <  4:385074/149‑1 (MQ=255)
                                                                                                aaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACgcgc                                                        <  4:127530/149‑1 (MQ=255)
                                                                                                       aaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAgg                                                  <  4:208610/149‑1 (MQ=255)
                                                                                                        aGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGc                                                <  7:264254/149‑1 (MQ=255)
                                                                                                           ccTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCGGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGa                                                                >  4:88414/1‑131 (MQ=255)
                                                                                                            cTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGcccc                                             >  2:243850/1‑148 (MQ=255)
                                                                                                                 aaGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCccatccat                                       <  5:238157/149‑1 (MQ=255)
                                                                                                                   gggATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATgg                                     <  6:96465/149‑1 (MQ=255)
                                                                                                                      aTGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGc                                  >  4:107832/1‑149 (MQ=255)
                                                                                                                      aTGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCACGCTCCGGCCTTGCTACCGGGAGCACCCGCTGCAGGCCCCAACCAAGGGGc                                  >  2:3496/1‑149 (MQ=14)
                                                                                                                          cccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATGGATGGGGCCtgc                              <  6:369718/149‑3 (MQ=255)
                                                                                                                          cccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCa                              <  2:5384/149‑1 (MQ=255)
                                                                                                                          cccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCCGCCATGCTACCGGGAGCACCCGCAGCCCGCCCCCTCCAAGGGGccg                                >  2:386059/1‑146 (MQ=14)
                                                                                                                          cccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCAGGCTGCGGCCGTGGTCCCCGGGGCCCGCCCGGCCGGCCCCCACCCAGGGGGcg                                 >  5:266985/1‑146 (MQ=11)
                                                                                                                           ccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATGGATGGGGCCtgca                             <  7:192550/149‑4 (MQ=255)
                                                                                                                           ccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAg                             >  3:43604/1‑149 (MQ=255)
                                                                                                                           ccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAg                             <  1:243850/149‑1 (MQ=255)
                                                                                                                           ccTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATACGGCCCCTAAGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCACCCGGCAGCACCCCCTCCAGGCCCCATCCATTGGGccgaac                             >  5:409083/1‑145 (MQ=14)
                                                                                                                            cTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCCgc                            >  5:313959/1‑149 (MQ=255)
                                                                                                                            cTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAgc                            <  1:173064/149‑1 (MQ=255)
                                                                                                                                aCGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAgcg                           <  3:335490/147‑1 (MQ=255)
                                                                                                                                 cggcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCt                         <  3:107832/148‑1 (MQ=255)
                                                                                                                                       gttccTTTTTGCCACGACAAAAAGGGACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACggg                  >  5:107770/4‑149 (MQ=17)
                                                                                                                                        ttccTTTTTGCCACGACAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGc                 >  2:281866/3‑149 (MQ=255)
                                                                                                                                         tccTTTTTGCCACGACAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGc                 >  1:144483/2‑148 (MQ=255)
                                                                                                                                         cccTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCACGCTCCGGCCTTTCTACCGGGAGCACGCGCAGCCGGCCCCATCCATGGGGCCTCAGCCCTTGACGGGc                 >  1:87784/1‑148 (MQ=17)
                                                                                                                                                 gTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGc         <  4:43604/148‑1 (MQ=255)
                                                                                                                                                    gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGt     >  2:189928/1‑149 (MQ=255)
                                                                                                                                                       gCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa  <  1:384885/149‑2 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CCACGAGGCCCCAAGGTGCTCAAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAT  >  NC_002947/2087077‑2087373

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: