Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,0860AC41.7% 117.2 / 22.2 72Y580D (TAC→GAC) ggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (37/5);  new base C (0/30);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.98e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGC  >  NC_000913/3584938‑3585235
                                                                                                                                                     |                                                                                                                                                     
tattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcc                                                                                                                                                     >  2:967880/1‑151 (MQ=255)
   tGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccg                                                                                                                                                  >  1:89477/1‑151 (MQ=255)
    gATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcccc                                                                                                                                                   <  2:859045/149‑1 (MQ=255)
     aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgcc                                                                                                                                                <  1:854980/151‑1 (MQ=255)
              gtggtgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCTGTTAAatca                                                                                                                                       <  2:271097/151‑1 (MQ=255)
               tggtgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAatcat                                                                                                                                      >  1:769572/1‑151 (MQ=255)
                      attcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCg                                                                                                                                >  2:1078469/1‑150 (MQ=255)
                        tcatccTCTGGAAACGAAAAACCCGCTGGAAAAATATCTACAGGATAACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAg                                                                                                                              <  2:1145508/150‑1 (MQ=255)
                          attcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAgcgc                                                                                                                           >  1:945750/1‑151 (MQ=255)
                           ttctgggggAACGAAAAACCCGGTGAAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGCCTCCCCCTAACAAAATACCCTTCTTCCCCGGAGGAGGGCCCTTACCTTTATTCCCCCGCCGTTAAATCATCCACCGAg                                                                                                                              <  2:131206/141‑1 (MQ=255)
                           ttcTCGGGGACCGCCAAACCCGCGGGAAAAATTTCTACGGGATTAAAAAAAAAAAAAGCCTACCTTCGCCTCCCCCTAAAAAAAGACCCCTCTCCCCACGAAGGGGGCCGCAAACCTTAGTACCCCGCCGTTAAATCATCCACCGAgc                                                                                                                             <  2:84319/148‑1 (MQ=255)
                                  ggAACGGAAAAACCGGGGGAAAAATATCTACAGGATAACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTGAAAAAATGCCCCTCTCCCCACGAGGGGGGCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGa                                                                                                                   <  2:1190346/151‑1 (MQ=255)
                                    aaCGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATg                                                                                                                 >  1:161755/1‑151 (MQ=255)
                                      cGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcg                                                                                                               >  2:178054/1‑151 (MQ=255)
                                         cAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCg                                                                                                            >  1:658682/1‑151 (MQ=255)
                                         aaaacCCGCTGGAAAAATATCAAAGGAATTAAAAAAAAAAAAGGGCTCCCTTCGGCTGCCCCTAAAAAAACACCCCTCTCCCCACGAGGGGGCCCGCTACCCTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCg                                                                                                            <  1:17550/150‑1 (MQ=255)
                                          aaacCCGCTGGAAAAATATCAACAGAATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCAAAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                                <  2:1120946/146‑1 (MQ=255)
                                          aaacCCGCTGAAAAATTACCTACAGGATTACAAAAAAAAAAGGCCTCCCTCCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAGGGGGCCCGCTACCTTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                           <  2:215214/151‑1 (MQ=255)
                                              ccGTTGGATAAATATCAAAAAGATAAAAAAAAAAAAAGGTCTCCCTTCGGCTTCCCCTAACAAAATGCCCCCCTCCCCACGAAGGGGCCCCCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                        <  1:1115102/150‑1 (MQ=255)
                                              ccGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                        >  1:633245/1‑150 (MQ=255)
                                               cGCTGAAAAAATATCAACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTCCCCCAAAAAAAACCCCCCTCTCCCCACGAGGGGGCCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                                                      <  1:1059227/151‑1 (MQ=255)
                                                gCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                                     >  2:414836/1‑151 (MQ=255)
                                                     aagaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                                >  2:651972/1‑151 (MQ=255)
                                                     aagaaTATCTACAGGATAACAAAAAAAAAAAGCCTCCCTCCGGCTCGCCCGGAAAAAATACCTCTCTCCCCACAAAGGGGGCCGCTAACTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                 <  2:825896/150‑1 (MQ=255)
                                                     aaaaaTATATAAAGAATTAAAAAAAAAAAAAGCCTCCCTTGGCCTTCCCCTAAAAAAATCCCCCCCTCCCCACAAAGAGGCCCGCAACCCTAAGCCCCCCGCCGTTAAATCATCCACCTAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                                <  1:1056050/148‑1 (MQ=255)
                                                        aaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                             >  1:178550/1‑151 (MQ=255)
                                                         atatCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                            >  1:799114/1‑151 (MQ=255)
                                                             cTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                        >  2:204396/1‑151 (MQ=255)
                                                              tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                        >  2:380799/1‑150 (MQ=255)
                                                               aCAGAATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGACAAAATACCCCTCTCCCCACAAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGcc                                                                                         <  2:35148/148‑1 (MQ=255)
                                                                cAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                     >  2:656169/1‑151 (MQ=255)
                                                                  gtaTTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCAACCCTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                     <  2:1103844/147‑1 (MQ=255)
                                                                     ttACAAAAAAAAAAAGCCTCCCTCCGGCTTCCCCTGAAAAAATCCCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCa                                                                                   <  2:766318/148‑1 (MQ=255)
                                                                       aaaaaaaaaaaaaGGCTCCCTTCGCCTTCCCCTAAAAAAATAGCCCTCTCCCCACGAGGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                              <  1:1136358/148‑1 (MQ=255)
                                                                         gAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                            >  1:715455/1‑151 (MQ=255)
                                                                            aTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                         >  2:461853/1‑151 (MQ=255)
                                                                              aaaaaaGGCTCCTTCCGCCTCCCCCTAAAAAAATACCCCTCCCCCCACAAAGAGGCCCGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                       <  2:949858/151‑1 (MQ=255)
                                                                              aaaaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                       >  2:514281/1‑151 (MQ=255)
                                                                                   aGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                                  >  2:102501/1‑151 (MQ=255)
                                                                                     gCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGcc                                                                >  2:157533/1‑151 (MQ=255)
                                                                                       tACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                               >  1:1068716/1‑150 (MQ=255)
                                                                                             cGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                        >  2:931453/1‑151 (MQ=255)
                                                                                                 ttGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                    >  1:13120/1‑151 (MQ=255)
                                                                                                    ccctaaaaaaaacccccTCTCCCCAGGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctct                                                      <  2:193909/132‑1 (MQ=255)
                                                                                                    cccTAAAAAAATACCCCTCCCCCCACAAGGAGGGCCGCACTCTTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                    <  1:699107/148‑1 (MQ=255)
                                                                                                       tGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                              >  2:14179/1‑151 (MQ=255)
                                                                                                           aaaagccccccCTCCCCAGAAAGAGGCCCGTTACCCTTGCTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACt                                          <  1:1078415/144‑1 (MQ=255)
                                                                                                            aaaTAGCCCTCTCCCCACAAGGGGGGCCGCTATCTTTAATCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACtt                                         <  2:300330/151‑1 (MQ=255)
                                                                                                               tAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                      >  1:681966/1‑151 (MQ=255)
                                                                                                                acccccCTCCCCACGAAGGGGCCCGCAACCCTTATCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                      <  2:776910/148‑1 (MQ=255)
                                                                                                                        ccccgggAAGCGGCCCGCTACCTTTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCttt                               <  1:147387/143‑1 (MQ=255)
                                                                                                                         cccAGGAGGAGGCCCGCCACCCTTATCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAcc                                  <  1:532799/145‑1 (MQ=255)
                                                                                                                         cccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGc                            >  1:1111365/1‑151 (MQ=255)
                                                                                                                         cccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGc                            >  2:789796/1‑151 (MQ=255)
                                                                                                                         cccACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtt                                <  2:476810/147‑1 (MQ=255)
                                                                                                                          ccaggagggggcccccTACCCTTACTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                           <  1:597784/136‑1 (MQ=255)
                                                                                                                             cgagggggcccttcaaCCTTGGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                          <  2:899542/135‑1 (MQ=255)
                                                                                                                             cGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                        >  2:1204300/1‑151 (MQ=255)
                                                                                                                                   gggCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagcag                   <  2:878522/150‑1 (MQ=255)
                                                                                                                                    ggCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCt                 >  2:126159/1‑151 (MQ=255)
                                                                                                                                    ggCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCt                 >  2:556230/1‑151 (MQ=255)
                                                                                                                                    ggCCCCTAACCTTTTGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTc                         <  1:6683/144‑1 (MQ=255)
                                                                                                                                        gCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg               >  1:1126189/1‑149 (MQ=255)
                                                                                                                                        cgatatctttAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa              <  2:127080/143‑1 (MQ=255)
                                                                                                                                          tAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCg           >  2:131586/1‑151 (MQ=255)
                                                                                                                                           acccTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc            <  1:170056/147‑1 (MQ=255)
                                                                                                                                           acccTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc            <  2:945600/147‑1 (MQ=255)
                                                                                                                                             ccTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagcag                   <  1:695693/140‑1 (MQ=255)
                                                                                                                                              cTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATc       >  1:96598/1‑151 (MQ=255)
                                                                                                                                                ttAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATcc      >  1:431112/1‑151 (MQ=255)
                                                                                                                                                 tAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCg     >  2:710363/1‑151 (MQ=255)
                                                                                                                                                    tACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGc  >  1:11679/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGC  >  NC_000913/3584938‑3585235

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