Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,376,505	0	A	C	46.4%	87.2 / 12.7	57	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (23/7);  new base C (0/26);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 7.85e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAT  >  NC_000913/1376359‑1376650                                                                                                                                                        |                                                                                                                                                   tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCt                                                                                                                                                >  1:200927/1‑151 (MQ=255) tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCt                                                                                                                                                >  1:839677/1‑151 (MQ=255)     gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGcc                                                                                                                                            >  2:505368/1‑151 (MQ=255)        gATGACCTCGCCGCGTAGCCTGGCCGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GAt                                                                                                                                        <  2:257096/151‑1 (MQ=255)           gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGc                                                                                                                                      <  1:215872/150‑1 (MQ=255)             ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGccc                                                                                                                                    <  1:126444/150‑1 (MQ=255)                cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCccgcc                                                                                                                                 <  2:843552/150‑1 (MQ=255)                 gccgcGGAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGGCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGt                                                                                                                               >  2:806142/1‑151 (MQ=255)                      gTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACAt                                                                                                                          >  2:734857/1‑151 (MQ=255)                        aGCCTGGACGGCGACCGCAAAATGGAGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACAtcg                                                                                                                        >  2:483539/1‑151 (MQ=255)                             ggACGGCGACCGAAAAATGGTGAGTGATCCCCGGCGCGTACGGGAACAGGTCGCCGCGCGTCCCGGAACGGTGGTGTATCTCGACCCGTTAAC‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACCCCCTCGCC‑GATGCCCGCCGTTACATCGTCGaa                                                                                                                   <  2:303964/151‑1 (MQ=255)                                 ggCGACCGCAAAAGGGTGAGTGATCCCCTGCGCGTACGGGAACAGGTCGCCGCGCTTCCCGGAACGGTGGTGTATCCCGACCCGTTAAC‑CCGT‑T‑ATCAGG‑ATCATATGATC‑AACCCCCCCGCC‑GAGGCCCGCCGTTACATCGTCGAAAAc                                                                                                                <  1:825059/150‑1 (MQ=255)                                 ggCGACCGAAAAATGGTGGGTGTTTCCCTGCGCGTACGGGAACAGGCCGCCGCGCGTCCGGGACCGGTGGTGTATCTCGACCCGTAAAC‑CCGC‑T‑ATCAGG‑ATCATATGATA‑AACCCCCTCGCC‑GATGCCCGCCGTTACATCGTCGaaa                                                                                                                  <  2:403773/148‑1 (MQ=255)                                 ggCGACCGAAAAAGGGTGAGTGTTCCCCGGCGCGTAGGGGAACAGGTCGCCGCGCGTCCCGGAACCGTGGTGTATCTCGAACCGTAAAC‑CCGC‑T‑ATCAGG‑ATAAAATGATC‑AACCCCCTCGCC‑GATGCCCGCCGTTACAtc                                                                                                                         <  1:865936/141‑1 (MQ=255)                                   cGACCGCAAATTGGGGAGTGATTCCCTGGGCGTACTGGAACAGTCCCCCGCGCTTCCCGGACCGGTGGTGTTTCTC‑AACCCTTAAA‑ACGG‑TAACCAGA‑AAAACGTATTC‑AAC‑CCCCCGCCAGCGGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                             <  2:364338/151‑1 (MQ=255)                                         aaaaaGGGTGAGTGTTCCCCGGCGCGTACTGAAACAGGTCCCCGCGCGTCCCGGACCCGGGGTGTATCTCGACCCGTTAAC‑CCGC‑T‑ATCAGG‑ATCATATGATA‑AACCCCCCCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                           <  2:91525/146‑1 (MQ=255)                                          aaaaGGGTGAGTGATCCCCTGCGCGTACTGAAACAGGTCCCCGCGCGTCCCGGACCCGTGGTTTTTCTCGACCCTTAAAC‑CCGC‑T‑ATAAGG‑ATCATATAATC‑AACCCCCCCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                       <  2:671975/150‑1 (MQ=255)                                           aaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                      >  2:517376/1‑150 (MQ=255)                                                          cctgccgcgtcggggAGCAGGCCCCCGCGCGTGCGGGACCGGGGGTTTTTCCCG‑ACCGTTAAACCCGC‑T‑ATCAGG‑TTAATATGATC‑ACCCCCCGCCCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                      <  1:505079/138‑1 (MQ=255)                                                          ccccGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGACCGGTGGTGTATCTCGAACCGTTAAC‑CCGC‑T‑ATCAGG‑ATCATATGAAC‑ACCCCCCCCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACTGATTATCg                                                                                       <  1:56799/150‑1 (MQ=255)                                                          cccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                      >  2:97562/1‑151 (MQ=255)                                                          cccGGCGCGTCCGGGACCAGGCCGCCGCGCGTCCCGAACCCGTGGTGTATCTCGAACCGTTAAC‑CCGC‑T‑ATCAGA‑ATCATATGACA‑AACCCCCTCCCC‑GATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                        <  1:264089/149‑1 (MQ=255)                                                           ccggggcttccTGAAACAGGTCGCCGCGCGTACCGGAACGGTGGTTTATCCCGACCTTTTAAC‑CCCC‑T‑ATCAGA‑ATCATATGATC‑AACCCCCCCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                         <  1:651414/137‑1 (MQ=255)                                                                gcgTACTGGAACGGGTCGCCGCGCGTCCCGGACCGGTGGTTTACCTCG‑ACCCTTAAA‑CCCCTT‑ATCAGG‑ATAATATAATC‑AACCCCCCCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                                 <  2:552526/150‑1 (MQ=255)                                                                       ggAAGAGGTCGCCGCGCTTACGGGAACGGTGGTTTCTCCCGACCCTTTAAC‑CCGC‑T‑A‑CAGGTTTCATAGAACC‑AACCCCCCCGCG‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                           <  2:384168/149‑1 (MQ=255)                                                                         aacgggccccggcgcgTACGGAACCCGGGGTGTACCCCAACCCTTTAAC‑CCGC‑T‑ACCGGG‑ATAAAAT‑ATCAAACCCCCCCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                            <  1:95812/135‑1 (MQ=255)                                                                         aaCAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                       >  1:814056/1‑151 (MQ=255)                                                                         aaCAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                       >  2:31816/1‑151 (MQ=255)                                                                           cAGGTCCCCGCGCGTACCGAAACCGTGGTGTATCCCGAACCTTTAAC‑CCCC‑T‑ATCAGA‑AAAATATAATC‑ACCCCCCCCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                      <  1:388784/150‑1 (MQ=255)                                                                              gTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                   >  1:306240/1‑150 (MQ=255)                                                                                cgccccGCGTCCGGGACCGGTGGTGTATCTCGACCCTTAAAC‑CCTT‑T‑AACAGG‑ATAATATGATC‑AACCCCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                  <  1:817894/149‑1 (MQ=255)                                                                                 gccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                               >  1:27209/1‑151 (MQ=255)                                                                                  ccggGCGTACCGGACCGGTGTTTTCTCTCGAACCGTAAAC‑CCCC‑T‑ATCAGG‑ATAATATAACC‑ACCCCCCCCGCG‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                   <  2:710043/146‑1 (MQ=255)                                                                                  ccgcGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAC‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                    <  2:734460/145‑1 (MQ=255)                                                                                         tACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                           >  1:511406/1‑147 (MQ=255)                                                                                          aCCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                       >  1:699150/1‑150 (MQ=255)                                                                                             ggAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                   >  2:636278/1‑151 (MQ=255)                                                                                             gaacccGTGGTGTATCTCGAACCGTTAAC‑CCGC‑T‑ATCAGA‑ATCATATGATC‑AACACCCCCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                     <  2:914004/145‑1 (MQ=255)                                                                                                           tctcGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                     >  1:206782/1‑151 (MQ=255)                                                                                                            cccGAACCGTTAAC‑CCGC‑T‑ATCAGG‑ATCATATAACC‑AACCCCCCCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                      <  2:161333/147‑1 (MQ=255)                                                                                                               gAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                 >  2:766214/1‑151 (MQ=255)                                                                                                                aaCCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                 >  1:257365/1‑150 (MQ=255)                                                                                                                  ccGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                               >  2:307026/1‑150 (MQ=255)                                                                                                                           ccGC‑T‑ATCAGA‑ATCATATGATC‑ACCCCCCCCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                        <  2:856351/149‑1 (MQ=255)                                                                                                                                t‑ataAGG‑ATAATATAATA‑AACCCCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                    <  2:435823/149‑1 (MQ=255)                                                                                                                                    cAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg               >  2:844402/1‑151 (MQ=255)                                                                                                                                     aGG‑ATATTATTACC‑AACCCCCCCGCC‑GGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                <  1:324904/149‑1 (MQ=255)                                                                                                                                     aGG‑ATAATTTGATA‑AACACCCTCGCC‑GATCGCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt                 <  2:39107/148‑1 (MQ=255)                                                                                                                                       gattCATATAACC‑AACCCCCCCGCC‑GATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt                 <  1:2873/144‑1 (MQ=255)                                                                                                                                            atatGACC‑AACCCCCCCGCC‑GTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg         <  1:266650/150‑1 (MQ=255)                                                                                                                                             tatGACC‑AACCCCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc            <  2:840907/146‑1 (MQ=255)                                                                                                                                              atGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc      >  2:264497/1‑151 (MQ=255)                                                                                                                                              atGACC‑AACCCCCACGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt          <  1:223326/147‑1 (MQ=255)                                                                                                                                               tGATC‑ACCCCCCGCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc      <  2:906182/150‑1 (MQ=255)                                                                                                                                                tgacaaaCCCCCGCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt          <  2:908442/141‑1 (MQ=255)                                                                                                                                                 aTC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat    >  2:556456/1‑150 (MQ=255)                                                                                                                                                  tC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat  >  2:139692/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                   TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAA‑CCGC‑T‑ATCAGG‑ATCATATGATC‑AACACCCTCGCC‑GATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAT  >  NC_000913/1376359‑1376650

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