Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	44.4%	96.2 / 18.5	63	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (30/5);  new base C (0/28);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 5.91e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537968‑4538253                                                                                                                                                       |                                                                                                                                                   aaaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTc                                                                                                                                         >  2:34763/1‑151 (MQ=255)  aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTca                                                                                                                                        <  2:116975/151‑1 (MQ=255)       gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAAcc                                                                                                                                   >  1:709476/1‑151 (MQ=255)       gCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAAcc                                                                                                                                   <  2:769660/151‑1 (MQ=255)               cGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCAcc                                                                                                                           >  2:53331/1‑151 (MQ=255)                      tGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATCATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGtt                                                                                                                    >  2:135949/1‑151 (MQ=255)                       ggTGATGAGCCGGGGGCAAGTTTATTCCATTTA‑AAATTATTCCC‑GGTAAAATAAAGTTAAAATCCGG‑CATCCTTTCGCAATCCTGTTT‑TGG‑CTC‑CGCAATAAAAAAC‑CAAGT‑TTTTAC‑CACT‑TTT‑ATCCACAACCGCCGCACCAGCCg                                                                                                                      <  2:944083/148‑1 (MQ=255)                       ggTAATGAGACGGAAGCAA‑CTTTTTCCATTTTAAAATTATACCG‑GGTAAAATAAATTTAAAATACGC‑CATCCTTTCCCAACCCTGGTT‑TGG‑CTT‑CGCCATTAAAAAG‑CAGGG‑TTTTAT‑CCCT‑TTT‑ATTCACAACCGCCGCACCAGCCg                                                                                                                      <  2:855257/148‑1 (MQ=255)                        gTGATGAGACGGGAGCAGGCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTcc                                                                                                                  >  2:266129/1‑151 (MQ=255)                          gatgaGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGt                                                                                                                >  1:857055/1‑151 (MQ=255)                                  gggAGCAAGTTTTTTCCATTTT‑AAATTTTTCCC‑GGTAAAATCAAGTTAAAATCCGG‑CATCCGTCCACAACCCGTTTT‑TGG‑TTT‑CGCCATAAAAAAC‑CCGGT‑TTTTAT‑CACC‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                          <  1:478607/149‑1 (MQ=255)                                  gggAGAAAGCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGCT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTCTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGc                                                                                                        >  2:890183/1‑151 (MQ=255)                                      gCAAGCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                    >  1:230565/1‑151 (MQ=255)                                        aaGCTTATTCCATTTA‑AAATTATTCCC‑GGTAAAATAAAGTTAAAATCCGG‑CATCCGTCCGCAATCCTGTTT‑GGG‑TTT‑CGCCTTAAATAAC‑CCAGT‑TTTTAC‑CCCT‑TTT‑ATTCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                    <  1:6145/149‑1 (MQ=255)                                         aagttttttcATTTT‑AAATTTATACCGGGTAAAATCAATTCAAAATCCGG‑CTCCCTTTCGAAATCCTGTTT‑TGG‑CTC‑CCCCATAAAAAAG‑CCGGT‑TTTTAC‑CCCT‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAt                                                                                                     <  1:1076685/139‑1 (MQ=255)                                          gCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                >  2:304071/1‑151 (MQ=255)                                            ttATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc                                                                                                >  2:1121175/1‑149 (MQ=255)                                                   aTTTT‑TAATTTTTGCC‑GGGTAAAAAAAGGTAAAAAACGGCCACCCTTGCGAACTCTTTGTT‑GGG‑CTT‑CGCCATTAA‑AAG‑CCGGGTTTTTAT‑CCCC‑TTT‑GTTCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGc                                                                                             <  2:357190/146‑1 (MQ=255)                                                      tt‑AAATTATTCCC‑GGTGAAATAAAGTTAAAATACGG‑CATCCGTCCGCAATCCTGGTT‑TGG‑CTT‑CCCCATTAATAAG‑CCAGT‑TTTTAT‑CCCC‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcc                                                                                     <  1:1112444/150‑1 (MQ=255)                                                         aaaTTATTCCC‑GGTAAAATCAAGTTAAAATCCGG‑CATCCGTTCGCAACCCTGTTT‑TGG‑TTT‑CGCCTTTAATAAC‑CCAGG‑TTTTAT‑CCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt                                                                                  <  1:845082/151‑1 (MQ=255)                                                          aaTTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTg                                                                                 >  2:562824/1‑151 (MQ=255)                                                             tattcCC‑GGTAAAAAAAAGTAAAAATGGGG‑CACCCGTTCGAAATCCTGTTT‑TGA‑TTC‑CGCCATTAAAAAG‑CCGGT‑TTTTAT‑CCCC‑TTT‑TTCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                               <  1:507595/150‑1 (MQ=255)                                                                taCC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                           >  1:549011/1‑151 (MQ=255)                                                                       tGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                     >  2:1052394/1‑151 (MQ=255)                                                                            tCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                 >  1:297109/1‑150 (MQ=255)                                                                              aatttaaaaTCCGG‑CCCCCTTCCGCATTCCT‑GTT‑TGG‑TTTCCGCCATAAAAAAC‑CCAGT‑TTTTAC‑CCCT‑TTT‑ATTCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                               <  2:932071/144‑1 (MQ=255)                                                                               aGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                               >  1:793714/1‑149 (MQ=255)                                                                                  tCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                           >  1:167529/1‑150 (MQ=255)                                                                                   cAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGaca                                                          >  2:285984/1‑150 (MQ=255)                                                                                     aaTACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                        >  2:744481/1‑150 (MQ=255)                                                                                           gccttccttccgaaaTCCGGGTT‑GGG‑TTC‑CGCAATTAAGAAG‑CCGGT‑TTTTAT‑CACG‑TTT‑ATCAACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                   <  2:411582/137‑1 (MQ=255)                                                                                           g‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                 >  1:254045/1‑151 (MQ=255)                                                                                              acccttcccaaacccTGGTT‑TGG‑C‑T‑CCCCAATAAAAAGCCCAGT‑TTTTAT‑CCCC‑TTT‑ACCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                <  2:792312/137‑1 (MQ=255)                                                                                                ccGTCCCCACTCCTGGTT‑TGC‑CTC‑CCCCTTTAATAAG‑CCAGT‑TTTTAT‑CCCC‑TTT‑ATTCACATCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                <  1:610553/148‑1 (MQ=255)                                                                                                 cGTTCCCAAACCTGGTTGTGG‑CTC‑CCCAATAAAAAAC‑CCGGT‑TTTTAC‑CCCC‑TTT‑TTCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACt                                                 <  2:290952/147‑1 (MQ=255)                                                                                                    tCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTga                                           >  2:260333/1‑149 (MQ=255)                                                                                                       cAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaagaa                                       >  2:591536/1‑150 (MQ=255)                                                                                                           ccTGTTT‑TGG‑CTA‑CCCCATAAAAAGG‑CCAGT‑TTTTTC‑CCTT‑TTT‑CTCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                  <  1:105094/151‑1 (MQ=255)                                                                                                              gttt‑tGG‑TTC‑CGCCATTAAAAAC‑CCGGT‑TTTTAT‑CCCC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                  <  1:940251/146‑1 (MQ=255)                                                                                                                tt‑tGGCCTT‑CGCAATTAAAAAC‑CCTGT‑TTTTAC‑CCCC‑TTT‑ATCCCCACCCGCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtc                                 <  1:885271/148‑1 (MQ=255)                                                                                                                tt‑tGG‑TTT‑CCCAATAAAAAAG‑CCGGT‑TTTTAC‑CCCT‑TTT‑ATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttc                              <  1:519806/150‑1 (MQ=255)                                                                                                                tt‑tGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                             >  1:648092/1‑151 (MQ=255)                                                                                                                tt‑tGG‑CTC‑CGCCATTAAAACG‑CCAGT‑TTTTAC‑CCCC‑TTT‑AACCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGaaa                                      <  2:948543/142‑1 (MQ=255)                                                                                                                 t‑tGG‑CTT‑CGCCATTAAAAAG‑CCAGT‑TTTTAC‑CCCT‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                             <  1:387432/150‑1 (MQ=255)                                                                                                                 t‑tGG‑ATT‑CGCCAATAAACAC‑GCATG‑TTTTAT‑CCCC‑TTTAACCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                               <  2:793510/149‑1 (MQ=255)                                                                                                                   tGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                           >  1:828042/1‑151 (MQ=255)                                                                                                                    gg‑CTT‑CCCCATAAATAAG‑CCGGT‑TTTTAT‑CCCT‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTg                            <  1:861444/149‑1 (MQ=255)                                                                                                                    gg‑CCT‑CCCCATTAAAAAG‑CCCGT‑TTTTATCCCCC‑TTT‑TCCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                             <  2:958987/149‑1 (MQ=255)                                                                                                                       cTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttt                        >  1:186332/1‑151 (MQ=255)                                                                                                                           cGCCATAAATAGG‑CCAGT‑TTTTCT‑CCCT‑TTT‑ATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                      <  2:259868/150‑1 (MQ=255)                                                                                                                              cATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                  >  2:599848/1‑151 (MQ=255)                                                                                                                              cATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                  >  1:1008781/1‑151 (MQ=255)                                                                                                                                 aaataaC‑CCAGG‑TTTTAT‑CCCT‑TTA‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaa                <  2:56289/149‑1 (MQ=255)                                                                                                                                   aaaaG‑CCGGT‑TTTTAC‑CCCT‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt             <  2:856841/149‑1 (MQ=255)                                                                                                                                    taaG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt             >  1:1072599/1‑150 (MQ=255)                                                                                                                                    aaaaaccctg‑TTTTAA‑CCCCTTTT‑ATCAACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa               <  2:109362/141‑1 (MQ=255)                                                                                                                                     aaG‑CCAGG‑GTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt             >  2:455857/1‑149 (MQ=255)                                                                                                                                         ccAGT‑TTTTAT‑CCCT‑TTT‑ATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGa        <  1:741157/151‑1 (MQ=255)                                                                                                                                         ccAGT‑TTTTAT‑CCCT‑TTT‑ATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg         <  1:421339/150‑1 (MQ=255)                                                                                                                                         ccAGT‑TTTTAT‑CCCC‑TTT‑ATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATa            <  1:853927/147‑1 (MQ=255)                                                                                                                                             g‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCAt    >  2:191299/1‑151 (MQ=255)                                                                                                                                              tttttAC‑CCCT‑TTT‑ATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCa     <  1:255379/149‑1 (MQ=255)                                                                                                                                                tttAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTg  >  1:700125/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                   AAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTT‑AAATTATTACC‑GGTGAAATCAAGTTCAAATACGG‑CATCCGTTCGCAATCCTGGTT‑TGG‑CTT‑CGCCATTAATAAG‑CCAGG‑TTTTAT‑CACC‑TTT‑ATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTG  >  NC_000913/4537968‑4538253

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑