Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,0860AC45.6% 105.8 / 21.8 69Y580D (TAC→GAC) ggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (30/7);  new base C (0/31);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 7.44e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584937‑3585236
                                                                                                                                                        |                                                                                                                                                       
ttattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAc                                                                                                                                                       >  1:1077604/1‑151 (MQ=255)
ttattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAc                                                                                                                                                       >  1:543468/1‑151 (MQ=255)
  attGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAccc                                                                                                                                                     >  2:285834/1‑151 (MQ=255)
    tGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcccc                                                                                                                                                    <  2:695890/150‑1 (MQ=255)
    tGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcccc                                                                                                                                                    <  2:55147/150‑1 (MQ=255)
            aCAATGGCGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGGCGAATATCTACAGCATTACGAAATAGAAAAGGCTACCGGCGGCTTGGCCTGCCAAAATAGCCCACTGCCCACGAAGAGGGCCGCAAACCTTAGTACCCCGCCGTTa                                                                                                                                             >  2:211596/1‑149 (MQ=255)
             cAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAat                                                                                                                                          >  2:357997/1‑151 (MQ=255)
                      aaTTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCg                                                                                                                                 >  2:1141189/1‑151 (MQ=255)
                        ttctttctCTGGAACCGGAAAACCCGCTGGAAAAATATCTACAGTATTACAAAAAAAAAAAGCCTCCCTTCGGCTTCCCCTGAAAAAATCCCCCTCTCCCCACGAAGGGGCCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGa                                                                                                                                <  1:259431/150‑1 (MQ=255)
                        ttctttctCGGGGAACGAAAACCCCGCTGAAAAAATATCTACAGAATTAAAAAAAAAAAAGGTCTACTTTCTACTCCCCCTAAAAAAACACCCCTTCCCCCACAAGGGGGCCCACTAACTTTACTCCCCCGCCGTTAAATCATCCAc                                                                                                                                   <  2:553097/147‑1 (MQ=255)
                             ttctcgggAAAGGAAAAACCCGTTGAAAAATTATCTACAGGATTACAAAAAAAAAAAGCTTCCCTCCGGCTTGCCCTGAAAAAATCCCCCTCTCCCCACAAAGAGGCCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAgc                                                                                                                              <  1:318850/143‑1 (MQ=255)
                             tctcTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCgg                                                                                                                          >  1:285835/1‑151 (MQ=255)
                                  gggAACGGCAAACCCGCGGGAAAAATATCTACAGGATTAAAAAAAAAAAAAGGCTACCTTCGGCTTCCCCTGAAAAAATAGCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                                     <  1:351767/151‑1 (MQ=255)
                                  gggAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                                     >  1:742275/1‑151 (MQ=255)
                                  gggAACGGAAAACCCGTGGGAAAAATATCTAAAGGATTAAAAAAAAAAAAAGCCTCCCTCCGGTTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTACCCTAAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCg                                                                                                                      <  1:1091419/150‑1 (MQ=255)
                                  gggAACGGAAAACCCGCTGAAAGAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTCCGGCTTGCCCTAAAAAAATAGCCCTCTCCCCACGAAGAGGCCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCgg                                                                                                                          <  1:1104609/146‑1 (MQ=255)
                                     aaCGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATg                                                                                                                  >  1:317192/1‑151 (MQ=255)
                                      acgcaaaaaCCGGTGAAAAAATATCTACAGGATTACAAAATAAAAAAGGCTACCTTCGGCTTCCCCTAAAAAAATAGCCCTCCCCCCACGAAGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                                 <  1:944465/146‑1 (MQ=255)
                                      aCGGCAAACCCGCTGGAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTAAAAAAATCCCCCTCTCCCCACGAGGAGGGCCGCTATCTTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGTAt                                                                                                                   <  2:571779/149‑1 (MQ=255)
                                           aaacCCGTTGGAAAAATATCTACAGGATTAAAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTAAAAAAATCCCCCTCTCCCCCCGAAGGGGGCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcgc                                                                                                               <  2:344072/148‑1 (MQ=255)
                                           aaaaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCAGATGCGCCGt                                                                                                            >  2:186050/1‑151 (MQ=255)
                                           aaaaCCGAGGAAAAAATATCTAAAGAATAACAAAAAAAAAAGGGCCCCCTTCGGCTTCCCCTGAAAAAATCCCCCTCTCCCCACAAAGGGGGCCCCTAACTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                                 <  2:728220/146‑1 (MQ=255)
                                                ccGCTGAAAGAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                         <  1:44206/150‑1 (MQ=255)
                                                     ggAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                                                   >  2:529269/1‑151 (MQ=255)
                                                       aagaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                                 >  1:417675/1‑151 (MQ=255)
                                                        aaaaTATCTACAGAATTACAAAAAAAAAAGGGCTCCCTCCGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGCCCCCAACCTTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGt                                                                                                <  1:54654/149‑1 (MQ=255)
                                                          aaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                              >  1:330738/1‑151 (MQ=255)
                                                           atatCTACAGGATTACAAAAAAAAAAAGGCTACCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAGGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                             <  1:1117614/151‑1 (MQ=255)
                                                              tCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGcc                                                                                          >  1:1166313/1‑151 (MQ=255)
                                                                tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                         >  1:457988/1‑150 (MQ=255)
                                                                 acaaattcaCGAAATAAAAAAGGCTCCATCCGGGTTCCCCCGAAAAAATCCCCCCCGCCCAACAAGAGGGGCGGTTAACGTTAGACCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGcc                                                                                          <  1:123879/140‑1 (MQ=255)
                                                                 aCAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGACAAAATACCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa                                                                                        <  1:387180/150‑1 (MQ=255)
                                                                  cAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                      >  1:742371/1‑151 (MQ=255)
                                                                  cAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                      >  1:938573/1‑151 (MQ=255)
                                                                    ggATTACGAAAAAAAAAAGGCTCCCTCCGGCTTCCCCGAACAAAATACCCCTCTCCCCACAAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     <  1:65915/150‑1 (MQ=255)
                                                                    gattaaaaaaaaaaaaaaGGCTCCCTTCGCCTTCCCCTGAAAAAATAGCCCCCTGCCAACAAAGTGGGCCGCTAGCTTGAGCGCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     <  1:426152/142‑1 (MQ=255)
                                                                            aaaaaaaaaaGGCTCCCTGGGGCTGCCCCTAAAAAAATAGCCCCCTCCCCACGAGGAGGGCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                                                <  1:793252/147‑1 (MQ=255)
                                                                            aaaaaaaaaaGGCCTCCTTTCGCTTCCCCCGGAAAAAAGCGCCCTCTCCCCACGAGGGGGGCCTCTACCCTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                             <  2:745095/151‑1 (MQ=255)
                                                                                   aaaGGCTACCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                     <  1:191418/151‑1 (MQ=255)
                                                                                    aaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                     >  1:1033714/1‑150 (MQ=255)
                                                                                         tcccTCCGGCTTCCCCTAACAAAATACCCCTCTCCCCACGAAGGGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                                  <  1:494799/146‑1 (MQ=255)
                                                                                         tACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  2:47919/1‑150 (MQ=255)
                                                                                             ttCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  2:723570/1‑151 (MQ=255)
                                                                                             ttCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  2:170730/1‑151 (MQ=255)
                                                                                             ttCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  1:1020540/1‑151 (MQ=255)
                                                                                                ggCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                        >  1:202901/1‑151 (MQ=255)
                                                                                                gccTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTACCCTAAGCCCCCCGCCGTTAAATCATCCACCAAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGcc                                                          <  2:292023/147‑1 (MQ=255)
                                                                                                    tGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCTGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                     <  1:814007/150‑1 (MQ=255)
                                                                                                       cctaaaCAAAATCCCCCTCTCCCCAAGAAGAGGCCCGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCa                                                   <  2:67175/145‑1 (MQ=255)
                                                                                                       cccTGAAAAAATCGCCCTCTCCCCAGGAAGAGGCCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                      <  1:470644/147‑1 (MQ=255)
                                                                                                       cccTGAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                     <  1:1128379/148‑1 (MQ=255)
                                                                                                                     ccttctcCCCACGAAGAGGGCGCTTATTCTGCGCACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         <  2:883430/143‑1 (MQ=255)
                                                                                                                     cccTCTCCCCCAGAAGGGGCCCGCTAACTTTATCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTg                                      <  2:1102370/149‑1 (MQ=255)
                                                                                                                      ccTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAcc                                   >  2:448856/1‑151 (MQ=255)
                                                                                                                      ccTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAcc                                   <  1:826851/151‑1 (MQ=255)
                                                                                                                            cccAGGAGGGGGCCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGc                             <  2:902637/151‑1 (MQ=255)
                                                                                                                             ccACGAGGGGGGCCGCTAACCTTGGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            <  1:305381/151‑1 (MQ=255)
                                                                                                                             ccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            >  2:727476/1‑151 (MQ=255)
                                                                                                                             ccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            >  2:83753/1‑151 (MQ=255)
                                                                                                                               aCGAAGGGGCCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTc                          <  1:8548/151‑1 (MQ=255)
                                                                                                                               aCGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTc                          >  2:923409/1‑151 (MQ=255)
                                                                                                                                 gAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                         >  2:464672/1‑150 (MQ=255)
                                                                                                                                  aaGAGGGCCGCTAACCTTGTACCCCCGCCGTTAAATCATCCACCGAGCGCGGGCCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            <  1:920696/146‑1 (MQ=255)
                                                                                                                                  aaGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                         <  1:565134/149‑1 (MQ=255)
                                                                                                                                     gggggCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcag                       <  2:242788/147‑1 (MQ=255)
                                                                                                                                       ggcccacTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCagcagc                   <  1:1134727/145‑1 (MQ=255)
                                                                                                                                          ccGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg                >  1:641035/1‑151 (MQ=255)
                                                                                                                                                   ttAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATcc       >  1:852937/1‑151 (MQ=255)
                                                                                                                                                    tAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTa          <  1:733005/147‑1 (MQ=255)
                                                                                                                                                        aCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTCCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  1:1086779/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                       
TTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584937‑3585236

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: