Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,0860AC44.4% 126.7 / 17.2 81Y580D (TAC→GAC) ggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (37/8);  new base C (0/36);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.42e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCG  >  NC_000913/3584938‑3585232
                                                                                                                                                      |                                                                                                                                                  
tattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTAcc                                                                                                                                                  >  2:537402/1‑151 (MQ=255)
tattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTAcc                                                                                                                                                  >  2:1241528/1‑151 (MQ=255)
 attGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTAccc                                                                                                                                                 <  1:250659/151‑1 (MQ=255)
 attGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTAcc                                                                                                                                                  <  2:1074825/150‑1 (MQ=255)
   tGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCccg                                                                                                                                               >  2:1136046/1‑151 (MQ=255)
    gATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCccgc                                                                                                                                              <  1:1295676/151‑1 (MQ=255)
    gATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCccg                                                                                                                                               <  2:740672/150‑1 (MQ=255)
      tAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCccgccg                                                                                                                                            >  2:285653/1‑151 (MQ=255)
       aaaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCccgccg                                                                                                                                            <  1:1099109/150‑1 (MQ=255)
            cAGTGGCGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAat                                                                                                                                      >  2:630735/1‑151 (MQ=255)
              gtggtgAAATTCTTTCTCTGGGAACGGCAAAACCGCGGGAAAAATATCTACAGAATTACAAAAAAAAAAAGCCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCCTT‑A‑GTCCCCCGCCGTTaaa                                                                                                                                       <  2:389997/148‑1 (MQ=255)
                     aaTTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCg                                                                                                                             >  1:168970/1‑151 (MQ=255)
                        tcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAgc                                                                                                                          >  1:1260038/1‑151 (MQ=255)
                          attcTGGGGGACCGGCAAAACCGCAGAAAAAGTATCTCCAGGATTACAAAAAAAAAACGCCCCCCTTCTGCTGCCCCTCACAAAACGGCCCTCTTCCACCAAGGAGGGCCGTTACCTTC‑A‑CTCCCCCGCCGTTAAATCATCCACCGAgc                                                                                                                          <  2:1262164/149‑1 (MQ=255)
                            tctcTGGGAACGAAAAACCCGGGGGAAAAATATCTACAGGATTACGAAAAAAAAAAGCCTCCCTTCGCCTTCCCCTAAAAAAATACCCCTCTCCCCACGAGGAGGCCCACTACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAgcg                                                                                                                         <  1:460727/148‑1 (MQ=255)
                             ctctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCggg                                                                                                                     >  1:921836/1‑151 (MQ=255)
                              tctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGt                                                                                                                    >  1:1239471/1‑151 (MQ=255)
                               ctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTc                                                                                                                   >  1:252562/1‑151 (MQ=255)
                                tGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCg                                                                                                                  >  2:14841/1‑151 (MQ=255)
                                    aaCGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATg                                                                                                              >  1:1249091/1‑151 (MQ=255)
                                     aCGGAAAACCCGGTGAAAAAATATCTACAGAATTACAAAAAAAAAAAGTCTACCTTCGGCTCCCCCTAAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTAACCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGAt                                                                                                               <  1:1093031/149‑1 (MQ=255)
                                          aaacCCGCGGAAAAAATATCTAAAGGATAACAAAAAAAAAAAGCCTCCCTTCGGCTCCCCCTAAAAAAAAACCCCTCTCCCCACGAAGGGGCCCGCAACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATg                                                                                                              <  2:746186/145‑1 (MQ=255)
                                           aacccGCGGGAAAAATACCAAAAGGATAAAAAAAAAAAAAAGGCTCCTTTCGGTTTCCCCTGAAAAAATACCCCTCTCCCCAAGAAGGGGGCCCCTAACCTA‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcgc                                                                                                           <  2:558974/144‑1 (MQ=255)
                                             aCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTAc                                                                                                      >  1:407915/1‑150 (MQ=255)
                                               cGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                                                   >  1:603273/1‑151 (MQ=255)
                                                gCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                                  >  2:492513/1‑151 (MQ=255)
                                                   ggAAAAATATCTAAAGGATAACAAAAAAAAAAAGGCTCCCTGCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTAACCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa                                                                                                <  2:1246124/150‑1 (MQ=255)
                                                       gaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTc                                                                                           >  1:401481/1‑151 (MQ=255)
                                                         atatCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAAGGGGCCCGCAACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                         <  2:86190/151‑1 (MQ=255)
                                                              tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                     >  1:374660/1‑150 (MQ=255)
                                                               aCAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTAAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCaa                                                                                   <  1:541214/151‑1 (MQ=255)
                                                               aCAGAATTAAAAAAAAAAAAGGGCTCCTTGCGGCTTCCCCTAAAAAACTCCCCCTC‑TCCCCGGAAGGGGGCCCCCTACCTTAA‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCaa                                                                                   <  1:1160092/151‑1 (MQ=255)
                                                                 agaatatcaaaaaaaaaaaGTTTACTTTGGGCTTCCCTAAAAAAAATCCCCCCCTCCCCACAAAGAGACCCGAAAACTTT‑A‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa                                                                                    <  2:561258/139‑1 (MQ=255)
                                                                  ggATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                 >  1:841249/1‑150 (MQ=255)
                                                                  gaataacaaaaaaaaaaaGCCTCCCTTGGCCTTGCCCTAAAAAAATCCCCCTCTCCCCACAAAGGGGGCCGCTACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                  <  1:1288697/141‑1 (MQ=255)
                                                                  gaataacaaaaaaaaaaaGCCTCCCTGGGTCTTGCCCTGAAAAAATACCCCTCTCCCCACAAGGGGGCCCGCTACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa                                                                                    <  1:1157916/139‑1 (MQ=255)
                                                                      tACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                                            >  1:1215694/1‑151 (MQ=255)
                                                                       aaaaaaaaaaaaaGCCTCCCTGCGCCTCCCCCTAAAAAAATCCCCCTCTCCCCACGAAGGGGGCCCCTACCCTT‑A‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATaa                                                                             <  1:818974/146‑1 (MQ=255)
                                                                          aaaaaaaaaaGCCTCCCTCCGGTTTCCCCCAAAAAAATCCCCCTTTCCCCAGGAGGGGGGCCCCTATCCTT‑A‑GACCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                           <  1:776379/148‑1 (MQ=255)
                                                                          aaaaaaaaAGGGTTGCCTGCGGCTGGCCTGGAAAAAAGCCCCCTCCCCCCAGGAGGGGGGCGGCAAGCCGT‑A‑GGCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                        <  1:387087/151‑1 (MQ=255)
                                                                            aaaaaaaaGGCTCCCTTGGGCTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTAACCTA‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                      <  2:1155420/149‑1 (MQ=255)
                                                                            aTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                      >  2:586407/1‑151 (MQ=255)
                                                                             aaaaaaaGGCTCCCTTCGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTACCTTT‑A‑TTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                     <  2:799448/150‑1 (MQ=255)
                                                                              aaaaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                    >  1:9774/1‑151 (MQ=255)
                                                                                aaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTgtgt                                                                   >  1:914139/1‑150 (MQ=255)
                                                                                 aaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                 >  1:301768/1‑151 (MQ=255)
                                                                                      ctcccTTCGGCTTGCCCTGAAAAAATACCCCTCCCCCCACGAAGAGGGCCGCTACCCTT‑A‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                               <  1:324741/145‑1 (MQ=255)
                                                                                      cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                            >  1:280476/1‑151 (MQ=255)
                                                                                       acccTCCGGCTTGCCCTAACAAAACCCCCCTCTCCCCCGGAGGAGGCCCGCT‑ACCTT‑AGGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                              <  1:704371/146‑1 (MQ=255)
                                                                                        aCCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCAt                                                          >  1:1071503/1‑151 (MQ=255)
                                                                                         ccTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                                         >  1:718532/1‑151 (MQ=255)
                                                                                         ccTTCGGCTTGCCCTAAAAAAATACCCCCCTCCCCACGAAGGGGGCCCCTAACCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                                         <  2:1073968/151‑1 (MQ=255)
                                                                                          cTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                        >  2:492694/1‑151 (MQ=255)
                                                                                              ggCTTCCCCTAAAAAAATCCCCCTCTCCCCAGGAGGGGGCCCGCAACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                        <  1:817142/147‑1 (MQ=255)
                                                                                               gCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctct                                                   >  1:1115427/1‑151 (MQ=255)
                                                                                                cTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                  >  1:786224/1‑151 (MQ=255)
                                                                                                   gCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                  <  1:1134357/148‑1 (MQ=255)
                                                                                                   ccccTGAAAAAATGCCCCTCTCCCCACGAAGAGGGCCCCAACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCa                                               <  2:655739/150‑1 (MQ=255)
                                                                                                   ccccTAAAAAAACAGCCCTCTCCCCCCGAGGGGGCCCGCTACCTTG‑A‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                <  1:585334/149‑1 (MQ=255)
                                                                                                    cccTGAAAAAAGAGCCCTCTCCCCACAAAGAGGCCCCCTACCCTT‑A‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCa                                               <  1:350777/150‑1 (MQ=255)
                                                                                                      cggaaaaaaTGGCCCTCTCCCCACGAAGTGGGCCCCGACCTTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcg                                            <  1:461798/146‑1 (MQ=255)
                                                                                                          cAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCAc                                        >  2:347130/1‑151 (MQ=255)
                                                                                                             aaTACCCCTCTCCCCCCGAAGGGGGCCGCTACCCTG‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                     <  1:362709/151‑1 (MQ=255)
                                                                                                              atcccccTCTCCCCACAAAGAGGTCCGCAAACCTG‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTc                                    <  2:1252884/147‑1 (MQ=255)
                                                                                                              aTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTc                                    >  2:174756/1‑151 (MQ=255)
                                                                                                                   ccTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAcc                               >  2:806949/1‑151 (MQ=255)
                                                                                                                         cccACGAAGGGGGCCGCTACCTTT‑A‑CTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtttt                           <  1:504382/149‑1 (MQ=255)
                                                                                                                          ccACGAAGAGGCCCGCTACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTg                          <  1:552313/149‑1 (MQ=255)
                                                                                                                           cACGAAGGGGCCCGCTACCCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                       <  2:1084793/151‑1 (MQ=255)
                                                                                                                               aaGGGGCCCGCAACCCTT‑A‑CTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                       <  1:35663/147‑1 (MQ=255)
                                                                                                                               aaGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                     >  1:379111/1‑149 (MQ=255)
                                                                                                                               aaGAGGCCCGCTACCCTT‑A‑GCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                        <  1:1007644/146‑1 (MQ=255)
                                                                                                                                agagGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc                  >  1:870744/1‑151 (MQ=255)
                                                                                                                                      ccGCTACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCtt             <  2:642173/150‑1 (MQ=255)
                                                                                                                                      ccGCTACCCTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg            <  1:891430/151‑1 (MQ=255)
                                                                                                                                       cGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg            >  1:630165/1‑150 (MQ=255)
                                                                                                                                           acccTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc         <  1:630685/147‑1 (MQ=255)
                                                                                                                                           acccTT‑A‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt       <  2:417245/149‑1 (MQ=255)
                                                                                                                                           aaCCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt       >  1:1054825/1‑151 (MQ=255)
                                                                                                                                              cTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATc    >  2:942946/1‑151 (MQ=255)
                                                                                                                                                  a‑GTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCg  <  2:1209252/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑A‑GTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCG  >  NC_000913/3584938‑3585232

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