Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,662	0	T	A	45.0%	80.5 / 28.0	62	E384D (GAA→GAT)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (7/26);  new base A (27/0);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 8.34e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508519‑3508806                                                                                                                                                  |                                                                                                                                                    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)          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)                                                                                                                                                  |                                                                                                                                                    GTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGG‑AAAAATTCACACACGGTA‑GCCCCCA‑CATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAT‑AAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508519‑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: ‑