Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,261,155	A→T	D186V (GAT→GTT)	yjbM →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,261,155	0	A	T	96.6%	297.8 / 0.7	88	D186V (GAT→GTT)	yjbM	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (46/39);  total (48/40)
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.23e-01

ACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTT  >  NC_000913/4261059‑4261252                                                                                                 |                                                                                                  aCCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTcc                                                                                                <  2:421806/100‑1 (MQ=255) aCCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTcc                                                                                                >  1:771930/1‑100 (MQ=255)   cTGGGAGAGATGAAGCGTCTGGCTGATGTGGGCGAAATGTATTTTGATTGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGGACGGTCGTTcc                                                                                                <  2:1076630/98‑1 (MQ=255)   cTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTg                                                                                              >  1:378252/1‑100 (MQ=255)    tGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGa                                                                                             <  2:1290994/100‑1 (MQ=255)    tGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGa                                                                                             <  1:1866982/100‑1 (MQ=255)    tGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGa                                                                                             <  1:1805694/100‑1 (MQ=255)      ggAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGata                                                                                           >  1:1280874/1‑100 (MQ=255)      ggAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGata                                                                                           >  2:1902592/1‑100 (MQ=255)       gagagaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGatat                                                                                          <  2:1293200/100‑1 (MQ=255)           gaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTa                                                                                      >  2:1383393/1‑100 (MQ=255)           gaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTa                                                                                      <  2:78105/100‑1 (MQ=255)             tGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTGACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAg                                                                                    >  2:1432772/1‑100 (MQ=255)             tGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAg                                                                                    >  1:57504/1‑100 (MQ=255)                aGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCaa                                                                                 >  1:1880343/1‑100 (MQ=255)                 gCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCaaa                                                                                >  1:321919/1‑100 (MQ=255)                  cGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCaaaa                                                                               >  2:320720/1‑100 (MQ=255)                   gTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAAt                                                                              >  1:29072/1‑100 (MQ=255)                      tGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCg                                                                           <  2:965605/100‑1 (MQ=255)                       ggCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCga                                                                          <  1:1432781/100‑1 (MQ=255)                        gCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCgag                                                                         <  2:61619/100‑1 (MQ=255)                        gCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCgag                                                                         >  2:908984/1‑100 (MQ=255)                           gATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCATTCCTGATATTTTAAGCAAAATGCGAGAtg                                                                      <  2:1088896/100‑1 (MQ=255)                            aTGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGAtgc                                                                     <  1:897361/100‑1 (MQ=255)                             tGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGAtgct                                                                    >  2:1991461/1‑100 (MQ=255)                               ggggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGAtgctgc                                                                  <  1:1350198/100‑1 (MQ=255)                               ggggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGAtgctgc                                                                  <  1:1465426/100‑1 (MQ=255)                                gggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGca                                                                 <  1:524131/100‑1 (MQ=255)                                 ggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGcac                                                                >  2:625060/1‑100 (MQ=255)                                  gagaAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGcacc                                                               <  1:837167/100‑1 (MQ=255)                                  gagaAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGcacc                                                               <  2:2008177/100‑1 (MQ=255)                                   agaAATGAATTTTGATGGCTATCTGGATAATATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCCGATATTTTAAGCAAAATGCGAGATGCTGcacca                                                              <  2:1809803/100‑1 (MQ=255)                                      aaTGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaa                                                           >  1:395746/1‑100 (MQ=255)                                        tGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaatt                                                         >  1:1883309/1‑100 (MQ=255)                                         gAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaatta                                                        >  1:187495/1‑100 (MQ=255)                                         gAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaatta                                                        <  1:1161139/100‑1 (MQ=255)                                         gAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaatta                                                        >  2:775315/1‑100 (MQ=255)                                          aaTTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaa                                                       >  2:923800/1‑100 (MQ=255)                                           aTTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaat                                                      >  2:23180/1‑100 (MQ=255)                                            ttttGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaatt                                                     <  2:771927/100‑1 (MQ=255)                                               tgtttgCTATCTGGATTCGCTGTCACATCTTGGGAATGAAGGATCGCTCTTTCCTGAGATTTTAAGCAAGATGCGAGATGCTGCACCACaat                                                          <  2:1937284/87‑1 (MQ=255)                                                gATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCt                                                 <  1:320721/100‑1 (MQ=255)                                                  tGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTc                                               <  2:1402748/100‑1 (MQ=255)                                                      tATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGa                                           <  2:343810/100‑1 (MQ=255)                                                       aTCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGaa                                          >  2:1222637/1‑100 (MQ=255)                                                        tCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAg                                         <  1:1797911/100‑1 (MQ=255)                                                         cTGGATAAGATGTCACACCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGt                                        >  2:1506346/1‑100 (MQ=255)                                                          tGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTc                                       >  2:446965/1‑100 (MQ=255)                                                            gATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGt                                     >  1:1983360/1‑100 (MQ=255)                                                             aTAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGtt                                    <  2:1880325/100‑1 (MQ=255)                                                             aTAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGtt                                    <  2:80810/100‑1 (MQ=255)                                                              tAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGttt                                   <  1:1999147/100‑1 (MQ=255)                                                               aaGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTg                                  >  2:970631/1‑100 (MQ=255)                                                                aGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGa                                 >  1:819060/1‑100 (MQ=255)                                                                aGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGa                                 >  2:2035239/1‑100 (MQ=255)                                                                 gATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAc                                <  2:29072/100‑1 (MQ=255)                                                                   tGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAccc                              >  2:1597612/1‑100 (MQ=255)                                                                   tGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAccc                              >  2:710008/1‑100 (MQ=255)                                                                     tCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCaa                            >  1:312915/1‑100 (MQ=255)                                                                      cacaTCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAc                           >  1:1489009/1‑100 (MQ=255)                                                                      cacaTCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAc                           >  1:1281016/1‑100 (MQ=255)                                                                       acaTCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAcc                          >  1:1087112/1‑100 (MQ=255)                                                                       acaTCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAcc                          >  1:1087127/1‑100 (MQ=255)                                                                           cTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCaa                      >  2:1214721/1‑100 (MQ=255)                                                                             ggTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGATGTCGTTTGACCCAACCTCAAAg                    <  2:1185575/100‑1 (MQ=255)                                                                              gTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCAGAATTAATTAGCTTCGCGAAGTCGTTTGATCCAACCTCaaatg                   <  1:1383402/100‑3 (MQ=255)                                                                              aaaatatcaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAgg                   >  2:1055551/9‑100 (MQ=255)                                                                               tgaatgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGa                  >  2:667578/1‑100 (MQ=255)                                                                                gaatgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGaa                 >  2:903246/1‑100 (MQ=255)                                                                                 aatgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAag                <  2:1677193/100‑1 (MQ=255)                                                                                   tgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAagag              >  2:431150/1‑100 (MQ=255)                                                                                   tgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAagag              <  1:1902610/100‑1 (MQ=255)                                                                                    gaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAgaga             >  1:572951/1‑100 (MQ=255)                                                                                     aaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAt            <  1:1991479/100‑1 (MQ=255)                                                                                     aaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAt            >  1:1422522/1‑100 (MQ=255)                                                                                      aGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAtt           <  1:1796971/100‑1 (MQ=255)                                                                                        gAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTaa         >  1:70816/1‑100 (MQ=255)                                                                                        gAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTaa         >  2:143787/1‑100 (MQ=255)                                                                                         aaCGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTaaa        <  1:908988/100‑1 (MQ=255)                                                                                          aCGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTaaaa       <  1:775318/100‑1 (MQ=255)                                                                                          aCGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTaaaa       >  2:556671/1‑100 (MQ=255)                                                                                          aCGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTaaaa       <  1:904839/100‑1 (MQ=255)                                                                                           cGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAAt      <  2:1280868/100‑1 (MQ=255)                                                                                            gCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATa     >  1:1685807/1‑100 (MQ=255)                                                                                             cTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATAc    >  1:1094428/1‑100 (MQ=255)                                                                                             cTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATAc    >  1:686287/1‑100 (MQ=255)                                                                                              tCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACt   <  2:395745/100‑1 (MQ=255)                                                                                               cGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACtt  >  1:1930823/1‑100 (MQ=255)                                                                                                 |                                                                                                  ACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTT  >  NC_000913/4261059‑4261252

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