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

TTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584937‑3585236
                                                                                                                                                        |                                                                                                                                                       
ttattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑c                                                                                                                                                       >  1:1077604/1‑151 (MQ=255)
ttattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑c                                                                                                                                                       >  1:543468/1‑151 (MQ=255)
  attGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑ccc                                                                                                                                                     >  2:285834/1‑151 (MQ=255)
    tGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑cccc                                                                                                                                                    <  2:695890/150‑1 (MQ=255)
    tGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑cccc                                                                                                                                                    <  2:55147/150‑1 (MQ=255)
            aCAATGGCGAAATTCATTCTCTGGGAACGGCAAAA‑CCGCTGGGCGAATATCTACAGCATTACGAAATAGAAAAGGCTACCGGCGGCTT‑GGCCTGCCAAAATAGCCCACTGCCCACGAAGA‑GGGCCGCAAACCTTAGTA‑CCCCGCCGTTa                                                                                                                                             >  2:211596/1‑149 (MQ=255)
             cAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAat                                                                                                                                          >  2:357997/1‑151 (MQ=255)
                      aaTTCATTCTCTGGGAACGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCg                                                                                                                                 >  2:1141189/1‑151 (MQ=255)
                        ttctttctCTGGAACCGGAAAAC‑CCGCTGGAAAAATATCTACAGTATTACAAAAAAAAAAAGCCTCCCTTCGGCTT‑CCCCTGAAAAAATCCCCCTCTCCCCACGAAGG‑GGCCCGCTAACCTTAGTC‑CCCCGCCGTTAAATCATCCACCGa                                                                                                                                <  1:259431/150‑1 (MQ=255)
                        ttctttctCGGGGAACGAAAACC‑CCGCTGAAAAAATATCTACAGAATTAAAAAAAAAAAAGGTCTACTTTCTACTC‑CCCCTAAAAAAACACCCCTTCCCCCACAAGGG‑GGCCCACTAACTTTACTC‑CCCCGCCGTTAAATCATCCAc                                                                                                                                   <  2:553097/147‑1 (MQ=255)
                             ttctcgggAAAGGAAAAACCCGTTGAAAAATTATCTACAGGATTACAAAAAAAAAAAGCTTCCCTCCGGCTT‑GCCCTGAAAAAATCCCCCTCTCCCCACAAAGA‑GGCCCGCTACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAgc                                                                                                                              <  1:318850/143‑1 (MQ=255)
                             tctcTGGGAACGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCgg                                                                                                                          >  1:285835/1‑151 (MQ=255)
                                  gggAACGGCAAAC‑CCGCGGGAAAAATATCTACAGGATTAAAAAAAAAAAAAGGCTACCTTCGGCTT‑CCCCTGAAAAAATAGCCCTCTCCCCACGAAGA‑GGCCCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                                     <  1:351767/151‑1 (MQ=255)
                                  gggAACGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg                                                                                                                     >  1:742275/1‑151 (MQ=255)
                                  gggAACGGAAAAC‑CCGTGGGAAAAATATCTAAAGGATTAAAAAAAAAAAAAGCCTCCCTCCGGTTT‑CCCCTAAAAAAATACCCCTCTCCCCACGAAGA‑GGGCCGCTACCCTAAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCg                                                                                                                      <  1:1091419/150‑1 (MQ=255)
                                  gggAACGGAAAAC‑CCGCTGAAAGAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTCCGGCTT‑GCCCTAAAAAAATAGCCCTCTCCCCACGAAGA‑GGCCCGCTACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCgg                                                                                                                          <  1:1104609/146‑1 (MQ=255)
                                     aaCGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATg                                                                                                                  >  1:317192/1‑151 (MQ=255)
                                      acgcaaaaa‑CCGGTGAAAAAATATCTACAGGATTACAAAATAAAAAAGGCTACCTTCGGCTT‑CCCCTAAAAAAATAGCCCTCCCCCCACGAAGA‑GGGCCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                                 <  1:944465/146‑1 (MQ=255)
                                      aCGGCAAAC‑CCGCTGGAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTT‑GCCCTAAAAAAATCCCCCTCTCCCCACGAGGA‑GGGCCGCTATCTTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGTAt                                                                                                                   <  2:571779/149‑1 (MQ=255)
                                           aaac‑CCGTTGGAAAAATATCTACAGGATTAAAAAAAAAAAAAGGCTCCCTTCGGCTT‑GCCCTAAAAAAATCCCCCTCTCCCCCCGAAGG‑GGGCCGCTACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcgc                                                                                                               <  2:344072/148‑1 (MQ=255)
                                           aaaa‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCAGATGCGCCGt                                                                                                            >  2:186050/1‑151 (MQ=255)
                                           aaaa‑CCGAGGAAAAAATATCTAAAGAATAACAAAAAAAAAAGGGCCCCCTTCGGCTT‑CCCCTGAAAAAATCCCCCTCTCCCCACAAAGG‑GGGCCCCTAACTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                                 <  2:728220/146‑1 (MQ=255)
                                                ccGCTGAAAGAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTT‑CCCCTAAAAAAATACCCCTCTCCCCACGAAGA‑GGGCCGCTACCTTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                         <  1:44206/150‑1 (MQ=255)
                                                     ggAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                                                   >  2:529269/1‑151 (MQ=255)
                                                       aagaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                                 >  1:417675/1‑151 (MQ=255)
                                                        aaaaTATCTACAGAATTACAAAAAAAAAAGGGCTCCCTCCGGCTT‑GCCCTGAAAAAATACCCCTCTCCCCACGAAGA‑GGCCCCCAACCTTTATTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGt                                                                                                <  1:54654/149‑1 (MQ=255)
                                                          aaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc                                                                                              >  1:330738/1‑151 (MQ=255)
                                                           atatCTACAGGATTACAAAAAAAAAAAGGCTACCTTCGGCTT‑CCCCTGAAAAAATACCCCTCTCCCCACGAGGA‑GGGCCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                             <  1:1117614/151‑1 (MQ=255)
                                                              tCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGcc                                                                                          >  1:1166313/1‑151 (MQ=255)
                                                                tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                         >  1:457988/1‑150 (MQ=255)
                                                                 acaaattcaCGAAATAAAAAAGGCTCCATCCGGGTT‑CCCCCGAAAAAAT‑CCCCCCCGCCCAACAAGAGGGGCGGTTAACGTTAGAC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGcc                                                                                          <  1:123879/140‑1 (MQ=255)
                                                                 aCAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTT‑CCCCTGACAAAATACCCCTCTCCCCACGAAGG‑GGCCCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa                                                                                        <  1:387180/150‑1 (MQ=255)
                                                                  cAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                      >  1:742371/1‑151 (MQ=255)
                                                                  cAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc                                                                                      >  1:938573/1‑151 (MQ=255)
                                                                    ggATTACGAAAAAAAAAAGGCTCCCTCCGGCTT‑CCCCGAACAAAATACCCCTCTCCCCACAAAGA‑GGCCCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     <  1:65915/150‑1 (MQ=255)
                                                                    gattaaaaaaaaaaaaaaGGCTCCCTTCGCCTT‑CCCCTGAAAAAATAGCCCCCTGCCAACAAAGT‑GGGCCGCTAGCTTGAGCG‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     <  1:426152/142‑1 (MQ=255)
                                                                            aaaaaaaaaaGGCTCCCTGGGGCTG‑CCCCTAAAAAAATAGCCCCCTCCCCACGAGGA‑GGGCCGCTAACCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                                                <  1:793252/147‑1 (MQ=255)
                                                                            aaaaaaaaaaGGCCTCCTTTCGCTTCCCCCGGAAAAAAGCGCCCTCTCCCCACGAGGG‑GGGCCTCTACCCTTATTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                             <  2:745095/151‑1 (MQ=255)
                                                                                   aaaGGCTACCTTCGGCTT‑CCCCTGAAAAAATACCCCTCTCCCCACGAAGA‑GGGCCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                     <  1:191418/151‑1 (MQ=255)
                                                                                    aaGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                                                     >  1:1033714/1‑150 (MQ=255)
                                                                                         tcccTCCGGCTT‑CCCCTAACAAAATACCCCTCTCCCCACGAAGG‑GGGCCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                                  <  1:494799/146‑1 (MQ=255)
                                                                                         tACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  2:47919/1‑150 (MQ=255)
                                                                                             ttCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  2:723570/1‑151 (MQ=255)
                                                                                             ttCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  2:170730/1‑151 (MQ=255)
                                                                                             ttCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  1:1020540/1‑151 (MQ=255)
                                                                                                ggCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                        >  1:202901/1‑151 (MQ=255)
                                                                                                gccTT‑CCCCTAAAAAAATACCCCTCTCCCCACGAAGA‑GGGCCGCTACCCTAAGCC‑CCCCGCCGTTAAATCATCCACCAAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGcc                                                          <  2:292023/147‑1 (MQ=255)
                                                                                                    t‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCTGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                     <  1:814007/150‑1 (MQ=255)
                                                                                                       cctaaaCAAAATCCCCCTCTCCCCAAGAAGA‑GGCCCGCAACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCa                                                   <  2:67175/145‑1 (MQ=255)
                                                                                                       cccTGAAAAAATCGCCCTCTCCCCAGGAAGA‑GGCCCGCTACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                      <  1:470644/147‑1 (MQ=255)
                                                                                                       cccTGAAAAAATACCCCTCTCCCCACGAAGA‑GGCCCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                                                     <  1:1128379/148‑1 (MQ=255)
                                                                                                                     ccttctcCCCACGAAGA‑GGGCGCTTATTCTGCGCA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         <  2:883430/143‑1 (MQ=255)
                                                                                                                     cccTCTCCCCCAGAAGG‑GGCCCGCTAACTTTATCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTg                                      <  2:1102370/149‑1 (MQ=255)
                                                                                                                      ccTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAcc                                   >  2:448856/1‑151 (MQ=255)
                                                                                                                      ccTCTCCCCACGAAGG‑GGCCCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGAcc                                   <  1:826851/151‑1 (MQ=255)
                                                                                                                            cccAGGAGGG‑GGCCCGCTACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGc                             <  2:902637/151‑1 (MQ=255)
                                                                                                                             ccACGAGGG‑GGGCCGCTAACCTTGGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            <  1:305381/151‑1 (MQ=255)
                                                                                                                             ccACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            >  2:727476/1‑151 (MQ=255)
                                                                                                                             ccACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            >  2:83753/1‑151 (MQ=255)
                                                                                                                               aCGAAGG‑GGCCCGCTAACCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTc                          <  1:8548/151‑1 (MQ=255)
                                                                                                                               aCGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTc                          >  2:923409/1‑151 (MQ=255)
                                                                                                                                 gAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                         >  2:464672/1‑150 (MQ=255)
                                                                                                                                  aaGA‑GGGCCGCTAACCTT‑GTACCCCCGCCGTTAAATCATCCACCGAGCGCGGGCCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            <  1:920696/146‑1 (MQ=255)
                                                                                                                                  aaGA‑GGCCCGCTACCCTTAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTcc                         <  1:565134/149‑1 (MQ=255)
                                                                                                                                     gggggCCGCTACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcag                       <  2:242788/147‑1 (MQ=255)
                                                                                                                                       ggcccacTACCCTTAGCC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCagcagc                   <  1:1134727/145‑1 (MQ=255)
                                                                                                                                          ccGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTg                >  1:641035/1‑151 (MQ=255)
                                                                                                                                                   ttAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATcc       >  1:852937/1‑151 (MQ=255)
                                                                                                                                                    tAGTC‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTa          <  1:733005/147‑1 (MQ=255)
                                                                                                                                                        a‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTCCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  1:1086779/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                       
TTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAA‑CCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTT‑GCCCTGACAAAATAGCCCTCTTCCCACGAAGA‑GGGCCGCTAACCTTAGTA‑CCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  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: ‑