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

CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCT  >  NC_000913/1376355‑1376635                                                                                                                                                          |                                                                                                                                    cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAaca                                                                                                                                     <  1:967367/151‑1 (MQ=255)  gCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAaca‑c                                                                                                                                   >  2:291272/1‑151 (MQ=255)     tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑ccc                                                                                                                                 <  1:1072460/150‑1 (MQ=255)      aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTc                                                                                                                               <  2:238010/151‑1 (MQ=255)       ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTc                                                                                                                               <  2:629354/150‑1 (MQ=255)         gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCTTCGcc                                                                                                                            >  1:755465/1‑151 (MQ=255)            gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGa                                                                                                                          >  2:289501/1‑150 (MQ=255)            gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGAt                                                                                                                         <  1:262936/151‑1 (MQ=255)            gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GCGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGAt                                                                                                                         >  2:663340/1‑151 (MQ=255)            gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGAGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGAt                                                                                                                         <  1:1103948/151‑1 (MQ=255)                aCCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGccc                                                                                                                     >  2:1104736/1‑151 (MQ=255)                 ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCccg                                                                                                                    <  1:952278/151‑1 (MQ=255)                    cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCccgccg                                                                                                                 >  2:732934/1‑151 (MQ=255)                    cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCccgcc                                                                                                                  <  2:753204/150‑1 (MQ=255)                    cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCccgc                                                                                                                   <  2:166807/149‑1 (MQ=255)                     gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGt                                                                                                                >  1:780850/1‑151 (MQ=255)                            aGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                         >  1:880622/1‑151 (MQ=255)                            aGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                         >  1:778912/1‑151 (MQ=255)                                gggggggcgcccaaaaagTGGTGAGTAATCCCCGGCGCGTACGGGAACGGGCCGCCGCGGGTCCCGGA‑CCC‑GGGGTGTATCCGGAACCGTAAAC‑CCGCTATAAGGATCATATTATAACCC‑CCCTCGCCTATGCCCGCCGTTACATCGTCGa                                                                                                     <  2:839064/137‑1 (MQ=255)                                    cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                 >  2:99971/1‑151 (MQ=255)                                      gcgcCCGAAAAAGGTTGGGGGATCCCCTGCGCGTGCGGGAACAGGCCGCCGCGCGTA‑CGGA‑ACCGGTGTTGTTTCTCGACCCTTAAAC‑CCGCTATAAGAATAATATGATCACCC‑CCCCCGCCGATGCCCGCCGTTACATCGTCGaaa                                                                                                   <  2:970669/147‑1 (MQ=255)                                       cgcccaaaaaaTGGTGAGTTTTTCCCTGCGCGTAGTGGAACGGGTCGCCGCGCGTCCCGGA‑CCC‑GTGGTTTATCCCAACCCTTTAAC‑CCGCTATCAGTATCATATGATCAACC‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                               <  1:507278/143‑1 (MQ=255)                                         aCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                            >  1:461943/1‑151 (MQ=255)                                          ccGAAAAGTGGTGGGTGTTCCCCTGCGCGTATTGGAACGGGTCCCCGCGCGTCCCGGA‑CCC‑GTGGTGTATCTCAACCCGTAAAC‑CCCTTATCGGGATAATATGATCACCC‑CCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                           <  1:706730/151‑1 (MQ=255)                                              aaaaTGGTGAGTGATCCCCTGCGCGTACTGGAACGGGCCGCCGCGGGTCCCGGA‑CCC‑GTGGTTTTTCTCAACCCTTTAAC‑CCGCTATCAGGATCATATGACCAACC‑CCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                        <  1:205234/150‑1 (MQ=255)                                              aaaaGGGGGGGTGTTCCCCTGCGGGTACGGAAACGGGTCGCCGCGCGTCCGGGA‑CCC‑GTGGTTTATCCCGACCCTTAAAC‑CCGCTATCAGGATAATATGACAACCC‑CCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                        <  1:909114/150‑1 (MQ=255)                                               aaggggTGGGTTATCCCCGGCGCTTACGGGAACAGCGCCCCGCGCTTCCCGAG‑CCC‑GTGGTGTACCTCGAACCATTAAC‑CCGTTATCAGGATAATATAAACACCC‑CCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                      <  2:850018/147‑1 (MQ=255)                                                aaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                     >  2:342577/1‑151 (MQ=255)                                                    gtgagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                 >  2:697801/1‑151 (MQ=255)                                                         tgtttCCCTGCGCGTACGGGACGGGGTCGCCGCGCTTACCGGA‑CCC‑GGGGTGTATCCCGACCCTTAAAC‑CCGTTATCAGGATAAAATGATCACCC‑CCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                             <  2:557346/147‑1 (MQ=255)                                                              cccggcgtttACGGGAACAGGTCCCCGCGCTTCCCGGA‑CCG‑GTGTTTTACCTCGACCCTTAAAC‑CCGTTACCAGGATAATATGATCAACC‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                        <  2:410019/141‑1 (MQ=255)                                                              cccTGCGCGTTCGGGAACAGGTCGCCGCGCGTCCCGGA‑CCC‑GTGGTGTATCCCGACCCTTTAAC‑CCGATATAAGTATAATATAAACACCC‑CCCGCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                       <  1:214010/151‑1 (MQ=255)                                                              cccGGCGCGTAGTGAAACAGGTCCCCGCGCGCACCGAA‑CCC‑GTGGTGTATCTCGACCCTTTAAC‑CCGTTACCAGGTTCATATGACAACCC‑CCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                          <  2:450116/148‑1 (MQ=255)                                                               ccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                      >  1:567865/1‑151 (MQ=255)                                                                        aCTGAAACGGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGACCCTTAAAC‑CCGCTATCAGGATCATATGATCAACC‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                              <  1:1017190/150‑1 (MQ=255)                                                                         cTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                            >  2:595316/1‑151 (MQ=255)                                                                           ggAACGGGGCGCCGCGCGTACCGGA‑CCG‑GTGGTTTACCTC‑AACCCTTAAACCCTTTATCAGAATCAAATGATCAACC‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCt                                                               <  2:967827/146‑1 (MQ=255)                                                                                gggcccccgcGCGT‑CCGGACACC‑GTGGTGTATCTCGACCCGTAAAC‑CCGCTATAAGAATAATATGACAAACA‑CCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                     <  2:963219/145‑1 (MQ=255)                                                                                 ggTCGCCGCGCTTCCCGGA‑CCC‑GTGGTGTATCTCGACCCTTTAAC‑CCGCTATCAGGATAATATGATAAACC‑CCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                     <  1:251114/150‑1 (MQ=255)                                                                                    cgccgcGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                  >  2:300622/1‑150 (MQ=255)                                                                                    cccggcgcgTACCGGA‑CCG‑GTGGTTTTTCTCGAACCGTTAAC‑CCGCTAAAAGGATAATATAATCAACC‑CCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                  <  2:611886/146‑1 (MQ=255)                                                                                      ccgcGCGTACCGAA‑CCG‑GTGGTGTATCTCGAACCGTAAAC‑CCGCTATCAGGATCATATAATAACCC‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                               <  2:380882/151‑1 (MQ=255)                                                                                          gctgcccGGA‑CCG‑GTGGTTTACCCCGAACCGTAAAC‑CCGCTATAAGGATCAAATGATCACCC‑CCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                           <  2:1024424/146‑1 (MQ=255)                                                                                               ccGGA‑CCC‑GTGGTTTATCTCGACCCTTAAAC‑CCGCTATCAGGATCATATAATCAACC‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                      <  1:785664/151‑1 (MQ=255)                                                                                                 ggA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                    >  2:1041677/1‑151 (MQ=255)                                                                                                    accc‑GTGGTGTATCCCGACCCTTAAAC‑CCGCTATCAGGATAATATGACCACCC‑CCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                     <  1:599196/146‑1 (MQ=255)                                                                                                        ctgggggtttcacGAAGCTTTAAC‑CCGCTATAGGAATAATA‑GGTCAACCCCCCTCGCCAATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                    <  1:1101700/134‑1 (MQ=255)                                                                                                                  cgcaacccGTTAAC‑CCGTTATGATGATCATTTAATCAACC‑CCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                          <  1:970800/140‑1 (MQ=255)                                                                                                                    cGACCCGTTAAC‑CCGCTATCAGGATCATTTCATCAACC‑CCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                     <  1:525586/149‑1 (MQ=255)                                                                                                                    cGACCCGTAAAC‑CCGCTATCAGGATCATATGACCAACC‑CCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                      <  1:291328/148‑1 (MQ=255)                                                                                                                     gAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                  >  2:1053172/1‑151 (MQ=255)                                                                                                                       aCCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                >  2:263418/1‑151 (MQ=255)                                                                                                                        ccGTAAAC‑CCGCTATCAGGATCATATGATCAACC‑CCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                <  1:958178/150‑1 (MQ=255)                                                                                                                          gTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc             >  1:696502/1‑151 (MQ=255)                                                                                                                           ttAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc            >  1:293758/1‑151 (MQ=255)                                                                                                                              aacccGCTTTAATGATCATATGACCAACC‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc          <  1:707853/148‑1 (MQ=255)                                                                                                                                 ccGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACAGCGa        >  2:386254/1‑150 (MQ=255)                                                                                                                                    cTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg    >  1:312052/1‑151 (MQ=255)                                                                                                                                      aTCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt  >  1:411902/1‑151 (MQ=255)                                                                                                                                         aGGATCATATGACCAACC‑CCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc      <  1:408730/144‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                    CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGA‑ACC‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACA‑CCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCT  >  NC_000913/1376355‑1376635

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: ‑