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

ATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584936‑3585236
                                                                                                                                                           |                                                                                                                                                      
attattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTa                                                                                                                                                        <  2:656301/151‑1 (MQ=255)
attattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTa                                                                                                                                                        <  2:1063009/151‑1 (MQ=255)
 ttattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAACG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCATCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTa                                                                                                                                                        >  2:1204460/1‑150 (MQ=255)
      gATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCccg                                                                                                                                                   <  1:721048/150‑1 (MQ=255)
       aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCccgcc                                                                                                                                                 <  1:1019488/151‑1 (MQ=255)
          aaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGtt                                                                                                                                              >  2:437896/1‑151 (MQ=255)
                    tgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATcc                                                                                                                                     >  1:67011/1‑150 (MQ=255)
                     gAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCAc                                                                                                                                   >  2:930382/1‑151 (MQ=255)
                      aaaTTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCAcc                                                                                                                                  >  1:441174/1‑151 (MQ=255)
                          tcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAgc                                                                                                                              >  2:620514/1‑151 (MQ=255)
                             ttcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGcgcg                                                                                                                           >  1:250943/1‑151 (MQ=255)
                              tctcTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑AGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTCAATCATACACCGAGCGCgg                                                                                                                          >  1:299229/1‑151 (MQ=255)
                                tctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGt                                                                                                                        >  1:1130045/1‑151 (MQ=255)
                                            aaaccccgtgaaaaaaTACCTAAAGGATTAC‑AAAAAAAAAAGG‑GCCCCCTCCGCTTCCCCTTAAAAAAACACCCCTCTCCCCACGAGGGGGG‑CCC‑CAACCC‑TATGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGcc                                                                                                              <  1:1049832/135‑1 (MQ=255)
                                            aaacCCGGTGAAAAAATATCGAAAGGATTAAAAAAAAAAAAAGG‑CTCCCTTCGGCTTCCCC‑AAAAAAAATCCCCCCCTCCCCACGAAGGGGC‑CCG‑TTCACC‑TTGTTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                                 <  1:1048555/146‑1 (MQ=255)
                                            aaacCCGCGGAAAGAATATCTACAGAATTACAAAAAAAAAAAGC‑CTCCCTCCGGCTTCCCC‑TAAAAAAATAGCCCTCTCCCCACGAAGAGGC‑CCG‑CTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGcc                                                                                                              <  2:301437/149‑1 (MQ=255)
                                            aaaaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            >  2:676281/1‑151 (MQ=255)
                                            aaaaCCGCTGGAAAAATATCTACAGGATTACAAAAAAAAAAAGG‑CTCCCTTCGGCTTGCCC‑TGAAAAAATAGCCCTCTCCCCACGAAGGGGG‑CCG‑CTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            <  2:1228348/151‑1 (MQ=255)
                                                ccGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                         >  1:81120/1‑150 (MQ=255)
                                                  gCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACaa                                                                                                        >  2:1078367/1‑149 (MQ=255)
                                                  gCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                                      >  2:1011978/1‑151 (MQ=255)
                                                    tGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                                     >  2:712877/1‑150 (MQ=255)
                                                       aagaaTATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                                 >  2:247073/1‑151 (MQ=255)
                                                       aaaaaTATCTACGGAATTCCAAAAAAAAAATGT‑CTACCTTCTGTTCGCCG‑TGAAAAAATCCCCCTCTCCCCAGGAAGAGTG‑CCC‑TTACGC‑TTACCCCTCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  1:237079/147‑1 (MQ=255)
                                                       aaaaaTATCTACAGGATAACAAAAAAAAAAAGG‑CTCCCTGCGGTTTCCCC‑TAAAAAAATCGCCCTCTCCCCACGAAGGGGG‑CCG‑CAACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  1:813206/147‑1 (MQ=255)
                                                       aaaaaTATCAACAGGTTAACAAAAAAAAAAAGC‑CTCCCTTCGGTTCCTCC‑AAAAAAAAACCCCCTCTCCCCACGAAGAGGG‑CCC‑CTACCC‑TTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  1:1068228/147‑1 (MQ=255)
                                                       aaaaaTATCAAAAGGATTAAAAAAAAAAAAGGA‑ATCCCTTGGGCTTCCCA‑AAAAAAAACCCCCCTCTCCCCAGGAGGGGGG‑CCG‑CTACCA‑TTGGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                                     <  2:868520/144‑1 (MQ=255)
                                                          aagaccaccAGAATAACAAAAAAAAAAAGG‑CTTCCGTCGGTTGGCCC‑TAAAAAAATAGCCCGCTCCCCACG‑GGAGGGCCCG‑CTACCC‑TTGCTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTc                                                                                               <  1:849050/142‑1 (MQ=255)
                                                            tatCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCgg                                                                                            >  2:426965/1‑151 (MQ=255)
                                                             atCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGc                                                                                           >  1:1113488/1‑151 (MQ=255)
                                                              tCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGcc                                                                                          >  1:126013/1‑151 (MQ=255)
                                                                tACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                         >  2:714369/1‑150 (MQ=255)
                                                                 aCAGGATTAAAAAAAAAAAAAGG‑CTCCCTTCGGCTCCCCC‑TAACAAAATCCCCCTCTCCCCACGAGGGGGC‑CCG‑CTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                         <  1:700176/149‑1 (MQ=255)
                                                                 aCAGGATGACAAAAAAAAAAAGG‑TGCCCTTCGCCTTCCCC‑GAAAAAAATACCCCTCTCCCCGGGAGGGGGC‑CCA‑ATACCC‑GTATCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCaa                                                                                       <  1:231223/151‑1 (MQ=255)
                                                                   aGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     >  1:713287/1‑151 (MQ=255)
                                                                      aTTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCAt                                                                                   >  1:386936/1‑150 (MQ=255)
                                                                       ttaaaaaaaaaaaaaGG‑CTACCTTCGGTTCCCCC‑GGAAAAAAAACCCCCCTCCCCACGAAGAGGC‑CCG‑CAACCT‑TTATTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCag                                                                                   <  2:1110965/144‑2 (MQ=255)
                                                                       ttACAAAAAAAAAAAGG‑CTCCCTCCGGCTTGCCC‑TGAAAAAATACCCCTCTCCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATa                                                                                  <  2:362402/150‑1 (MQ=255)
                                                                           aaaaaaaaaaaGC‑CTCCCTTCGGCTTCCCC‑TGAAAAAATACCCCTCTCCCCACGAAGAGGC‑CCG‑CAACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                             <  2:156888/150‑1 (MQ=255)
                                                                           aaaaTAAAAAGGG‑CTCCCTTGGTCTTCCCG‑TAAAAAATTCGCCCTCTCCCCACAAAGGGCC‑CCG‑CTACCC‑TTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGc                                                                              <  2:1160121/149‑1 (MQ=255)
                                                                            aaaaaaaaaaGGCCTCCCTTCGCCTTCCCC‑TAAAAAAATACCCCTCTCCCCACGAGGGGGC‑CCG‑CTAACC‑TTCGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                             <  1:915600/151‑1 (MQ=255)
                                                                            aaaaaaaaaaGG‑CTCCCTCCGGCTTCCCC‑TAAAAAAATACCCCTCCCCCCACAAAGGGGC‑CCG‑CTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGc                                                                              <  2:799646/149‑1 (MQ=255)
                                                                            aaaTAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  2:551950/1‑151 (MQ=255)
                                                                             aaTAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGTGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  1:224183/1‑150 (MQ=255)
                                                                               tAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                          >  2:135723/1‑150 (MQ=255)
                                                                               aaaaaaaGG‑CTCCCTTCGGCTTGCCC‑TAAAAAAATCGCCCTCTCCCCACGAAGAGGC‑CCG‑CTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                         <  2:409167/150‑1 (MQ=255)
                                                                                      gg‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                                  >  1:1230063/1‑151 (MQ=255)
                                                                                       g‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGACCGGCCCAACCATAATGCTTTGCGCACTGcc                                                                 >  1:1207999/1‑151 (MQ=255)
                                                                                         cTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  1:1157047/1‑151 (MQ=255)
                                                                                         cTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  2:299771/1‑151 (MQ=255)
                                                                                            ccTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                                             >  2:990760/1‑151 (MQ=255)
                                                                                               tCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGCTGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  1:1072625/1‑150 (MQ=255)
                                                                                                 ggtttcccc‑TAAAAAAATACCCCTTTCCCCACAAAGAGGC‑CCC‑TTACCC‑TTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                        <  2:138258/145‑1 (MQ=255)
                                                                                                 ggCTTGCCC‑AAAAAAAAGAGCCCTCTCCCCAAAAAGGGGC‑CCG‑CTAACC‑GTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                         <  2:606422/150‑1 (MQ=255)
                                                                                                   cTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                      >  2:907581/1‑151 (MQ=255)
                                                                                                      gCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                    >  1:304707/1‑150 (MQ=255)
                                                                                                       ccc‑TGAAAAAATACCCCTCTCCCCACGAAG‑GGG‑CCC‑CAAACCTTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                      <  1:747572/147‑1 (MQ=255)
                                                                                                        cccaaaaaaaaccccccTCTCCCCACGAAGAGGC‑CCG‑CTAACC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                         <  2:1123277/130‑1 (MQ=255)
                                                                                                            gACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGcc                                              >  2:776702/1‑151 (MQ=255)
                                                                                                                 aaagccccTCCCCCCACGAAGAGGC‑CCT‑CTACCC‑TTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         <  2:771654/146‑1 (MQ=255)
                                                                                                                     ccccTCTCCCCACGAAGAGGG‑CCG‑CTACCC‑TTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                       <  1:246610/148‑1 (MQ=255)
                                                                                                                     ccccTCTCCCCAAGAGGAGGG‑CGG‑CAACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                       <  1:67193/148‑1 (MQ=255)
                                                                                                                        ctctTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCt                                  >  2:26970/1‑151 (MQ=255)
                                                                                                                          ctccccACGAAGAGGC‑CCG‑CTACCC‑TTACTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCt                                  <  1:920607/146‑1 (MQ=255)
                                                                                                                           ttCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtttt                               >  1:108004/1‑151 (MQ=255)
                                                                                                                               cACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                           >  2:169919/1‑151 (MQ=255)
                                                                                                                                caaaggagGG‑CCCCCTAACCGTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCttt                                <  2:479013/142‑1 (MQ=255)
                                                                                                                                   aaGGGGC‑CCG‑CTACCT‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            <  2:1213412/146‑1 (MQ=255)
                                                                                                                                    agagGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc                      >  2:1029165/1‑151 (MQ=255)
                                                                                                                                     gagGC‑CCG‑CTACCC‑TTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCca                        <  2:970975/148‑1 (MQ=255)
                                                                                                                                        gg‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCt                  >  1:577623/1‑151 (MQ=255)
                                                                                                                                         g‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCtt                 >  2:433928/1‑151 (MQ=255)
                                                                                                                                                         gTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCggg    >  2:887822/1‑151 (MQ=255)
                                                                                                                                                          aaCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGc   <  2:441195/150‑1 (MQ=255)
                                                                                                                                                           aCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCt  >  1:211351/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                      
ATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGG‑CTACCTTCGGCTTGCCC‑TGACAAAATAGCCCTCTTCCCACGAAGAGGG‑CCG‑CTAACC‑TTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584936‑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: ‑