Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 85,081 T→G 11.9% intergenic (+59/+71) aroE → / ← PP_0075 shikimate 5‑dehydrogenase/choline sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294785,0810TG11.9% 148.4 / 5.7 65intergenic (+59/+71)aroE/PP_0075shikimate 5‑dehydrogenase/choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base T (29/27);  new base G (2/6);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTG  >  NC_002947/84933‑85224
                                                                                                                                                      |                                                                                                                                                
gcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgt                                                                                                                                                  <  3:129905/149‑1 (MQ=255)
    cTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTa                                                                                                                                              <  2:129857/149‑1 (MQ=255)
         gCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGcc                                                                                                                                          >  7:118724/1‑148 (MQ=255)
            ttcttTATCTGGCGTGGCGTGAGGACGGATACGGCGACGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATc                                                                                                                                      >  5:88103/1‑149 (MQ=255)
             tcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCCTTCCCACAAAGTACGTCGTACGCCGATcc                                                                                                                                     >  3:275465/1‑149 (MQ=255)
              cttTATCTGTCCTGGCGTGAGGACGGATACGGCGCAGGTGTTGGCTTAACTGCGCCGGAAACTCGCTCTGTGATTAGCTGGACGGGGCCGCTTGGCTGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATc                                                                                                                                      <  1:347342/147‑1 (MQ=37)
                ttATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTAGTACGCCGATCcact                                                                                                                                  >  5:40492/1‑149 (MQ=255)
                 tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCctct                                                                                                                                  <  8:118724/148‑1 (MQ=255)
                  aTCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATAGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                >  7:224894/1‑149 (MQ=255)
                   tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                <  3:109705/148‑1 (MQ=255)
                   tCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                <  6:84710/148‑1 (MQ=255)
                    cTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTGCCACAAAGTACGTCGTACGCCGCTCCACTGTgg                                                                                                                              >  1:180265/1‑149 (MQ=255)
                     tggcgtggcgTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTAGTACGCCGATACACTGTgg                                                                                                                              >  4:285166/1‑148 (MQ=255)
                        cgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGc                                                                                                                          <  7:159158/149‑1 (MQ=255)
                          tggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACgg                                                                                                                        <  6:331809/149‑1 (MQ=21)
                          tggcgtgAGGCCGGATACGGCGCAGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCGTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACgg                                                                                                                        <  8:282473/149‑1 (MQ=21)
                           ggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGc                                                                                                                       >  1:330804/1‑149 (MQ=255)
                           ggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTAGCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGc                                                                                                                       <  5:264148/149‑1 (MQ=255)
                            gcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCATTCCGACAAAGTACGTCGAACGCCGATCCTCGGTGGGAGCgg                                                                                                                        >  4:336079/1‑147 (MQ=255)
                              gtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                <  2:129907/137‑1 (MQ=255)
                              gtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                >  1:129907/1‑137 (MQ=255)
                              gtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGc                                                                                                                          <  2:55197/143‑1 (MQ=255)
                              gtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGc                                                                                                                          >  1:55197/1‑143 (MQ=255)
                              gtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCtt                                                                                                                    <  7:285670/149‑1 (MQ=255)
                              gtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCtt                                                                                                                    <  8:224894/149‑1 (MQ=255)
                              gtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGGCGTACGCCGATCCTCTGTGGGAGCGGCCtt                                                                                                                    <  3:259710/149‑1 (MQ=255)
                                gAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCCATCCTCTGTGCGAGCGGCCTAGc                                                                                                                 >  7:85084/1‑149 (MQ=255)
                                  ggCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtcg                                                                                                                                                >  6:389440/1‑117 (MQ=255)
                                  ggCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtcg                                                                                                                                                <  5:389440/117‑1 (MQ=255)
                                            ggCGCCGGTGTTGGGTGAACTGCGCGGGCAACTGGGTGGGGGCTGAGCTGGCGGGGGCGGCTTGGGGGCCCTTTCGCGACGCAAGGCCGTTCCCACAAAGTACGTGGTACGCCGATCCTCTGTGGGTGCGGCCTTGCGTCGCGATggg                                                                                                      <  4:275465/148‑1 (MQ=37)
                                                        ggCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCCTTAACACAAAATAACGAGTACGCCGCAACACTGTGGGACCGGCCATGCGTCGCGATTGGCTGCAAACCCGAc                                                                                         >  6:245530/1‑149 (MQ=255)
                                                         gCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGAcc                                                                                        >  2:23535/1‑149 (MQ=255)
                                                          cTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGc                                                                                                                          <  3:44414/115‑1 (MQ=255)
                                                          cTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGc                                                                                                                          >  4:44414/1‑115 (MQ=255)
                                                           ttaccTGGGCGGGAAACTCGCGGGGGGCTGCGGGGGCGGGGGCCGGTGGGGGGCCCTATCGCGACCAATGGCCGTTCCCACAACTACGGTCGTACGCCTTCCCTGTGTGGGCGGGGCGTTGGGTCGCGATGGGCTGCAAAGCAGACCCa                                                                                      <  2:180265/145‑1 (MQ=255)
                                                             aaCTGCGCCGGCAACTCGCTGGGGGCCGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAg                                                                                              <  4:64161/139‑1 (MQ=255)
                                                             aaCTGCGCCGGCAACTCGCTCGGGGCCGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAgctggg                                                                                      >  3:64161/1‑142 (MQ=255)
                                                                  cgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATc                                                                                 <  1:23535/147‑1 (MQ=255)
                                                                     cGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCTTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAAc                                                                            <  1:51444/149‑1 (MQ=255)
                                                                         aaCTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCTATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCa                                                                        <  3:285166/149‑1 (MQ=255)
                                                                         aaCTCGCGCGGGGCTGCGCGGGCCGGGGCCGCTTGGCGGCCCTAGCGCGACCCAAGGCCGTTCCCACAAAGTACGGCGTACGCCTATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCa                                                                        <  1:155058/149‑1 (MQ=38)
                                                                             cgctggGGGCTGAGCTGGCCGGGGCGGCTTGGCGGCCATATCGGGACGCAAGGCCGTTCCCAAAAAGTACGTCTTACGGCGATCCTGTGTGGGAGGGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTc                                                                     <  8:85084/148‑1 (MQ=38)
                                                                               gtcgGGGCTGCGGGGCCGGGGGCCGTTTGGGGGCCCTAGCGCGACCCATGGCCGTCCCCACAAAGTACGGCGTCCGCCTATCCTCTGGGGGCGGGGCCTTGGGCCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTc                                                                  <  6:40492/148‑1 (MQ=255)
                                                                               ctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTc                                                                  <  3:281082/149‑1 (MQ=255)
                                                                                   gggCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACTCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGTAGGGGCGTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAAc                                                              <  6:88103/149‑1 (MQ=255)
                                                                                       tGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGa                                                          >  5:194854/1‑149 (MQ=255)
                                                                                        gggctggctgGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGAt                                                         >  6:301308/3‑149 (MQ=255)
                                                                                           cTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCgg                                                      >  1:76918/1‑149 (MQ=255)
                                                                                              gcgggggCCGCTTGGCGGCCCTAGCGCCCCGCAAGGCCGTTCCCACAACTTACGTCGTACGCCGATCCTCTGTGGCAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                   <  2:330804/146‑1 (MQ=38)
                                                                                                        cTTGGCGGCCCTATCGCGACGCAAGGCCATTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAcc                                         <  3:336079/149‑1 (MQ=255)
                                                                                                               gcCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCa                                  <  2:76918/149‑1 (MQ=255)
                                                                                                                 ccTATCGCGACGCAAGGCCGTTACCACAAAGTACGTCGTACCCCGATACTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTGCTCAAACCCGATAGGGCACAACTCCAccccgtaccgc                                >  8:5431/1‑142 (MQ=38)
                                                                                                                    aTCGCGGCGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAAc                                               >  1:224634/1‑131 (MQ=255)
                                                                                                                    aTCGCGACGCAAGGCCGTTCCCACAAAGTACGGCGTACGCCTATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAAc                                               <  2:224634/131‑1 (MQ=255)
                                                                                                                    aTCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTc                             >  8:29258/1‑149 (MQ=21)
                                                                                                                     tCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTc                             <  7:281902/148‑1 (MQ=255)
                                                                                                                     tCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCt                            >  1:310656/1‑149 (MQ=255)
                                                                                                                       gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                     >  8:121003/1‑122 (MQ=255)
                                                                                                                       gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTg                           <  6:194854/148‑1 (MQ=255)
                                                                                                                       gcgACGCAAGGCCGTTCCCACAAAGTACGCCGTACGCCGATCCTCTGTGGCAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTg                           <  5:245530/148‑1 (MQ=255)
                                                                                                                             cAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGGCGTCGCGATGGGCTGCAAAGCAGCCCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCCCCCCTTCCAGCTTCTGCAActc                     >  1:112470/1‑149 (MQ=255)
                                                                                                                              aaGGCCGTTCACACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCACAATCCCAACCTCATTCCTCAAACCGGATCGCGCACAACTCCACCCCTTCAAGCTTCTGCAActcct                   >  7:230042/1‑149 (MQ=255)
                                                                                                                                    gTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCa             >  3:289569/1‑149 (MQ=255)
                                                                                                                                        ccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCg                                                       >  6:282553/1‑103 (MQ=255)
                                                                                                                                         cacaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc        >  3:407778/1‑149 (MQ=255)
                                                                                                                                         cacaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc        >  2:78209/1‑149 (MQ=255)
                                                                                                                                         cacaAAGTACGTCGTACGCCGAGCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTGCACCCCTTCCAGCTTCTGCAACTCCTCGCTCACCTGc        >  7:260958/1‑149 (MQ=255)
                                                                                                                                               gTACGTCGTACGCCTTTCCTGTGTTGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTg  <  1:220480/149‑1 (MQ=255)
                                                                                                                                                      tcgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                          >  1:305632/1‑120 (MQ=255)
                                                                                                                                                      tcgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                          <  2:305632/120‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTG  >  NC_002947/84933‑85224

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