Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	43.5%	88.5 / 14.9	62	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (8/27);  new base G (27/0);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.53e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508514‑3508806                                                                                                                                                  |                                                                                                                                                         aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GCGG‑gaaa                                                                                                                                                       >  1:720316/1‑147 (MQ=255)      gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTcacac                                                                                                                                               <  2:106847/150‑1 (MQ=255)          ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGG‑GCGG‑GAAAAATcaca                                                                                                                                                >  2:586530/1‑145 (MQ=255)            aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑g                                                                                                                                       <  1:539117/151‑1 (MQ=255)             aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑GAAAAATTAAACACGGt                                                                                                                                          >  1:1223189/1‑148 (MQ=255)              aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GGGG‑GAAAAATTACACACGGGA‑Gcc                                                                                                                                     >  1:1136513/1‑151 (MQ=255)              aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGG‑GCGG‑gaaaaataaaaaaacgga                                                                                                                                         >  1:670088/1‑135 (MQ=255)               gCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑Gccc                                                                                                                                    <  2:630950/151‑1 (MQ=255)                  ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCca                                                                                                                                 <  2:797985/151‑1 (MQ=255)                     gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACTGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑GAAAAATTACACACGGGA‑Gccccccc                                                                                                                                >  2:1051279/1‑147 (MQ=255)                     gagaTCAATATGCAGATGCTCATGGGCGAAAGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGCCGCACCGGAAAAAACTcacacagcatcccc                                                                                                                                      >  2:952270/1‑138 (MQ=255)                      agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGG‑GCGG‑GAAGAATcacttacagttagcccc                                                                                                                                   >  2:73884/1‑132 (MQ=255)                        aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGG‑GCGG‑GAAAAATCAAACCAAGTA‑Tccccccc                                                                                                                                >  2:670872/1‑144 (MQ=255)                         tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GGGG‑GAAAAATAACACACGGGA‑Gccccc                                                                                                                                  >  2:968811/1‑143 (MQ=255)                                gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑GAAAAATCACACACGGGA‑GCCCCCC‑CATCACAAAACgtg                                                                                                                  >  2:724233/1‑151 (MQ=255)                                gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑GAAAAATCAAACACGGAA‑ACCCCCC‑CCTCCAAAAGCgt                                                                                                                   >  2:825216/1‑150 (MQ=255)                                  aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGC‑GCGG‑GAAAAATTACACACGGTA‑GCCCCCA‑CATCAGAAAACgggg                                                                                                                 >  2:1113457/1‑150 (MQ=255)                                   gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCgtggtg                                                                                                               <  1:865023/151‑1 (MQ=255)                                    aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTgc                                                                                                              <  1:612638/151‑1 (MQ=255)                                      gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GCGG‑GAAAAATCAAACACCGGA‑GCCCCCC‑CAACACAAAACgggggggg                                                                                                             >  2:88434/1‑146 (MQ=255)                                            gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTACAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑GAAAAATCACACACCGTA‑GCCCCCA‑CAACAGAAAACGTGGGGCGCACAcc                                                                                                       >  2:1185537/1‑150 (MQ=255)                                               cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGc                                                                                                   <  2:95156/151‑1 (MQ=255)                                                    gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGt                                                                                              <  1:822896/151‑1 (MQ=255)                                                      gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                                                             <  2:445875/150‑1 (MQ=255)                                                        aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                                                                                           <  1:357272/150‑1 (MQ=255)                                                           ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                       <  1:836455/151‑1 (MQ=255)                                                           ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGG‑GCGG‑GAAAATACCCAAACGGCA‑GCCCCCC‑CACCAAAAAGGGAGGGGCGCCGACCCAGCTGCGTGGcccgc                                                                                       >  1:1068923/1‑147 (MQ=255)                                                             gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatcat                                                                                     <  1:330905/151‑1 (MQ=255)                                                               cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATcc                                                                                   <  2:735020/151‑1 (MQ=255)                                                                ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑GAAAAATCACACCCGGGA‑GCCCCCA‑CATCAGAAAGCGTGGGGCGGACACCGAGCTGCGGAGCGAACATccc                                                                                  >  2:45935/1‑151 (MQ=255)                                                                 gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GGGG‑GAAAAATCACACACCGTA‑GCCCCCC‑CACCAAAAAACGAGGGGGGCACACCGACCTGCGCGGCGATCATccc                                                                                  >  2:213148/1‑150 (MQ=255)                                                                     aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑GAAAAATTACACACGGGA‑GGCCCCA‑CATCAGAAAACGGGGGGGGCACACCGCGCCGCGGAGCGATCATCCCCTGCt                                                                             >  2:28805/1‑151 (MQ=255)                                                                      aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑GAAAAATTACACACGGGA‑GCCCCCC‑CCTCAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCtt                                                                                >  2:978363/1‑147 (MQ=255)                                                                       aaCTCATAGGGTGTCTCCCGCAAACCGGAGGGAGAAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑GAAAAATCAAAAACCGTA‑GCCCCCC‑CATCAAAAAACGGGGGGCGCACACCGAGCTGCGGAGCGATCATcccct                                                                                >  2:220750/1‑144 (MQ=255)                                                                           cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑aaa                                                                      <  2:998589/151‑1 (MQ=255)                                                                           cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑aa                                                                       <  1:352949/150‑1 (MQ=255)                                                                               gggTGTCGCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGGCCATTTTGTGGCAGTGACGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCa                                                                   <  2:648378/150‑1 (MQ=255)                                                                                 gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGTGCGG‑GCGC‑GAAAAATGAAACACCGGA‑CCCCCCACCACAAGAAATGGGGGT‑CGCCCACCGAACGGGGTAGCGATCATCCCTTGCTGGT‑aa                                                                       >  2:868660/1‑144 (MQ=255)                                                                                   gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACt                                                              <  2:1136573/151‑1 (MQ=255)                                                                                     ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GGGG‑GAAAAATTACACACGGGG‑ACCCCCA‑CATCAAAAAGCGAGGGGCGCACACCCAGCCGGGGAGCGGGCAACCCCCGCTGAT‑AAACCAGCACCgg                                                            >  2:636413/1‑151 (MQ=255)                                                                                       ccTGTAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTg                                                             <  1:708264/148‑1 (MQ=255)                                                                                           cAAACCGAAGGGAGGAAAGAAAAGAGCGACCATTTTGTGTCGGTGGCGC‑GCGG‑AAAAAGTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AACCCAGCACTGGcaccac                                                      <  1:42281/151‑1 (MQ=255)                                                                                            aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGG‑GCGG‑GAAAAATAAAACACCGGA‑GGCCCCA‑CAACAGAAAAAGTGGGGGGCACACCCAACTGCGTAGGGAACACACCTTGCTGAT‑AAAACAGCACTAGcgcc                                                        >  2:1227629/1‑148 (MQ=255)                                                                                            aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGC‑GCGG‑AAAAAATTACACAACGTA‑CACCCCC‑CACCCGAAAACGCGGGGCGCACACCCGGCCTCGTAGCGATCATCCCCTGCTGTT‑AAACAAACACTGGgacccc                                                      >  2:1110124/1‑148 (MQ=255)                                                                                              aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GCGG‑GAAAATACAAACACGGGA‑GGTTCCC‑CCACCGAAAGCGTGGGGCGCACAACGAGCCGCGGAGCCAGCACCCCTTTCTGTTAAAACCAGCCCGGCCACCCCt                                                     >  1:661228/1‑150 (MQ=255)                                                                                                   agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGC‑GGGG‑GAAAAATTACACACCGTATCCCCCCC‑CATCA‑AAAGCGGGGTGCGCACACCGAGCCTCGTAGCCATAACTCCCTCCTTCA‑AAAACAACCCCGGCACCACTTCGCggg                                              >  2:481854/1‑151 (MQ=255)                                                                                                    gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGt                                             <  1:658315/151‑1 (MQ=255)                                                                                                      gaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTAt                                           <  2:597576/151‑1 (MQ=255)                                                                                                           aagaaaagaGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTAtggtgg                                      <  2:162099/151‑1 (MQ=255)                                                                                                             gaaaagaGCGACCATTTTGTGGCGGTGGGGC‑GCGG‑GAAAAATCACACACCGGA‑GCCCCCA‑CATCAGAAAACGTGGGGCGCACACCCAGGTGCGTAGCGATCAACCCCTTCTGGT‑AAAACAGCACTGGCACCACCTCGCGGGGAtggtgg                                      >  2:106635/1‑149 (MQ=255)                                                                                                              aaaagaGCGACCATTTTGTGGCGGTGGGGG‑GCGG‑GAAAAATAAAAAACCGGA‑ACCCCCC‑CCTCAAAAAACGGGGGGGGCACCCCCAGCTGGGGAGCGATCAACCCCTTCTGAT‑AAAACAACACCGGCACCACCTCCCGGGGAtggggggt                                    >  1:1087291/1‑148 (MQ=255)                                                                                                                 agagCGACCATTTTGTGGCGGTGGGGC‑GCGG‑GAAAAATCAAACACGGGA‑GCCCCCA‑CCCCAGAAAACGGGGGGCGCCCACCGAGCCGCGGAGGGATAATCCCCTTCTGAT‑AAAACAGCACCGGGCCCACCTCGCGGGGATGGGGGCtg                                   >  2:1029898/1‑148 (MQ=255)                                                                                                                 agagCGACCATTTTGTGGCGGGGGTGC‑GCGG‑AAAAATTTCCACCCAGCT‑GCCCCCC‑CAATCTCAAGAGAGGGGCGCCCACCGAGCCGTGAAACGAGAAACCCCTTCTGCT‑AAAACAGCAATCGCACCACCTCCCGGGGAgggtgggtctgg                                >  1:366420/1‑145 (MQ=255)                                                                                                                  gagCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAAATTACACACGGGA‑GCCCCCA‑CATCAGAAAACGGGGGGGGCACACCCAGCCGCGGAGCGATCATCCCCTGCTGAT‑AAAACAGCACTGGCACCACCTCCCGGGGATGGGGGGtgtg                                 >  2:781200/1‑149 (MQ=255)                                                                                                                     cGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                             <  2:725118/150‑1 (MQ=255)                                                                                                                          aTTTTGTGGCGGTGGCGC‑GCGG‑GAAAAATTACACACGGCA‑AGCCCCC‑CCTCAGAAAAAGGGGGGGGCCCACCCCGCTGCGTAGCGGTCCTCCCTTTCTGGT‑AAACCAACACCGGCACCACCTCGCGGGGATGGGGGCTGTGGCCCAtag                         >  1:280762/1‑147 (MQ=255)                                                                                                                                tggcggtggcgC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                 >  1:380294/1‑151 (MQ=255)                                                                                                                                 ggcggtggcgC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                 <  1:724201/150‑1 (MQ=255)                                                                                                                                 ggcggtggcgC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                 <  2:200062/150‑1 (MQ=255)                                                                                                                                 ggcggtggcgC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt                >  1:117119/1‑151 (MQ=255)                                                                                                                                  gcggtggcaC‑GCGG‑AAAAATTCACACACGGTA‑ACCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATTCCTTGCTGAT‑AAACCACCACTGGCACCACGTCGCGGGTATGGTGGCTGTGTCCAATGGTCGGAt                   <  1:43371/147‑1 (MQ=255)                                                                                                                                             gcgG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATcc     >  1:925132/1‑151 (MQ=255)                                                                                                                                              cgG‑AAAAATCCACACACGGTA‑GCCCCCA‑CATCAGGCAGCGGGGTGCGCACACCGAGCTGCGCAGCGATCATCCCTTCCTGAT‑AACCCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGATCGGCTCGTTGCCGCGATcc     >  1:799707/1‑150 (MQ=255)                                                                                                                                                g‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc  <  2:662002/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                         ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508514‑3508806

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