Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,096,0900AG40.3% 117.3 / 20.0 72I83V (ATT→GTT) hisFimidazole glycerol phosphate synthase, catalytic subunit with HisH
Reads supporting (aligned to +/- strand):  ref base A (4/39);  new base G (29/0);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 3.57e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTG  >  NC_000913/2095945‑2096238
                                                                                                                                                 |                                                                                                                                                    
gCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGtt                                                                                                                                                    >  1:785435/1‑148 (MQ=255)
gCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTaa                                                                                                                                                  >  1:571256/1‑150 (MQ=255)
   ggCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTaaagtt                                                                                                                                              >  1:515334/1‑147 (MQ=255)
      aaaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTg                                                                                                                                           <  1:807760/151‑1 (MQ=255)
      aaaaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGtt                                                                                                                                                    >  1:746316/1‑142 (MQ=255)
        aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGa                                                                                                                                         <  1:1008441/151‑1 (MQ=255)
        aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGa                                                                                                                                         <  1:906791/151‑1 (MQ=255)
        aaCGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAATCTCTgg                                                                                                                                          >  1:112046/1‑150 (MQ=255)
         acgctacgctGAAGAAGGCGCTGACGAACTGGTGCTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGAGTTAGActcttg                                                                                                                                          >  2:893866/1‑147 (MQ=255)
          cgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTAGTCTCTgg                                                                                                                                          >  2:345884/1‑148 (MQ=255)
          cgctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAg                                                                                                                                       <  2:606442/151‑1 (MQ=255)
           gctacgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGTctct                                                                                                                                            >  2:45375/1‑145 (MQ=255)
               cgctGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGcc                                                                                                                                  <  1:78906/151‑1 (MQ=255)
                   gaagaaGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCcgc                                                                                                                                <  2:1190715/149‑1 (MQ=255)
                   gaagaaGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCcgc                                                                                                                                <  1:1084217/149‑1 (MQ=255)
                   gaagaaGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTGTTTGGGAGAAGcc                                                                                                                                  >  2:1000199/1‑147 (MQ=255)
                      gaagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGAAGATGCCGCGaaa                                                                                                                            >  2:349077/1‑150 (MQ=255)
                       aagGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAAt                                                                                                                          <  2:204143/151‑1 (MQ=255)
                        agGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAAtt                                                                                                                         <  2:575856/151‑1 (MQ=255)
                         ggCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTc                                                                                                                        <  1:171216/151‑1 (MQ=255)
                            gcTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGGAGAAGGCGCGAAAAtttttt                                                                                                                     >  2:894817/1‑151 (MQ=255)
                            gcTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCttt                                                                                                                     <  1:1205004/151‑1 (MQ=255)
                            gcTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCttt                                                                                                                     <  2:810157/151‑1 (MQ=255)
                             cTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTAGAATCTGGGAGATGCCGCGAAAAtttttt                                                                                                                     >  2:910942/1‑150 (MQ=255)
                              tGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTAATCTCTTGTAGAAGCAGCGGAAATTATTTcc                                                                                                                   >  1:319491/1‑151 (MQ=255)
                              tGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAATCTTTGGAAAAAGGCGCGCAAAATtttt                                                                                                                     >  1:197284/1‑149 (MQ=255)
                                aCGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCtt                                                                                                                 <  1:235989/151‑1 (MQ=255)
                                  gAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTg                                                                                                               <  1:154675/151‑1 (MQ=255)
                                  gAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTTAGGCTCTGGAAGATGCCGCGGAAAAtcttttctttt                                                                                                               >  2:921861/1‑150 (MQ=255)
                                   aaCTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGTAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTgg                                                                                                              <  2:377391/151‑1 (MQ=255)
                                   aaCTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTgg                                                                                                              <  2:769513/151‑1 (MQ=255)
                                       ggTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGGTTAGTCCTTGGTAGAAGCCCCGGAAAATTTTTCCTTTTGggcg                                                                                                          >  2:1133822/1‑151 (MQ=255)
                                       ggTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGgcgc                                                                                                           <  2:511203/150‑1 (MQ=255)
                                       ggTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGgcg                                                                                                            <  1:109127/149‑1 (MQ=255)
                                          gtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGa                                                                                                        <  1:1192182/150‑1 (MQ=255)
                                          gtTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGAGTAAGGCTCTTGAAGATTCCGCGAAAAAtcttttcttttgggggg                                                                                                         >  2:271599/1‑142 (MQ=255)
                                                cGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTAGAATCTGGGAGATTGCGCGAAAATTTTTTTCTTTTGGGCGGATaaaa                                                                                                   >  2:618110/1‑149 (MQ=255)
                                                cGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAAtt                                                                                                 <  2:677273/151‑1 (MQ=255)
                                                        ccGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCaa                                                                                         <  2:145659/151‑1 (MQ=255)
                                                         cGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAAACTTTTGGAGATGCCCCGGAAATTTTTTTCTTTTGGGCGGGTAAAATTTCCATcca                                                                                         >  2:833865/1‑148 (MQ=255)
                                                                aGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTg                                                                                 <  2:1024076/151‑1 (MQ=255)
                                                                      ggCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTgg                                                                           <  2:627874/151‑1 (MQ=255)
                                                                      ggCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTTTGGGAGAAGCCGCGAAAAATCTTTTCTTTTGCGCGGGTAAAATTTTCATCAACTCTCCTgcg                                                                               >  1:942843/1‑147 (MQ=255)
                                                                       gCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGGCTCTGGGAGATGCCGCGAAAAATCTTTTCTTTTGCGCGGGTAAAAATTTCCTCAAATCTCCTTCGCCgg                                                                           >  2:1183339/1‑150 (MQ=255)
                                                                         cGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTTCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTgg                                                                           >  1:625327/1‑148 (MQ=255)
                                                                            gtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAATATCCGGGAGATGCCGCGGAAATTTTTTTCTTTTGGGCGGGGAAAATTTTCATCAAAATCCCTGCGCttgtggacc                                                                     >  2:53615/1‑143 (MQ=255)
                                                                                       aaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAGTCTCTGGGAGATGCCGCGAAAATTTTTTCCTTTTGGGCGGGTAAAATTTCCATCAACTCTCCCGCGCTGGCAGACCCAACaattatt                                                          >  2:446443/1‑151 (MQ=255)
                                                                                            gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACt                                                       <  1:1071204/149‑1 (MQ=255)
                                                                                            gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCg                                                     <  1:45376/151‑1 (MQ=255)
                                                                                              gTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGcc                                                   <  1:695688/151‑1 (MQ=255)
                                                                                                aTCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTATGTCTCTGGGAGATGCCCCGTAAATTTTTTTCTTTTGCGCGGGTAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACCTTAATTACTCGcc                                                   >  2:254590/1‑149 (MQ=255)
                                                                                                 tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAATCTCTGGGAGATGCCGCGAAAATTTTTTCCTTTTGCGCGGGTAAAATTTCCATCAACTCTCCTTCGCTGGCAGACCCAACATTAATTAATCGCCTgg                                                >  2:640148/1‑151 (MQ=255)
                                                                                                    cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTg                                                 <  1:173373/147‑1 (MQ=255)
                                                                                                    cgcgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCg                                             <  2:892921/151‑1 (MQ=255)
                                                                                                            ggAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCttt                                     <  2:1039350/151‑1 (MQ=255)
                                                                                                            ggAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTTTCTGGAAGATGCCGCGAAAAATCTTTCCTTTTGCGCGGATAAAAATTCCATCAACTCTCCTTCGCTGGCAGACCCAACATTAATTTATCGCCCGGCCGATCGCttt                                     >  1:1034770/1‑151 (MQ=255)
                                                                                                             gAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTg                                    <  1:1000248/151‑1 (MQ=255)
                                                                                                              aaGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTGAGTCTGCGGGGGATTACGCGAAAATTCTTTCCTTTGGCGCGGGTAAAATTTCCATCAACACTCCTTGGCTCGCCGACCCCACATTCATTTCTTCCCTTGCCgcccccttt                                     >  2:706420/1‑141 (MQ=255)
                                                                                                               aGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGc                                  <  1:572835/151‑1 (MQ=255)
                                                                                                                  gATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTTAGTCTCTGGAAGAAGCCGCGGAAATTTTTTCCTTTGGGGCGGGTAAAAATTTCATCAACTCTCCTGTGGTGGCAGACCCCACATTAATTACTCGCCCGGCCGATctgttt                                     >  1:547013/1‑140 (MQ=255)
                                                                                                                   aTCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAATCTCTGGAAGAAGCCGCGGAAATTTTTTCCTTTTGCGCGGATAAAATTTCCCTCCACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTg                               >  1:313138/1‑150 (MQ=255)
                                                                                                                   aTCGACATTCCGTTTTGTGTGGCGGGAGGGGTTAATTCTCCGGAAGAAGCCGCGTAAATTCTTTCCTTTTGCGCGGGTAAAAATTACCACAAATCTCCCTCGCCGGCAGACCCAACATTTATTACTCCCCTGGCCGCTCCCCTTgggggg                               >  1:1261776/1‑146 (MQ=255)
                                                                                                                     cGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCa                             <  1:893910/150‑1 (MQ=255)
                                                                                                                     cGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAg                            <  1:1267053/151‑1 (MQ=255)
                                                                                                                         aTTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTa                        <  2:133556/151‑1 (MQ=255)
                                                                                                                            ccGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtg                     <  2:1195322/151‑1 (MQ=255)
                                                                                                                             cGTTTTGTGTGGCGGGTGGGGTTAAATCTCTTGAAGATGCCGCCGAAATTCTTTTCTTTTGCGCGGATAAAATTTTCATCAACTCTCCTGCGCTTGCGGACCCAACATTAATTTATCCCCTGGCCGATTGCTTTTGCCTGCAGTGTAtt                      >  2:963939/1‑149 (MQ=255)
                                                                                                                             cGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATtgt                    <  2:430470/151‑1 (MQ=255)
                                                                                                                                 ttGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGt                 <  2:947726/150‑1 (MQ=255)
                                                                                                                                  tgtgtgGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTc                <  2:716285/150‑1 (MQ=255)
                                                                                                                                     gtgGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGt             <  1:501924/150‑1 (MQ=255)
                                                                                                                                               ggATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTg  >  1:848890/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCTGGCAAAACGCTACGCTGAAGAAGGCGCTGACGAACTGGTGTTCTACGATATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCGCGAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCAACATTAATTACTCGCCTGGCCGATCGCTTTGGCGTGCAGTGTATTGTGGTCGGTATTGATACCTG  >  NC_000913/2095945‑2096238

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