Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,673,783	G→A	10.2%	intergenic (+15/+33)	yphA → / ← yphB	DoxX family inner membrane protein/mutarotase superfamily protein, YphB family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,673,783	0	G	A	10.2%	175.4 / 8.8	69	intergenic (+15/+33)	yphA/yphB	DoxX family inner membrane protein/mutarotase superfamily protein, YphB family
Reads supporting (aligned to +/- strand):  ref base G (34/27);  new base A (4/3);  total (39/30)
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.58e-01

TACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTT  >  NC_000913/2673640‑2673931                                                                                                                                                |                                                                                                                                                     tACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATgcggcg                                                                                                                                                 >  1:206240/1‑149 (MQ=255)  aCTGGGATATGACCGGCGATGAGGTTGGGCCAAATATGATTAATATCTGGAAGAATGTCCGTATCGCTGGCGCGTTCTTGCTATTGTCAATTACCGGGCCGGGTGCAATTTATCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGt                                                                                                                                                <  1:860058/149‑1 (MQ=255)    tGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGa                                                                                                                                              >  1:214573/1‑149 (MQ=255)        atatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGc                                                                                                                                          >  1:336425/1‑149 (MQ=255)         tatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGcc                                                                                                                                         <  1:584735/149‑1 (MQ=255)         tatGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGcc                                                                                                                                         <  1:731950/149‑1 (MQ=255)           tGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                       <  1:625612/149‑1 (MQ=255)            gACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                      >  1:1006669/1‑149 (MQ=255)                     tGCGGTTGGGCCAAATCTGATTAATTTCTGGAAGAATGTCTGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGTGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATAAGGCGTTCACGCCGCATCCGGCAtc                                                                                                                             <  2:414539/149‑2 (MQ=38)                        ggTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTTCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaa                                                                                                                          <  1:912094/149‑1 (MQ=255)                         gTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTTCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTtaat                                                                                                                         >  2:770350/1‑149 (MQ=255)                              gCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAAtt                                                                                                                    >  1:352567/1‑149 (MQ=255)                                cAACTATGATTAATTTCTGGAAGAATGTCAGTATCGCCGGCGCGTTCGTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCTTTAGGATAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTAc                                                                                                                  <  1:761363/149‑1 (MQ=255)                                 aaaTATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGTATTTAATAAATTACa                                                                                                                 >  1:891409/1‑149 (MQ=255)                                   atatGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGc                                                                                                               <  1:589393/149‑1 (MQ=255)                                       gATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCAct                                                                                                           >  2:287687/1‑149 (MQ=255)                                        aTTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActc                                                                                                          >  2:851006/1‑149 (MQ=255)                                         ttAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCActct                                                                                                         >  2:13810/1‑149 (MQ=255)                                            aTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                                                                      >  2:460306/1‑149 (MQ=255)                                            aTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTAcc                                                                                                      >  1:1079493/1‑149 (MQ=255)                                            aTTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTAGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTAAAGCCACTCTAcc                                                                                                      >  1:1021694/1‑149 (MQ=255)                                                   gaagaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACg                                                                                               >  1:1106597/1‑149 (MQ=255)                                                     agaaTGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGAc                                                                                             <  2:336425/149‑1 (MQ=255)                                                         tGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCCGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATATAATAAATTACAGCCACTCTACCCGCAACGACAtca                                                                                         >  2:1129155/1‑148 (MQ=255)                                                           tCAGGATCGCTGGCGCGTTCTTGCTATTTGCTATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATctct                                                                                       <  2:28660/149‑1 (MQ=255)                                                             aGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGa                                                                                     >  1:115007/1‑149 (MQ=255)                                                              gTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGaa                                                                                    <  2:1067419/149‑1 (MQ=255)                                                                   gCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATTCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAggtgt                                                                               >  2:544536/1‑149 (MQ=255)                                                                        cgcgTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTCCCCGCACCGCCATCTCTGAAGTTGTTGAtt                                                                          >  2:773634/1‑149 (MQ=255)                                                                         gcgTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTc                                                                         >  2:857683/1‑149 (MQ=255)                                                                         gcgTGCTTGCTATTGGCAATTACCGGGCCGGGGGCAAGTTCTCTCGATCGGCGTTAGGAAAGTCGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTc                                                                         <  1:541369/149‑1 (MQ=255)                                                                               ttGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACTTCGCTGAAGTTGTTGATTCCCCTgg                                                                   >  1:329210/1‑149 (MQ=255)                                                                               ttGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGGGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTgg                                                                   <  1:113790/149‑1 (MQ=255)                                                                                   tATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGca                                                               >  2:390401/1‑148 (MQ=255)                                                                                    aTTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTa                                                              <  1:785083/149‑1 (MQ=255)                                                                                    aTTGGCAATTAACGGGCCGGGTGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTa                                                              <  2:9939/149‑1 (MQ=255)                                                                                     ttGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAATATGACGGATGCCGCGTGAACGCCTTATCCGGCATATTATAAATTACAGACACTCTACACGAAACGCAATCTCTGAAGTTGTTTATTCCCATGGTGCTAg                                                             >  1:105727/1‑149 (MQ=255)                                                                                     ttGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCCACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg                                                             >  1:749880/1‑149 (MQ=255)                                                                                     ttGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAg                                                             >  1:807983/1‑149 (MQ=255)                                                                                        gCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGCACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAgcgc                                                          >  2:76685/1‑149 (MQ=255)                                                                                          aaTTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACTCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAAGGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTTGCGCaa                                                        >  2:525856/1‑149 (MQ=255)                                                                                          aaTTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCaa                                                        >  2:560424/1‑149 (MQ=255)                                                                                           aTTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAAt                                                       <  1:671893/149‑1 (MQ=255)                                                                                            ttACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATg                                                      <  2:1079493/149‑1 (MQ=255)                                                                                             tACCGGGCAGGGGGCAATTGCTCTCGATCGGCGGTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAACTTACAGCCACTCTACCCTCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATga                                                     <  2:912666/149‑1 (MQ=255)                                                                                                cgggccgggGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGAt                                                  >  1:407321/1‑149 (MQ=255)                                                                                                      gggggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCAcccc                                            >  2:431727/1‑149 (MQ=255)                                                                                                       ggggCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCAccccc                                           >  2:321765/1‑149 (MQ=255)                                                                                                          gCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTc                                        <  1:460306/149‑1 (MQ=255)                                                                                                            aaTTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCGTATCCGGCATTTAATAAATTAAAGCAACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTACCCTGGTGCTAGCGCAATGAGATCAcccaattagg                                      >  2:547090/1‑142 (MQ=255)                                                                                                                  tctcGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgat                                <  1:857683/149‑1 (MQ=255)                                                                                                                     cGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGTATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgat                             <  2:891409/149‑1 (MQ=255)                                                                                                                      gATCGGCGTTATGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGAAACGACATCTCTGACGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgatg                            <  1:770350/149‑1 (MQ=255)                                                                                                                      gATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgatg                            >  2:1010602/1‑149 (MQ=255)                                                                                                                       aTCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCGCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCgatgatga                           <  1:797052/149‑1 (MQ=255)                                                                                                                          ggCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGatca                        <  2:105285/149‑1 (MQ=255)                                                                                                                            cGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCGTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGCCATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGatcatc                      >  2:149845/1‑149 (MQ=255)                                                                                                                              ttAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCg                    <  1:653922/149‑1 (MQ=255)                                                                                                                                aGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgc                  <  2:214573/149‑1 (MQ=255)                                                                                                                                 ggAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgcg                 >  2:233469/1‑149 (MQ=255)                                                                                                                                  attaaaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCGCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgcgc                >  1:1083705/6‑149 (MQ=38)                                                                                                                                  attaaaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGgcgc                >  1:670839/6‑149 (MQ=38)                                                                                                                                    aaGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCAt              <  1:1172035/149‑1 (MQ=255)                                                                                                                                     aaaTGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATAATCCGGCGCATg             >  1:1092687/3‑149 (MQ=38)                                                                                                                                      gATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATgg            <  1:996022/149‑1 (MQ=255)                                                                                                                                      gATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATgg            <  2:456760/149‑1 (MQ=255)                                                                                                                                          ccGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCa        >  2:153167/1‑149 (MQ=255)                                                                                                                                           cGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCACCGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTcct       >  2:723659/1‑147 (MQ=255)                                                                                                                                           cGGATAAGGCGTTCACGCCGCATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCag       >  1:388600/1‑148 (MQ=255)                                                                                                                                               tGCGGCGGGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCCACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCGAGCGCAATGAGCTCACCCCCTTCCGGGCGATGAGGATCATCCGGCGCATGGCTCATCGGt   <  2:1153362/149‑1 (MQ=255)                                                                                                                                                gcggcgTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGtt  >  1:176692/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     TACTGGGATATGACCGGCGATGCGGTTGGGCCAAATATGATTAATTTCTGGAAGAATGTCAGTATCGCTGGCGCGTTCTTGCTATTGGCAATTACCGGGCCGGGGGCAATTTCTCTCGATCGGCGTTAGGAAAGATGCCGGATGCGGCGTGAACGCCTTATCCGGCATTTAATAAATTACAGCCACTCTACCCGCAACGACATCTCTGAAGTTGTTGATTCCCCTGGTGCTAGCGCAATGAGATCACCCCCTTCCGGGCGATGATGATCATCCGGCGCATGGCTCATCGGTT  >  NC_000913/2673640‑2673931

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