Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,538,109	0	A	C	47.0%	92.2 / 23.2	66	G219G (GGT→GGG)	nanM	N‑acetylneuraminic acid mutarotase
Reads supporting (aligned to +/- strand):  ref base A (32/3);  new base C (0/31);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 9.39e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATT  >  NC_000913/4537962‑4538256                                                                                                                                                     |                                                                                                                                                    ccccGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCt                                                                                                                                                  >  1:480114/1‑151 (MQ=255)      cAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTAtt                                                                                                                                             >  2:339329/1‑151 (MQ=255)        aaaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTc                                                                                                                                            >  2:781775/1‑150 (MQ=255)         aaCCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTcac                                                                                                                                          >  2:500984/1‑151 (MQ=255)                     cGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACcgccgc                                                                                                                                 >  1:562681/1‑148 (MQ=255)                       ccATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCAcc                                                                                                                              >  1:645610/1‑149 (MQ=255)                          tCTGGTGATGAGACGGGAGCAAGTTTATCCCTTTTTAAATTATTCCCGGTAAAATCAAGTTCAAATCCGGCATCCTTCCGCAATCCTGTTTTGG‑CTTCCCCATTAATAAGCAAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGcc                                                                                                                          <  1:381413/150‑1 (MQ=255)                           cTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGt                                                                                                                        >  2:829917/1‑151 (MQ=255)                           cTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGAGAAATCAAGCTCAACTACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCg                                                                                                                         >  2:324416/1‑150 (MQ=255)                              gTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTcc                                                                                                                     >  2:858487/1‑151 (MQ=255)                                   gagaCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTAcc                                                                                                                >  2:259323/1‑151 (MQ=255)                                      aCGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAgg                                                                                                             >  2:32090/1‑151 (MQ=255)                                              aaGCTTATTCCATTTTAAATTATTCCGGGTAAAATCAAGTTCAAATACGGCTCCCGTTCGCATCCCTGTTTGGG‑CTTCGCCATAAAAAAGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGAtt                                                                                                       <  1:288893/149‑1 (MQ=255)                                                  ttATTCCATTTAAAATTATTCCCGGTAAAATCAATTTAAAATACGCCATCCTTCCCCAACCCTGTTTGGG‑TTCCGCCATAAAAAAGCCAGGTTTTACCCCCTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCa                                                                                                  <  1:579831/150‑1 (MQ=255)                                                     ttCCTTTTAAATTTATTCCCGGTAAAATCAATTTAAAATCCGGCCTCCTTTCCAAACCCTGGTTGGG‑TTTCCCTTTAAAAAAGCCAGTTTTTATCCCTTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg                                                                                               <  1:950953/150‑1 (MQ=255)                                                         aTTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                          >  1:259719/1‑151 (MQ=255)                                                         aTTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt                                                                                          >  1:495180/1‑151 (MQ=255)                                                                  tattaCCGGTAAAATCAATTTCAAATCCGCCACCCTTTCGCAATCCTGTTTTGG‑CTTCGCCATTAATAAGCCAGTTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCt                                                                                  <  2:220781/150‑1 (MQ=255)                                                                     taCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTgtgt                                                                              >  2:766382/1‑151 (MQ=255)                                                                             aaataaattaaaaaTCCGGACTCCGTTCCCATCCCTGGTTTGG‑CTCCCCAATTAATAAGCCGGTTTTTACCCCTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                        <  1:26547/138‑1 (MQ=255)                                                                             aaataaatTTAAAATCCGCCCTCCGTTCCCAACCCTGTTTTGG‑CTTCGCCTTAAAAAAGCCATGTTTTATCCCCTTTATCCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGAt                                                                      <  1:309564/151‑1 (MQ=255)                                                                             aaaTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGTTTTGG‑CTTCGCCATTAATAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAggg                                                                        <  2:168104/149‑1 (MQ=255)                                                                                tCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCa                                                                    >  1:843418/1‑150 (MQ=255)                                                                                 cAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa                                                                   >  2:547221/1‑150 (MQ=255)                                                                                  aagtaaaaaTACGGCACCGTTTCGC‑ACCCTGGTTTGGCCTTCGCTTTTAATACACCCGGTTTTCCCCCTTTTATCCAAAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaaa                                                                 <  2:921584/145‑1 (MQ=255)                                                                                  aaGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaaa                                                                  >  1:283407/1‑150 (MQ=255)                                                                                        aaaTCCGGCACCCTTTCGCAATCCTGTTTTGG‑CTTCGCCATAAATAACCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaa                                                            <  2:179547/150‑1 (MQ=255)                                                                                        aaaTACGGCACCCTTTCGAAATCCTGTTTTGG‑CTTCCCCATTAATAACCCAGTTTTTATCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGac                                                              <  2:599403/148‑1 (MQ=255)                                                                                          aTACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                           >  2:131226/1‑149 (MQ=255)                                                                                           tACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                           >  2:437833/1‑148 (MQ=255)                                                                                           tACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGacaac                                                           >  1:437520/1‑148 (MQ=255)                                                                                            aCGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGa                                                        >  1:120689/1‑150 (MQ=255)                                                                                             cGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAg                                                         >  1:30522/1‑148 (MQ=255)                                                                                             cGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaa                                                       >  2:208990/1‑150 (MQ=255)                                                                                              ggCATCCGTTCGCAACCCGGGTTTGG‑CTTCGCCATTAATAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGaaa                                                      <  2:482996/150‑1 (MQ=255)                                                                                                   ccGTTCGCAATCCTGTTTTGG‑TTTCGCCATAAATAAGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACtt                                                   <  2:274555/148‑1 (MQ=255)                                                                                                   ccGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTAtt                                                >  1:72137/1‑151 (MQ=255)                                                                                                        cGCAATCCTGTTTTGG‑CTTCGCCATTAAAAAGCCAGGTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTgaaga                                           <  2:271780/151‑1 (MQ=255)                                                                                                           aaTCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                        >  1:95090/1‑151 (MQ=255)                                                                                                              ccTGTTTTGG‑CTCCGCCATAAATAACCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                     <  2:932641/151‑1 (MQ=255)                                                                                                              ccTGTTTTGG‑CTCCCCCTTTAATAACCCAGTTTTTATCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaata                                       <  1:440861/149‑1 (MQ=255)                                                                                                              ccTGTTTGGG‑CTCCCCCATAAAAAGCCCAGTTTTTACCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaataat                                     <  1:260208/151‑1 (MQ=255)                                                                                                                 gttttGG‑TTCCGCCATTAATAAGCCAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtctt                                  <  1:707179/149‑1 (MQ=255)                                                                                                                 gttttGG‑CTTCGCCATTAATAAGCCAGGTTTTATCCCCTTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtct                                   <  2:598558/148‑1 (MQ=255)                                                                                                                  ttttGC‑CTTCCCCATTAATAAGCCAGTTTTTATCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAaat                                        <  1:778485/143‑1 (MQ=255)                                                                                                                   tttGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAAtcttct                                >  1:772144/1‑151 (MQ=255)                                                                                                                      gg‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCt                             >  2:913648/1‑151 (MQ=255)                                                                                                                       gtcTTCCCCATAAAAAACCCGGTTTTTATCCCCTTTATCCCCACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                              <  1:238383/148‑1 (MQ=255)                                                                                                                          ttCGCCATTAATAAGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCtt                            <  2:246700/149‑1 (MQ=255)                                                                                                                            ccccATAAATAAGCCAGTTTTTACCCCCTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGc                              <  1:724556/143‑1 (MQ=255)                                                                                                                            cGCAATTAATAAGCCAGTTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCttttt                         <  1:66237/150‑1 (MQ=255)                                                                                                                              cctttAAAAAGCCAGTTTTTATCCCCTTTTTCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                       <  2:500904/147‑1 (MQ=255)                                                                                                                              ccATTAATAAGCCAGTTTTTACCCCTTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                       <  1:905670/150‑1 (MQ=255)                                                                                                                                aTTAATAAGCAAGGTTTTTGCCCATTAATACACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                    <  2:562736/151‑1 (MQ=255)                                                                                                                                 ttAATAAGCCAGTTTTTACCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTc                     <  2:480161/149‑1 (MQ=255)                                                                                                                                 taaataaGCCAGTTTTTATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCa                    <  2:674018/148‑1 (MQ=255)                                                                                                                                   aataaGCCGGTTTTTACCCCCTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa                 <  1:385947/151‑1 (MQ=255)                                                                                                                                   aataaGCCAGTTTTTATCCCCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaa                  <  2:632396/150‑1 (MQ=255)                                                                                                                                   aaaaacccggtttttATCCCTTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCaaaa                 <  1:777389/140‑1 (MQ=255)                                                                                                                                   aaaaacccAGTTTTTATCCCTTTTATCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTg                       <  1:500935/139‑1 (MQ=255)                                                                                                                                     taaGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAAt                >  2:682233/1‑150 (MQ=255)                                                                                                                                         ccgttttttATCCCCTTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGt             <  1:688952/144‑1 (MQ=255)                                                                                                                                          cAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTg            >  2:252902/1‑149 (MQ=255)                                                                                                                                                 tATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAt   >  2:487715/1‑151 (MQ=255)                                                                                                                                                  aTCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtt  >  2:490407/1‑151 (MQ=255)                                                                                                                                                  aTCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGAtt  >  1:582603/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCCCGCAAAACCGCCAGCTACGCCATCTGGTGATGAGACGGGAGCAAGCTTATTCCATTTTAAATTATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGGTTTTATCACCTTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGACAACAGAAACTTATTGAAGAAATAATCTTCTGCTTTTTTGTCAAAATAGTGAGCATTGATT  >  NC_000913/4537962‑4538256

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