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

TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCG  >  NC_000913/3584938‑3585232
                                                                                                                                                      |                                                                                                                                                  
tattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcc                                                                                                                                                  >  2:537402/1‑151 (MQ=255)
tattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcc                                                                                                                                                  >  2:1241528/1‑151 (MQ=255)
 attGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAccc                                                                                                                                                 <  1:250659/151‑1 (MQ=255)
 attGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcc                                                                                                                                                  <  2:1074825/150‑1 (MQ=255)
   tGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccg                                                                                                                                               >  2:1136046/1‑151 (MQ=255)
    gATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgc                                                                                                                                              <  1:1295676/151‑1 (MQ=255)
    gATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccg                                                                                                                                               <  2:740672/150‑1 (MQ=255)
      tAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgccg                                                                                                                                            >  2:285653/1‑151 (MQ=255)
       aaaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgccg                                                                                                                                            <  1:1099109/150‑1 (MQ=255)
            cAGTGGCGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAat                                                                                                                                      >  2:630735/1‑151 (MQ=255)
              gtggtgAAATTCTTTCTCTGGGAACGGCAAAACCGCGGGAAAAATATCTACAGAATTACAAAAAAAAAAAGCCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTaaa                                                                                                                                       <  2:389997/148‑1 (MQ=255)
                     aaTTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCg                                                                                                                             >  1:168970/1‑151 (MQ=255)
                        tcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAgc                                                                                                                          >  1:1260038/1‑151 (MQ=255)
                          attcTGGGGGACCGGCAAAACCGCAGAAAAAGTATCTCCAGGATTACAAAAAAAAAACGCCCCCCTTCTGCTGCCCCTCACAAAACGGCCCTCTTCCACCAAGGAGGGCCGTTACCTTCACTCCCCCGCCGTTAAATCATCCACCGAgc                                                                                                                          <  2:1262164/149‑1 (MQ=255)
                            tctcTGGGAACGAAAAACCCGGGGGAAAAATATCTACAGGATTACGAAAAAAAAAAGCCTCCCTTCGCCTTCCCCTAAAAAAATACCCCTCTCCCCACGAGGAGGCCCACTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAgcg                                                                                                                         <  1:460727/148‑1 (MQ=255)
                             ctctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCggg                                                                                                                     >  1:921836/1‑151 (MQ=255)
                              tctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGt                                                                                                                    >  1:1239471/1‑151 (MQ=255)
                               ctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTc                                                                                                                   >  1:252562/1‑151 (MQ=255)
                                tGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCg                                                                                                                  >  2:14841/1‑151 (MQ=255)
                                    aaCGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATg                                                                                                              >  1:1249091/1‑151 (MQ=255)
                                     aCGGAAAACCCGGTGAAAAAATATCTACAGAATTACAAAAAAAAAAAGTCTACCTTCGGCTCCCCCTAAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGAt                                                                                                               <  1:1093031/149‑1 (MQ=255)
                                          aaacCCGCGGAAAAAATATCTAAAGGATAACAAAAAAAAAAAGCCTCCCTTCGGCTCCCCCTAAAAAAAAACCCCTCTCCCCACGAAGGGGCCCGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATg                                                                                                              <  2:746186/145‑1 (MQ=255)
                                           aacccGCGGGAAAAATACCAAAAGGATAAAAAAAAAAAAAAGGCTCCTTTCGGTTTCCCCTGAAAAAATACCCCTCTCCCCAAGAAGGGGGCCCCTAACCTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcgc                                                                                                           <  2:558974/144‑1 (MQ=255)
                                             aCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTAc                                                                                                      >  1:407915/1‑150 (MQ=255)
                                               cGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                                                   >  1:603273/1‑151 (MQ=255)
                                                gCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                                  >  2:492513/1‑151 (MQ=255)
                                                   ggAAAAATATCTAAAGGATAACAAAAAAAAAAAGGCTCCCTGCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa                                                                                                <  2:1246124/150‑1 (MQ=255)
                                                       gaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTc                                                                                           >  1:401481/1‑151 (MQ=255)
                                                         atatCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAAGGGGCCCGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                         <  2:86190/151‑1 (MQ=255)
                                                              tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                     >  1:374660/1‑150 (MQ=255)
                                                               aCAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTAAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCaa                                                                                   <  1:541214/151‑1 (MQ=255)
                                                               aCAGAATTAAAAAAAAAAAAGGGCTCCTTGCGGCTTCCCCTAAAAAACTCCCCCTCTCCCCGGAAGGGGGCCCCCTACCTTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCaa                                                                                   <  1:1160092/151‑1 (MQ=255)
                                                                 agaatatcaaaaaaaaaaaGTTTACTTTGGGCTTCCCTAAAAAAAATCCCCCCCTCCCCACAAAGAGACCCGAAAACTTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa                                                                                    <  2:561258/139‑1 (MQ=255)
                                                                  ggATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                 >  1:841249/1‑150 (MQ=255)
                                                                  gaataacaaaaaaaaaaaGCCTCCCTTGGCCTTGCCCTAAAAAAATCCCCCTCTCCCCACAAAGGGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                  <  1:1288697/141‑1 (MQ=255)
                                                                  gaataacaaaaaaaaaaaGCCTCCCTGGGTCTTGCCCTGAAAAAATACCCCTCTCCCCACAAGGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa                                                                                    <  1:1157916/139‑1 (MQ=255)
                                                                      tACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                                            >  1:1215694/1‑151 (MQ=255)
                                                                       aaaaaaaaaaaaaGCCTCCCTGCGCCTCCCCCTAAAAAAATCCCCCTCTCCCCACGAAGGGGGCCCCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATaa                                                                             <  1:818974/146‑1 (MQ=255)
                                                                          aaaaaaaaaaGCCTCCCTCCGGTTTCCCCCAAAAAAATCCCCCTTTCCCCAGGAGGGGGGCCCCTATCCTTAGACCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg                                                                           <  1:776379/148‑1 (MQ=255)
                                                                          aaaaaaaaAGGGTTGCCTGCGGCTGGCCTGGAAAAAAGCCCCCTCCCCCCAGGAGGGGGGCGGCAAGCCGTAGGCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                        <  1:387087/151‑1 (MQ=255)
                                                                            aaaaaaaaGGCTCCCTTGGGCTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTAACCTAAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                      <  2:1155420/149‑1 (MQ=255)
                                                                            aTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                      >  2:586407/1‑151 (MQ=255)
                                                                             aaaaaaaGGCTCCCTTCGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTACCTTTATTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                     <  2:799448/150‑1 (MQ=255)
                                                                              aaaaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                    >  1:9774/1‑151 (MQ=255)
                                                                                aaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTgtgt                                                                   >  1:914139/1‑150 (MQ=255)
                                                                                 aaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                 >  1:301768/1‑151 (MQ=255)
                                                                                      ctcccTTCGGCTTGCCCTGAAAAAATACCCCTCCCCCCACGAAGAGGGCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTg                                                               <  1:324741/145‑1 (MQ=255)
                                                                                      cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                            >  1:280476/1‑151 (MQ=255)
                                                                                       acccTCCGGCTTGCCCTAACAAAACCCCCCTCTCCCCCGGAGGAGGCCCGCTACCTTAGGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                              <  1:704371/146‑1 (MQ=255)
                                                                                        aCCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCAt                                                          >  1:1071503/1‑151 (MQ=255)
                                                                                         ccTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                                         >  1:718532/1‑151 (MQ=255)
                                                                                         ccTTCGGCTTGCCCTAAAAAAATACCCCCCTCCCCACGAAGGGGGCCCCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                                         <  2:1073968/151‑1 (MQ=255)
                                                                                          cTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                        >  2:492694/1‑151 (MQ=255)
                                                                                              ggCTTCCCCTAAAAAAATCCCCCTCTCCCCAGGAGGGGGCCCGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCg                                                        <  1:817142/147‑1 (MQ=255)
                                                                                               gCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctct                                                   >  1:1115427/1‑151 (MQ=255)
                                                                                                cTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                  >  1:786224/1‑151 (MQ=255)
                                                                                                   gCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                  <  1:1134357/148‑1 (MQ=255)
                                                                                                   ccccTGAAAAAATGCCCCTCTCCCCACGAAGAGGGCCCCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCa                                               <  2:655739/150‑1 (MQ=255)
                                                                                                   ccccTAAAAAAACAGCCCTCTCCCCCCGAGGGGGCCCGCTACCTTGAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                <  1:585334/149‑1 (MQ=255)
                                                                                                    cccTGAAAAAAGAGCCCTCTCCCCACAAAGAGGCCCCCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCa                                               <  1:350777/150‑1 (MQ=255)
                                                                                                      cggaaaaaaTGGCCCTCTCCCCACGAAGTGGGCCCCGACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcg                                            <  1:461798/146‑1 (MQ=255)
                                                                                                          cAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCAc                                        >  2:347130/1‑151 (MQ=255)
                                                                                                             aaTACCCCTCTCCCCCCGAAGGGGGCCGCTACCCTGAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                     <  1:362709/151‑1 (MQ=255)
                                                                                                              atcccccTCTCCCCACAAAGAGGTCCGCAAACCTGAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTc                                    <  2:1252884/147‑1 (MQ=255)
                                                                                                              aTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTc                                    >  2:174756/1‑151 (MQ=255)
                                                                                                                   ccTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAcc                               >  2:806949/1‑151 (MQ=255)
                                                                                                                         cccACGAAGGGGGCCGCTACCTTTACTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtttt                           <  1:504382/149‑1 (MQ=255)
                                                                                                                          ccACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTg                          <  1:552313/149‑1 (MQ=255)
                                                                                                                           cACGAAGGGGCCCGCTACCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                       <  2:1084793/151‑1 (MQ=255)
                                                                                                                               aaGGGGCCCGCAACCCTTACTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                       <  1:35663/147‑1 (MQ=255)
                                                                                                                               aaGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                     >  1:379111/1‑149 (MQ=255)
                                                                                                                               aaGAGGCCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                        <  1:1007644/146‑1 (MQ=255)
                                                                                                                                agagGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc                  >  1:870744/1‑151 (MQ=255)
                                                                                                                                      ccGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCtt             <  2:642173/150‑1 (MQ=255)
                                                                                                                                      ccGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg            <  1:891430/151‑1 (MQ=255)
                                                                                                                                       cGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg            >  1:630165/1‑150 (MQ=255)
                                                                                                                                           acccTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGc         <  1:630685/147‑1 (MQ=255)
                                                                                                                                           acccTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt       <  2:417245/149‑1 (MQ=255)
                                                                                                                                           aaCCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt       >  1:1054825/1‑151 (MQ=255)
                                                                                                                                              cTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATc    >  2:942946/1‑151 (MQ=255)
                                                                                                                                                  aGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCg  <  2:1209252/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCG  >  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: