Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	521,180	A→T	6.0%	A255A (GCA→GCT)	ybbP →	putative ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	521,180	0	A	T	6.0%	212.7 / 3.7	67	A255A (GCA→GCT)	ybbP	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base T (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAGTGGTTGATGGTGCTGACGCTTTCA  >  NC_000913/521042‑521321                                                                                                                                           |                                                                                                                                              tGAGAACTGGTTGTTACCTCAGCTTAAACCCGAACAAGGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCTACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCg                                                                                                                                      <  4:40390/148‑1 (MQ=255) tGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCgg                                                                                                                                     >  3:275912/1‑149 (MQ=255)  gagaAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGt                                                                                                                                    >  1:282293/1‑149 (MQ=255)   agaAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCtggcagtgg                                                                                                                                           >  7:226567/1‑141 (MQ=255)   agaAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCtggcagtgg                                                                                                                                           <  8:226567/141‑1 (MQ=255)       aTGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCgg                                                                                                                               >  4:78385/1‑149 (MQ=255)              ttACCTCAGCTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCGCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGAt                                                                                                                        <  3:88510/149‑1 (MQ=255)                aCCTCAGCTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATg                                                                                                                       >  6:155176/1‑148 (MQ=255)                 ccTCAGCTTAAACCCGAACAACGCTGGTACGGTCGGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGaa                                                                                                                     <  6:41811/149‑1 (MQ=255)                  cTCAGCTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAAt                                                                                                                    >  4:174639/1‑149 (MQ=255)                  cTCAGCTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAAt                                                                                                                    >  6:359829/1‑149 (MQ=255)                       cTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTa                                                                                                               >  6:222573/1‑149 (MQ=255)                        ttAAACCCGCACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCAGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTAc                                                                                                              >  5:312651/1‑149 (MQ=255)                           aaCCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGt                                                                                                           <  5:359829/149‑1 (MQ=255)                            aCCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTc                                                                                                          <  2:282293/149‑1 (MQ=255)                            aCCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTGGCTGCTGGCAGTGGCAGCgg                                                                                                                                     <  1:342165/122‑1 (MQ=255)                            aCCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTAGCTGCTGGCAGTGGCAGCgg                                                                                                                                     >  2:342165/1‑122 (MQ=255)                             cccGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCg                                                                                                         <  8:361816/149‑1 (MQ=255)                              ccGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCTGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCg                                                                                                         >  3:163494/1‑148 (MQ=255)                              ccGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCTGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCg                                                                                                         >  3:157930/1‑148 (MQ=255)                                gAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCg                                                                                                                                <  8:222683/123‑1 (MQ=255)                                gAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCg                                                                                                                                >  7:222683/1‑123 (MQ=255)                                    aacGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGc                                                                                                  >  2:291602/1‑149 (MQ=255)                                         tGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCTCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATTAATCATTACTGTCTCAGTCGCTACGa                                                                                             <  7:127618/149‑1 (MQ=255)                                              cGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCtgg                                                                                        <  1:291602/149‑1 (MQ=255)                                               ggTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCtggt                                                                                       >  8:326626/1‑149 (MQ=255)                                                   tGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGc                                                                                    >  6:372374/1‑148 (MQ=255)                                                     gcaCAGGCCGAAGGCGCGCGGGGGCGATCGATGGTACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGTTACGATCTGGTGGCGAt                                                                                 <  4:263148/147‑1 (MQ=255)                                                     gAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGAt                                                                                 >  1:175390/1‑149 (MQ=255)                                                      aaCAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATc                                                                                <  4:275912/149‑1 (MQ=255)                                                         aGGACGAAGGCGCGCTGGGGCGAGCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTc                                                                             <  1:67759/149‑1 (MQ=255)                                                           gACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCa                                                                            >  7:310456/1‑148 (MQ=255)                                                           gACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCa                                                                            >  3:194513/1‑148 (MQ=255)                                                            aCGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCaaa                                                                          <  3:174639/149‑1 (MQ=255)                                                            aCGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCaa                                                                           >  2:357722/1‑148 (MQ=255)                                                                   cgcgcTGGGGCGATCGATGGAACGCTCTCAACAGTTCCTGCTGCTTTCGGCGCGTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCGTCAAAACGCTgg                                                                   <  4:219589/149‑1 (MQ=255)                                                                   cgcgcTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCTGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTgg                                                                   <  4:163494/149‑1 (MQ=255)                                                                   cgcgcTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCTGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTgg                                                                   <  4:157930/149‑1 (MQ=255)                                                                    gcgcTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTggg                                                                  <  7:326626/149‑1 (MQ=255)                                                                          gggCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAgg                                                            <  3:78385/149‑1 (MQ=255)                                                                              gatcgatGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGa                                                        <  5:155176/149‑1 (MQ=255)                                                                                  gatGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAgcg                                                     >  1:78319/1‑148 (MQ=255)                                                                                   atGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAgcgc                                                    >  3:334921/1‑148 (MQ=255)                                                                                    tGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCaa                                                  <  4:334921/149‑1 (MQ=255)                                                                                     ggAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAAc                                                 >  7:21136/1‑149 (MQ=255)                                                                                        aCGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGc                                              <  8:21136/149‑1 (MQ=255)                                                                                        aCGCGCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGc                                              <  4:344718/149‑1 (MQ=255)                                                                                         cGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCg                                             >  8:131082/1‑149 (MQ=255)                                                                                              gCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTgg                                                                   >  2:45252/1‑122 (MQ=255)                                                                                              gCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTgg                                                                   <  1:45252/122‑1 (MQ=255)                                                                                              gCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGc                                        >  7:285349/1‑149 (MQ=255)                                                                                                aaCAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCt                                       >  3:193001/1‑148 (MQ=255)                                                                                                  cAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAgcg                                                     <  6:296485/132‑1 (MQ=255)                                                                                                  cAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAgcg                                                     >  5:296485/1‑132 (MQ=255)                                                                                                  cAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAAt                                    >  5:213901/1‑149 (MQ=255)                                                                                                   aGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAAtc                                   <  2:78319/149‑1 (MQ=255)                                                                                                        ctgctgCTTTCGGGTGTTCTGACCTTGCTGCTGGGAGTGGCATCGGTGGCGGTAGAGATTAATCATTACTGTCGCAGTCGCTATGATGTGGTGGCGTTCCTCAAATCGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAAtcgtcg                               <  6:312651/148‑1 (MQ=255)                                                                                                        ctgctgCATTCGGCGCTACTGACCTTGATGCTGGCAGTGGCAGCGGTGGCGGTAGCTAAGAATCATTACTGTCGCTGTCGCTACGATCTGGTGGCGATCCTCAACACGCTGGGCGCAGGGCGAGCGCAACTGCGTAAGCTAAGCGTCgg                              >  8:159055/1‑149 (MQ=255)                                                                                                           ctgctTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTAAGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGt                             <  4:356753/147‑1 (MQ=255)                                                                                                            tgctTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAg                          >  6:149875/1‑149 (MQ=255)                                                                                                              ctTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAGTg                        <  5:222573/149‑1 (MQ=255)                                                                                                                  cGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAGTGGtt                     >  2:264003/1‑148 (MQ=255)                                                                                                                     cgcTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAGTGGTTGATg                 <  1:264003/149‑1 (MQ=255)                                                                                                                             aCCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAACCGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAGTGGTTGATGGTGCTGAc         >  1:231224/1‑149 (MQ=255)                                                                                                                               cTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAGTGGTTGATGGTGCTGACGc       >  7:371178/1‑149 (MQ=255)                                                                                                                                ttGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGCGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAGTGGTTGATGGTGCTGACGCt      >  6:16146/1‑149 (MQ=255)                                                                                                                                    tgctgGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAGTGGtt                     <  3:180292/130‑1 (MQ=255)                                                                                                                                    tgctgGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAGTGGtt                     >  4:180292/1‑130 (MQ=255)                                                                                                                                      ctgGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAGTGGTTGATGGTGCTGACGCTTTCa  <  8:285349/147‑1 (MQ=255)                                                                                                                                           |                                                                                                                                              TGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGGTCTGGAACAGGACGAAGGCGCGCTGGGGCGATCGATGGAACGCTCGCAACAGTTCCTGCTGCTTTCGGCGCTTCTGACCTTGCTGCTGGCAGTGGCAGCGGTGGCGGTAGCGATGAATCATTACTGTCGCAGTCGCTACGATCTGGTGGCGATCCTCAAAACGCTGGGGGCAGGGCGAGCGCAACTGCGTAAGCTAATCGTCGGTCAGTGGTTGATGGTGCTGACGCTTTCA  >  NC_000913/521042‑521321

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