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

TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGC  >  NC_000913/3584938‑3585235
                                                                                                                                                     |                                                                                                                                                     
tattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTAcc                                                                                                                                                     >  2:967880/1‑151 (MQ=255)
   tGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCccg                                                                                                                                                  >  1:89477/1‑151 (MQ=255)
    gATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTAcccc                                                                                                                                                   <  2:859045/149‑1 (MQ=255)
     aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCccgcc                                                                                                                                                <  1:854980/151‑1 (MQ=255)
              gtggtgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCTGTTAAatca                                                                                                                                       <  2:271097/151‑1 (MQ=255)
               tggtgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAatcat                                                                                                                                      >  1:769572/1‑151 (MQ=255)
                      attcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCg                                                                                                                                >  2:1078469/1‑150 (MQ=255)
                        tcatccTCTGGAAACGAAAAACCCGCTGGAAAAATATCTACAGGATAACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTAACC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAg                                                                                                                              <  2:1145508/150‑1 (MQ=255)
                          attcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAgcgc                                                                                                                           >  1:945750/1‑151 (MQ=255)
                           ttctgggggAACGAAAAACCCGGTGAAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGCCTCCCCCTAACAAAATACCCTTCTTCCCCGGAGGAGGGCCCTTACCT‑TTATTCCCCCGCCGTTAAATCATCCACCGAg                                                                                                                              <  2:131206/141‑1 (MQ=255)
                           ttcTCGGGGACCGCCAAACCCGCGGGAAAAATTTCTACGGGATTAAAAAAAAAAAAAGCCTACCTTCGCCTCCCCCTAAAAAAAGACCCCTCTCCCCACGAAGGGGGCCGCAAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAgc                                                                                                                             <  2:84319/148‑1 (MQ=255)
                                  ggAACGGAAAAACCGGGGGAAAAATATCTACAGGATAACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTGAAAAAATGCCCCTCTCCCCACGAGGGGGGCCGCTACCT‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGa                                                                                                                   <  2:1190346/151‑1 (MQ=255)
                                    aaCGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATg                                                                                                                 >  1:161755/1‑151 (MQ=255)
                                      cGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcg                                                                                                               >  2:178054/1‑151 (MQ=255)
                                         cAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCg                                                                                                            >  1:658682/1‑151 (MQ=255)
                                         aaaacCCGCTGGAAAAATATCAAAGGAATTAAAAAAAAAAAAGGGCTCCCTTCGGCTGCCCCTAAAAAAACACCCCTCTCCCCACGAGGGGGCCCGCTACCC‑TAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCg                                                                                                            <  1:17550/150‑1 (MQ=255)
                                          aaacCCGCTGGAAAAATATCAACAGAATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCAAAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCC‑TTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                                <  2:1120946/146‑1 (MQ=255)
                                          aaacCCGCTGAAAAATTACCTACAGGATTACAAAAAAAAAAGGCCTCCCTCCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAGGGGGCCCGCTACCT‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                           <  2:215214/151‑1 (MQ=255)
                                              ccGTTGGATAAATATCAAAAAGATAAAAAAAAAAAAAGGTCTCCCTTCGGCTTCCCCTAACAAAATGCCCCCCTCCCCACGAAGGGGCCCCCTACCT‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                        <  1:1115102/150‑1 (MQ=255)
                                              ccGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                        >  1:633245/1‑150 (MQ=255)
                                               cGCTGAAAAAATATCAACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTCCCCCAAAAAAAACCCCCCTCTCCCCACGAGGGGGCCCGCTACCC‑TTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                                                      <  1:1059227/151‑1 (MQ=255)
                                                gCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                                     >  2:414836/1‑151 (MQ=255)
                                                     aagaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                                >  2:651972/1‑151 (MQ=255)
                                                     aagaaTATCTACAGGATAACAAAAAAAAAAAGCCTCCCTCCGGCTCGCCCGGAAAAAATACCTCTCTCCCCACAAAGGGGGCCGCTAACT‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                 <  2:825896/150‑1 (MQ=255)
                                                     aaaaaTATATAAAGAATTAAAAAAAAAAAAAGCCTCCCTTGGCCTTCCCCTAAAAAAATCCCCCCCTCCCCACAAAGAGGCCCGCAACCC‑TAAGCCCCCCGCCGTTAAATCATCCACCTAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                                <  1:1056050/148‑1 (MQ=255)
                                                        aaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                             >  1:178550/1‑151 (MQ=255)
                                                         atatCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                            >  1:799114/1‑151 (MQ=255)
                                                             cTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                        >  2:204396/1‑151 (MQ=255)
                                                              tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                        >  2:380799/1‑150 (MQ=255)
                                                               aCAGAATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGACAAAATACCCCTCTCCCCACAAAGAGGCCCGCTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGcc                                                                                         <  2:35148/148‑1 (MQ=255)
                                                                cAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                     >  2:656169/1‑151 (MQ=255)
                                                                  gtaTTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCAACCC‑TAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                     <  2:1103844/147‑1 (MQ=255)
                                                                     ttACAAAAAAAAAAAGCCTCCCTCCGGCTTCCCCTGAAAAAATCCCCCTCTCCCCACGAAGAGGCCCGCTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCa                                                                                   <  2:766318/148‑1 (MQ=255)
                                                                       aaaaaaaaaaaaaGGCTCCCTTCGCCTTCCCCTAAAAAAATAGCCCTCTCCCCACGAGGAGGCCCGCTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                              <  1:1136358/148‑1 (MQ=255)
                                                                         gAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                            >  1:715455/1‑151 (MQ=255)
                                                                            aTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                         >  2:461853/1‑151 (MQ=255)
                                                                              aaaaaaGGCTCCTTCCGCCTCCCCCTAAAAAAATACCCCTCCCCCCACAAAGAGGCCCGCAACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                       <  2:949858/151‑1 (MQ=255)
                                                                              aaaaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                       >  2:514281/1‑151 (MQ=255)
                                                                                   aGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                                  >  2:102501/1‑151 (MQ=255)
                                                                                     gCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGcc                                                                >  2:157533/1‑151 (MQ=255)
                                                                                       tACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                               >  1:1068716/1‑150 (MQ=255)
                                                                                             cGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                        >  2:931453/1‑151 (MQ=255)
                                                                                                 ttGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                    >  1:13120/1‑151 (MQ=255)
                                                                                                    ccctaaaaaaaacccccTCTCCCCAGGAAGGGGCCCGCTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctct                                                      <  2:193909/132‑1 (MQ=255)
                                                                                                    cccTAAAAAAATACCCCTCCCCCCACAAGGAGGGCCGCACTCT‑TTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                    <  1:699107/148‑1 (MQ=255)
                                                                                                       tGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcgc                                              >  2:14179/1‑151 (MQ=255)
                                                                                                           aaaagccccccCTCCCCAGAAAGAGGCCCGTTACCC‑TTGCTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACt                                          <  1:1078415/144‑1 (MQ=255)
                                                                                                            aaaTAGCCCTCTCCCCACAAGGGGGGCCGCTATCT‑TTAATCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACtt                                         <  2:300330/151‑1 (MQ=255)
                                                                                                               tAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                      >  1:681966/1‑151 (MQ=255)
                                                                                                                acccccCTCCCCACGAAGGGGCCCGCAACCC‑TTATCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                      <  2:776910/148‑1 (MQ=255)
                                                                                                                        ccccgggAAGCGGCCCGCTACCT‑TTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCttt                               <  1:147387/143‑1 (MQ=255)
                                                                                                                         cccAGGAGGAGGCCCGCCACCC‑TTATCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAcc                                  <  1:532799/145‑1 (MQ=255)
                                                                                                                         cccACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGc                            >  1:1111365/1‑151 (MQ=255)
                                                                                                                         cccACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGc                            >  2:789796/1‑151 (MQ=255)
                                                                                                                         cccACGAAGAGGCCCGCTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtt                                <  2:476810/147‑1 (MQ=255)
                                                                                                                          ccaggagggggcccccTACCC‑TTACTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                           <  1:597784/136‑1 (MQ=255)
                                                                                                                             cgagggggcccttcaaCC‑TTGGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                          <  2:899542/135‑1 (MQ=255)
                                                                                                                             cGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                        >  2:1204300/1‑151 (MQ=255)
                                                                                                                                   gggCCGCTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagcag                   <  2:878522/150‑1 (MQ=255)
                                                                                                                                    ggCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCt                 >  2:126159/1‑151 (MQ=255)
                                                                                                                                    ggCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCt                 >  2:556230/1‑151 (MQ=255)
                                                                                                                                    ggCCCCTAACCTTTTGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTc                         <  1:6683/144‑1 (MQ=255)
                                                                                                                                        gCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg               >  1:1126189/1‑149 (MQ=255)
                                                                                                                                        cgatatctttAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGa              <  2:127080/143‑1 (MQ=255)
                                                                                                                                          tAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCg           >  2:131586/1‑151 (MQ=255)
                                                                                                                                           accc‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc            <  1:170056/147‑1 (MQ=255)
                                                                                                                                           accc‑TAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc            <  2:945600/147‑1 (MQ=255)
                                                                                                                                             cc‑TTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagcag                   <  1:695693/140‑1 (MQ=255)
                                                                                                                                              c‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATc       >  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)
                                                                                                                                                     |                                                                                                                                                     
TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGC  >  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: ‑