Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	237,900	A→T	6.0%	D189V (GAC→GTC)	yafT →	lipoprotein YafT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	237,900	0	A	T	6.0%	202.9 / 3.4	67	D189V (GAC→GTC)	yafT	lipoprotein YafT
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base T (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

AACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGCCAGTACGACTGATATGACAGAGC  >  NC_000913/237755‑238049                                                                                                                                                  |                                                                                                                                                      aaCTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCt                                                                                                                                                   <  1:1020685/150‑1 (MQ=255)     ggATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTAtct                                                                                                                                               >  1:7184/1‑150 (MQ=255)        tGCAGGCACTGCGTTTTGTGGCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTccc                                                                                                                                            >  1:1310088/1‑150 (MQ=255)        tGCAGGCACTGCGTTTTGTGGCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTccc                                                                                                                                            <  2:458312/150‑1 (MQ=255)           aGGCACTGCGTTTTGTGGCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTgc                                                                                                                                         <  1:23996/150‑1 (MQ=255)              cACTGCGTTTTGTGGCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCt                                                                                                                                      >  1:448908/1‑150 (MQ=255)              cACTGCGTTTTGTGGCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCt                                                                                                                                      <  2:170510/150‑1 (MQ=255)                 tGCGTTTTGTGGCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTaac                                                                                                                                   >  1:1073284/1‑148 (MQ=255)                    gTTTTGTGGCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTgg                                                                                                                                <  2:29673/150‑1 (MQ=255)                           ggCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGttt                                                                                                                         >  1:678659/1‑150 (MQ=255)                             cTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTcc                                                                                                                       >  2:776790/1‑149 (MQ=255)                             cTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCa                                                                                                                       >  1:176919/1‑150 (MQ=255)                              tGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCAc                                                                                                                      >  2:423842/1‑150 (MQ=255)                                    aGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCtggt                                                                                                                <  1:319524/150‑1 (MQ=255)                                       aCATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGa                                                                                                             <  1:211809/150‑1 (MQ=255)                                        cATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGAt                                                                                                            <  1:423842/150‑1 (MQ=255)                                          tCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGAtgt                                                                                                          <  1:722601/150‑1 (MQ=255)                                           cAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGAtgtg                                                                                                         >  1:70984/1‑150 (MQ=255)                                             gAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGc                                                                                                       <  1:919053/150‑1 (MQ=255)                                             gAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGc                                                                                                       >  1:9927/1‑150 (MQ=255)                                                  cGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAg                                                                                                  >  1:940599/1‑150 (MQ=255)                                                   gATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAgg                                                                                                 >  2:940999/1‑150 (MQ=255)                                                    attattGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGt                                                                                                <  2:510329/150‑1 (MQ=255)                                                      tattGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGa                                                                                              >  1:893134/1‑150 (MQ=255)                                                          gTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTgg                                                                                          >  2:813119/1‑150 (MQ=255)                                                          gTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTgg                                                                                          >  1:249494/1‑150 (MQ=255)                                                           tttACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTggg                                                                                         <  2:355939/150‑1 (MQ=255)                                                            ttACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGc                                                                                        <  2:678659/150‑1 (MQ=255)                                                             tACCAGGATATGTTGCAGACAGTAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCt                                                                                       <  1:57681/150‑1 (MQ=255)                                                                 aGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGTCGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCt                                                                                   >  1:935184/1‑150 (MQ=255)                                                                 aGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGTCGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCt                                                                                   >  1:939929/1‑150 (MQ=255)                                                                        gTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGcc                                                                            <  1:823051/150‑1 (MQ=255)                                                                        gTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGcc                                                                            >  1:13902/1‑150 (MQ=255)                                                                         ttGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCg                                                                           <  1:85975/150‑1 (MQ=255)                                                                          tGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACg                                                                                                                            >  2:21629/1‑100 (MQ=255)                                                                          tGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACg                                                                                                                            <  1:21629/100‑1 (MQ=255)                                                                                cAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAAtt                                                                    >  1:1202466/1‑150 (MQ=255)                                                                                 aGGAAAATATGACTCTGGGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGGACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTa                                                                   <  2:925813/150‑1 (MQ=255)                                                                                   gAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATg                                                                 >  2:1087413/1‑150 (MQ=255)                                                                                   gAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTAGACGGTGAATTATg                                                                 >  1:1278645/1‑150 (MQ=255)                                                                                     aaaTATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGaa                                                               <  2:940599/150‑1 (MQ=255)                                                                                        tatGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAatat                                                            >  1:946978/1‑150 (MQ=255)                                                                                                  cgcTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTc                                                  >  1:446192/1‑150 (MQ=255)                                                                                                  cgcTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTc                                                  <  1:974813/150‑1 (MQ=255)                                                                                                         aTCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCAtt                                           >  2:1036180/1‑150 (MQ=255)                                                                                                           cAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGTCGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCTGTGAATTATGAATATAAAGTGCTTCCACCATTGc                                         <  2:939929/150‑1 (MQ=255)                                                                                                           cAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGTCGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCTGTGAATTATGAATATAAAGTGCTTCCACCATTGc                                         <  2:935184/150‑1 (MQ=255)                                                                                                           cAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGc                                         <  2:882862/150‑1 (MQ=255)                                                                                                           cAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGc                                         >  2:234541/1‑150 (MQ=255)                                                                                                             aCAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGccc                                       >  2:1131035/1‑150 (MQ=255)                                                                                                               aGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGa                                     <  1:131256/150‑1 (MQ=255)                                                                                                                gTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGAc                                    >  2:1375272/1‑150 (MQ=255)                                                                                                                    ggTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGa                                >  1:1067286/1‑150 (MQ=255)                                                                                                                     gTCCGACTATAAATATGACAAACTCACTGACGCTATCTCACTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGTCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGaa                               <  1:1379864/150‑1 (MQ=255)                                                                                                                     gTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGaa                               >  2:341624/1‑150 (MQ=255)                                                                                                                      tCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAAt                              <  1:1042104/150‑1 (MQ=255)                                                                                                                       ccGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATg                             >  2:906208/1‑150 (MQ=255)                                                                                                                       ccGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAAAAGGTGAGTGGGCTACCTGGTCGCCTGTGAATTATTAATATAAAGTTCTTCCACCATTGCCCGACAAGAATg                             <  2:603497/150‑1 (MQ=255)                                                                                                                        cGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGa                            >  2:1134579/1‑150 (MQ=255)                                                                                                                            tataAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGcc                        <  2:70984/150‑1 (MQ=255)                                                                                                                                 aaaTGACAAACTCACTGATGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGCCAGTAc                   >  2:512834/1‑150 (MQ=255)                                                                                                                                 aaaTGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGCCAGTAc                   <  1:1087413/150‑1 (MQ=255)                                                                                                                                     gACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGCCAGTACGACt               >  1:581088/1‑150 (MQ=255)                                                                                                                                          aCTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGCCAGTACGACTGatat          <  1:906208/150‑1 (MQ=255)                                                                                                                                           cTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGCCAGTACGACTGATATg         <  2:7184/150‑1 (MQ=255)                                                                                                                                           cTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGCCAGTACGACTGATATg         <  1:940999/150‑1 (MQ=255)                                                                                                                                                tGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGCCAGTACGACTGATATGACaga    <  2:1067286/150‑1 (MQ=255)                                                                                                                                                  aCGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGCCAGTACGACTGATATGACAGagc  >  1:565839/1‑150 (MQ=255)                                                                                                                                                  aCGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGGGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGCCAGTACGACTGATATGACAGagc  <  2:983982/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      AACTGGATGCAGGCACTGCGTTTTGTGGCTGCAAAAGGACATCAGAAAGCGATTATTGTTTACCAGGATATGTTGCAGACAGGAAAATATGACTCTGCGCTGAAATCAACAGTATGGTCCGACTATAAAAATGACAAACTCACTGACGCTATCTCCCTGCGCTACCTGGTACGTTTCACGCTGGTGGATGTGGCAACAGGTGAGTGGGCTACCTGGTCGCCGGTGAATTATGAATATAAAGTGCTTCCACCATTGCCCGACAAGAATGAAGCCAGTACGACTGATATGACAGAGC  >  NC_000913/237755‑238049

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