Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,620,676	C→T	6.0%	D358D (GAC→GAT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,676	0	C	T	6.0%	210.8 / 2.5	67	D358D (GAC→GAT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (39/24);  new base T (2/2);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

CCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTT  >  NC_000913/4620535‑4620822                                                                                                                                              |                                                                                                                                                   ccGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCTATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACAcgc                                                                                                                                             <  6:131815/149‑1 (MQ=255)      gcgACGTCACCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGaa                                                                                                                                        <  1:96681/149‑1 (MQ=255)        gACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGAGGACGACCTCCTGATCCTCACCGCTGACCACGGTTCCGATACGACCTGGACCGGTACTGACCACACGCGTGAaca                                                                                                                                      >  3:48007/1‑149 (MQ=255)          cGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAAcaca                                                                                                                                    <  7:94680/149‑1 (MQ=255)               ccGGTTATGCCGCGGTTCTGGAACTCTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCg                                                                                                                               >  1:87191/1‑149 (MQ=255)               ccGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGCCATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTcag                                                                                                                               >  1:54401/1‑147 (MQ=255)                cGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCgg                                                                                                                              >  3:224639/1‑149 (MQ=255)                      tGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTgg                                                                                                                        <  8:83146/149‑1 (MQ=255)                      tGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTg                                                                                                                         <  7:52507/148‑1 (MQ=255)                         cgcgGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGta                                                                                                                      >  4:55101/1‑148 (MQ=255)                                   gAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGa                                                                                                            >  6:53721/1‑148 (MQ=255)                                    aaCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGATCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGaaa                                                                                                          <  6:11873/149‑1 (MQ=255)                                    aaCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGaaa                                                                                                          <  6:136365/149‑1 (MQ=255)                                         ttCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTaa                                                                                                      >  7:209020/1‑148 (MQ=255)                                           cGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGCCCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAAc                                                                                                   >  5:141794/1‑149 (MQ=255)                                             aCCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGCCCACACGCGCGAACACATTCCGGTACTGGtat                                                                                                                     <  5:147907/129‑1 (MQ=255)                                             aCCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGCGAACACATTCCGGTACTGGtat                                                                                                                     >  6:147907/1‑129 (MQ=255)                                                             gCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATc                                                                                 <  2:87191/149‑1 (MQ=255)                                                               tGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGt                                                                               >  1:232040/1‑149 (MQ=255)                                                                   gTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGTCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGaaa                                                                           >  4:32147/1‑149 (MQ=255)                                                                      tctGCTGCGCGATGACGACATCCTGATCCTCACCGCTGCCCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGGATATGGCCCGAAAGTAAACCCGGGCTCACTGGGTCATCGTGAAACCt                                                                        >  1:180477/1‑149 (MQ=255)                                                                      tctGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCt                                                                        >  7:181476/1‑149 (MQ=255)                                                                           tgcGCGATGACGACATCCTGATCCTCACCGCTGACAACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTcgcg                                                                   <  5:53721/149‑1 (MQ=255)                                                                             cgcgATGACGACATCCTGATCCTCACCGCTGACCCCGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGTCCCGAAAGTAAAAACGGGCTCACTGGGTCATCGTGAAACCTTCGCGGa                                                                 >  6:179396/1‑149 (MQ=255)                                                                             cgcgATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTcgcg                                                                   <  1:233105/147‑1 (MQ=255)                                                                               cgATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGata                                                               >  3:138234/1‑149 (MQ=255)                                                                                     cgacATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTggg                                                                                      <  1:159045/120‑1 (MQ=255)                                                                                     cgacATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTggg                                                                                      >  2:159045/1‑120 (MQ=255)                                                                                      gccATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGc                                                          >  4:8945/3‑147 (MQ=255)                                                                                      gacATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGc                                                          >  6:226756/1‑147 (MQ=255)                                                                                        cATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGa                                                      >  2:24787/1‑149 (MQ=255)                                                                                         aTCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGAc                                                     >  2:91435/1‑149 (MQ=255)                                                                                           ccTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGActc                                                   >  8:20523/1‑149 (MQ=255)                                                                                            cTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGActct                                                  >  4:123227/1‑149 (MQ=255)                                                                                             tGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAacac                                                                                                                                     >  1:113888/1‑65 (MQ=255)                                                                                             tGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAacac                                                                                                                                     <  2:113888/65‑1 (MQ=255)                                                                                                 ccTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCa                                              >  5:46440/1‑148 (MQ=255)                                                                                                   tCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGaaa                                                                                                          <  3:173305/86‑1 (MQ=255)                                                                                                   tCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGaaa                                                                                                          >  4:173305/1‑86 (MQ=255)                                                                                                   tCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCaaaa                                           >  8:192787/1‑149 (MQ=255)                                                                                                    cACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCaaaa                                           >  1:75536/1‑148 (MQ=255)                                                                                                      ccGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAatat                                        >  4:125232/1‑149 (MQ=255)                                                                                                         cTGAGCACGGTTGCGATCAGACCTGGACCGGTACTGTCCACACGCGTGAACACTTTCCGGTACTGGTATATGGCCCGAAAGTAATACCGGGATCACTGGGTCATCTTGAAACCTTCGCGGATATCGGGCAGACTCTGGCAAAATAtttt                                     <  2:54401/149‑1 (MQ=255)                                                                                                         cTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATAtttt                                     <  6:79734/149‑1 (MQ=255)                                                                                                            aCCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTgg                                   <  3:55101/148‑1 (MQ=255)                                                                                                            aCCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTg                                    <  8:181476/147‑1 (MQ=255)                                                                                                             ccACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCAGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTACGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGt                                  >  4:45749/1‑148 (MQ=255)                                                                                                                 ggTTGCGATCCGACCTGGACCGGTACTGACCACACCCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTc                             <  4:48007/149‑1 (MQ=255)                                                                                                                  gTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCCCTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACttat                            >  1:183935/1‑147 (MQ=255)                                                                                                                       gATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGatat                       <  1:91435/149‑1 (MQ=255)                                                                                                                       gATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGatat                       >  1:221140/1‑149 (MQ=255)                                                                                                                        aTCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATg                      >  8:73959/1‑149 (MQ=255)                                                                                                                         tCCGACCTGGACCGGTACTGATCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCaaaa                                           >  8:58415/1‑127 (MQ=255)                                                                                                                         tCCGACCTGGACCGGTACTGATCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCaaaa                                           >  8:54896/1‑127 (MQ=255)                                                                                                                         tCCGACCTGGACCGGTACTGATCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCaaaa                                           <  7:54896/127‑1 (MQ=255)                                                                                                                         tCCGACCTGGACCGGTACTGATCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCaaaa                                           <  7:58415/127‑1 (MQ=255)                                                                                                                           cGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGa                    <  8:209020/148‑1 (MQ=255)                                                                                                                              ccTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAatat                >  3:64958/1‑149 (MQ=255)                                                                                                                               cTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCg                                                                                <  8:108255/84‑1 (MQ=255)                                                                                                                               cTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCg                                                                                >  7:108255/1‑84 (MQ=255)                                                                                                                               cTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATg               >  2:1285/1‑149 (MQ=255)                                                                                                                                 ggACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATgg              >  7:4791/1‑148 (MQ=255)                                                                                                                                 ggACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAACATGGc             >  2:102826/1‑149 (MQ=255)                                                                                                                                    ccGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGGATATGGCa            >  2:228756/1‑147 (MQ=255)                                                                                                                                          cTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCAt     <  2:232040/148‑1 (MQ=255)                                                                                                                                            gACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGtt  >  8:32461/1‑149 (MQ=255)                                                                                                                                             aCCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGtt  <  5:179396/148‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   CCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTT  >  NC_000913/4620535‑4620822

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