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

ATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  NC_000913/3584936‑3585236
                                                                                                                                                           |                                                                                                                                                      
attattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTa                                                                                                                                                        <  2:656301/151‑1 (MQ=255)
attattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTa                                                                                                                                                        <  2:1063009/151‑1 (MQ=255)
 ttattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAACGCTACCTTCGGCTTGCCCTGACAAAATAGCCATCTTCCCACGAAGAGGGCCGCTAACCTTAGTa                                                                                                                                                        >  2:1204460/1‑150 (MQ=255)
      gATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccg                                                                                                                                                   <  1:721048/150‑1 (MQ=255)
       aTAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgcc                                                                                                                                                 <  1:1019488/151‑1 (MQ=255)
          aaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGtt                                                                                                                                              >  2:437896/1‑151 (MQ=255)
                    tgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATcc                                                                                                                                     >  1:67011/1‑150 (MQ=255)
                     gAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCAc                                                                                                                                   >  2:930382/1‑151 (MQ=255)
                      aaaTTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCAcc                                                                                                                                  >  1:441174/1‑151 (MQ=255)
                          tcattcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAgc                                                                                                                              >  2:620514/1‑151 (MQ=255)
                             ttcTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGcgcg                                                                                                                           >  1:250943/1‑151 (MQ=255)
                              tctcTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCAGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTCAATCATACACCGAGCGCgg                                                                                                                          >  1:299229/1‑151 (MQ=255)
                                tctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGt                                                                                                                        >  1:1130045/1‑151 (MQ=255)
                                            aaaccccgtgaaaaaaTACCTAAAGGATTACAAAAAAAAAAGGGCCCCCTCCGCTTCCCCTTAAAAAAACACCCCTCTCCCCACGAGGGGGGCCCCAACCCTATGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGcc                                                                                                              <  1:1049832/135‑1 (MQ=255)
                                            aaacCCGGTGAAAAAATATCGAAAGGATTAAAAAAAAAAAAAGGCTCCCTTCGGCTTCCCCAAAAAAAATCCCCCCCTCCCCACGAAGGGGCCCGTTCACCTTGTTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc                                                                                                                 <  1:1048555/146‑1 (MQ=255)
                                            aaacCCGCGGAAAGAATATCTACAGAATTACAAAAAAAAAAAGCCTCCCTCCGGCTTCCCCTAAAAAAATAGCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGcc                                                                                                              <  2:301437/149‑1 (MQ=255)
                                            aaaaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            >  2:676281/1‑151 (MQ=255)
                                            aaaaCCGCTGGAAAAATATCTACAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTGAAAAAATAGCCCTCTCCCCACGAAGGGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGt                                                                                                            <  2:1228348/151‑1 (MQ=255)
                                                ccGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa                                                                                                         >  1:81120/1‑150 (MQ=255)
                                                  gCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACaa                                                                                                        >  2:1078367/1‑149 (MQ=255)
                                                  gCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt                                                                                                      >  2:1011978/1‑151 (MQ=255)
                                                    tGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                                     >  2:712877/1‑150 (MQ=255)
                                                       aagaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                                 >  2:247073/1‑151 (MQ=255)
                                                       aaaaaTATCTACGGAATTCCAAAAAAAAAATGTCTACCTTCTGTTCGCCGTGAAAAAATCCCCCTCTCCCCAGGAAGAGTGCCCTTACGCTTACCCCTCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  1:237079/147‑1 (MQ=255)
                                                       aaaaaTATCTACAGGATAACAAAAAAAAAAAGGCTCCCTGCGGTTTCCCCTAAAAAAATCGCCCTCTCCCCACGAAGGGGGCCGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  1:813206/147‑1 (MQ=255)
                                                       aaaaaTATCAACAGGTTAACAAAAAAAAAAAGCCTCCCTTCGGTTCCTCCAAAAAAAAACCCCCTCTCCCCACGAAGAGGGCCCCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAcc                                                                                                  <  1:1068228/147‑1 (MQ=255)
                                                       aaaaaTATCAAAAGGATTAAAAAAAAAAAAGGAATCCCTTGGGCTTCCCAAAAAAAAACCCCCCTCTCCCCAGGAGGGGGGCCGCTACCATTGGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc                                                                                                     <  2:868520/144‑1 (MQ=255)
                                                          aagaccaccAGAATAACAAAAAAAAAAAGGCTTCCGTCGGTTGGCCCTAAAAAAATAGCCCGCTCCCCACGGGAGGGCCCGCTACCCTTGCTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTc                                                                                               <  1:849050/142‑1 (MQ=255)
                                                            tatCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCgg                                                                                            >  2:426965/1‑151 (MQ=255)
                                                             atCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGc                                                                                           >  1:1113488/1‑151 (MQ=255)
                                                              tCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGcc                                                                                          >  1:126013/1‑151 (MQ=255)
                                                                tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                         >  2:714369/1‑150 (MQ=255)
                                                                 aCAGGATTAAAAAAAAAAAAAGGCTCCCTTCGGCTCCCCCTAACAAAATCCCCCTCTCCCCACGAGGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                                         <  1:700176/149‑1 (MQ=255)
                                                                 aCAGGATGACAAAAAAAAAAAGGTGCCCTTCGCCTTCCCCGAAAAAAATACCCCTCTCCCCGGGAGGGGGCCCAATACCCGTATCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCaa                                                                                       <  1:231223/151‑1 (MQ=255)
                                                                   aGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                     >  1:713287/1‑151 (MQ=255)
                                                                      aTTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCAt                                                                                   >  1:386936/1‑150 (MQ=255)
                                                                       ttaaaaaaaaaaaaaGGCTACCTTCGGTTCCCCCGGAAAAAAAACCCCCCTCCCCACGAAGAGGCCCGCAACCTTTATTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCag                                                                                   <  2:1110965/144‑2 (MQ=255)
                                                                       ttACAAAAAAAAAAAGGCTCCCTCCGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGGCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATa                                                                                  <  2:362402/150‑1 (MQ=255)
                                                                           aaaaaaaaaaaGCCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGAGGCCCGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                             <  2:156888/150‑1 (MQ=255)
                                                                           aaaaTAAAAAGGGCTCCCTTGGTCTTCCCGTAAAAAATTCGCCCTCTCCCCACAAAGGGCCCCGCTACCCTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGc                                                                              <  2:1160121/149‑1 (MQ=255)
                                                                            aaaaaaaaaaGGCCTCCCTTCGCCTTCCCCTAAAAAAATACCCCTCTCCCCACGAGGGGGCCCGCTAACCTTCGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCt                                                                             <  1:915600/151‑1 (MQ=255)
                                                                            aaaaaaaaaaGGCTCCCTCCGGCTTCCCCTAAAAAAATACCCCTCCCCCCACAAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGc                                                                              <  2:799646/149‑1 (MQ=255)
                                                                            aaaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  2:551950/1‑151 (MQ=255)
                                                                             aaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGTGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                                            >  1:224183/1‑150 (MQ=255)
                                                                               tAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg                                                                          >  2:135723/1‑150 (MQ=255)
                                                                               aaaaaaaGGCTCCCTTCGGCTTGCCCTAAAAAAATCGCCCTCTCCCCACGAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                                         <  2:409167/150‑1 (MQ=255)
                                                                                      ggCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGc                                                                  >  1:1230063/1‑151 (MQ=255)
                                                                                       gCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGACCGGCCCAACCATAATGCTTTGCGCACTGcc                                                                 >  1:1207999/1‑151 (MQ=255)
                                                                                         cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  1:1157047/1‑151 (MQ=255)
                                                                                         cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGccc                                                                >  2:299771/1‑151 (MQ=255)
                                                                                            ccTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                                             >  2:990760/1‑151 (MQ=255)
                                                                                               tCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGCTGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGc                                                           >  1:1072625/1‑150 (MQ=255)
                                                                                                 ggtttccccTAAAAAAATACCCCTTTCCCCACAAAGAGGCCCCTTACCCTTATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                                                        <  2:138258/145‑1 (MQ=255)
                                                                                                 ggCTTGCCCAAAAAAAAGAGCCCTCTCCCCAAAAAGGGGCCCGCTAACCGTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                         <  2:606422/150‑1 (MQ=255)
                                                                                                   cTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                      >  2:907581/1‑151 (MQ=255)
                                                                                                      gCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                                                    >  1:304707/1‑150 (MQ=255)
                                                                                                       cccTGAAAAAATACCCCTCTCCCCACGAAGGGGCCCCAAACCTTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCtt                                                      <  1:747572/147‑1 (MQ=255)
                                                                                                        cccaaaaaaaaccccccTCTCCCCACGAAGAGGCCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                         <  2:1123277/130‑1 (MQ=255)
                                                                                                            gACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGcc                                              >  2:776702/1‑151 (MQ=255)
                                                                                                                 aaagccccTCCCCCCACGAAGAGGCCCTCTACCCTTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                                         <  2:771654/146‑1 (MQ=255)
                                                                                                                     ccccTCTCCCCACGAAGAGGGCCGCTACCCTTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                       <  1:246610/148‑1 (MQ=255)
                                                                                                                     ccccTCTCCCCAAGAGGAGGGCGGCAACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                       <  1:67193/148‑1 (MQ=255)
                                                                                                                        ctctTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCt                                  >  2:26970/1‑151 (MQ=255)
                                                                                                                          ctccccACGAAGAGGCCCGCTACCCTTACTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCt                                  <  1:920607/146‑1 (MQ=255)
                                                                                                                           ttCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCtttt                               >  1:108004/1‑151 (MQ=255)
                                                                                                                               cACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                           >  2:169919/1‑151 (MQ=255)
                                                                                                                                caaaggagGGCCCCCTAACCGTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCttt                                <  2:479013/142‑1 (MQ=255)
                                                                                                                                   aaGGGGCCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt                            <  2:1213412/146‑1 (MQ=255)
                                                                                                                                    agagGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCcagc                      >  2:1029165/1‑151 (MQ=255)
                                                                                                                                     gagGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCca                        <  2:970975/148‑1 (MQ=255)
                                                                                                                                        ggCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCt                  >  1:577623/1‑151 (MQ=255)
                                                                                                                                         gCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCtt                 >  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)
                                                                                                                                                           |                                                                                                                                                      
ATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCT  >  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: