Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,771,984	A→C	intergenic (‑240/‑297)	yibI ← / → mtlA	DUF3302 family inner membrane protein/mannitol‑specific PTS enzyme: IIA, IIB and IIC components

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,771,984	0	A	C	100.0%	267.9 / NA	68	intergenic (‑240/‑297)	yibI/mtlA	DUF3302 family inner membrane protein/mannitol‑specific PTS enzyme: IIA, IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (31/37);  total (31/37)

GATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTG  >  NC_000913/3771836‑3772129                                                                                                                                                     |                                                                                                                                                  gATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTccac                                                                                                                                                 >  2:14881/1‑151 (MQ=255) gATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTccac                                                                                                                                                 >  2:190915/1‑151 (MQ=255) gATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTccac                                                                                                                                                 <  2:692488/151‑1 (MQ=255)  aTTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTccaca                                                                                                                                                <  2:216697/151‑1 (MQ=255)  aTTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTccac                                                                                                                                                 >  2:700228/1‑150 (MQ=255)       cGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAAc                                                                                                                                             <  2:806709/149‑1 (MQ=255)        gCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCAt                                                                                                                                          <  1:365566/151‑1 (MQ=255)           aTAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATcaa                                                                                                                                       <  2:80426/151‑1 (MQ=255)           aTAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATca                                                                                                                                        >  1:138912/1‑150 (MQ=255)           aTAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATca                                                                                                                                        >  2:472022/1‑150 (MQ=255)             aaCACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATcaaa                                                                                                                                      >  2:211985/1‑150 (MQ=255)                 ccAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATcaaaacaa                                                                                                                                  >  2:326839/1‑150 (MQ=255)                  cAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATcaaaacaaa                                                                                                                                 >  2:602536/1‑150 (MQ=255)                     aGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAAtg                                                                                                                             >  1:770378/1‑151 (MQ=255)                          aatTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGaca                                                                                                                        <  2:346078/151‑1 (MQ=255)                             ttttAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACact                                                                                                                      >  2:117719/1‑150 (MQ=255)                             ttttAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACact                                                                                                                      >  1:271049/1‑150 (MQ=255)                                   agGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTcaca                                                                                                               <  2:355328/151‑1 (MQ=255)                                    ggTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACAt                                                                                                              <  2:8731/151‑1 (MQ=255)                                          gaGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTaaa                                                                                                         <  1:193869/150‑1 (MQ=255)                                           aGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATg                                                                                                       <  2:354723/151‑1 (MQ=255)                                               gCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCa                                                                                                    <  1:223722/150‑1 (MQ=255)                                                 ttAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCAttt                                                                                                 >  2:97869/1‑151 (MQ=255)                                                 ttAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCAttt                                                                                                 >  1:297605/1‑151 (MQ=255)                                                  tAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCAtttt                                                                                                >  2:194982/1‑151 (MQ=255)                                                  tAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCAttt                                                                                                 <  2:821563/150‑1 (MQ=255)                                                   aGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCAttttt                                                                                               <  2:107985/151‑1 (MQ=255)                                                    gTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTa                                                                                              <  2:811145/151‑1 (MQ=255)                                                    gTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTa                                                                                              <  2:465170/151‑1 (MQ=255)                                                       aCATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGc                                                                                            <  2:155345/150‑1 (MQ=255)                                                       aCATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCg                                                                                           <  2:514559/151‑1 (MQ=255)                                                         aTAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGaa                                                                                         <  2:577474/151‑1 (MQ=255)                                                         aTAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGaa                                                                                         <  2:552142/151‑1 (MQ=255)                                                             cGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATc                                                                                     <  1:863/151‑1 (MQ=255)                                                               aTTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACGAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATcgc                                                                                   <  2:836558/151‑1 (MQ=255)                                                                ttGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATcgcc                                                                                  >  1:376105/1‑151 (MQ=255)                                                                 tGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATcgcc                                                                                  <  1:645853/150‑1 (MQ=255)                                                                 tGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATcgcc                                                                                  <  1:460405/150‑1 (MQ=255)                                                                    aTGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCt                                                                              <  2:830772/151‑1 (MQ=255)                                                                     tGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCtt                                                                             >  1:883011/1‑151 (MQ=255)                                                                     tGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCtt                                                                             >  2:542993/1‑151 (MQ=255)                                                                      gACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCttg                                                                            <  1:907835/151‑1 (MQ=255)                                                                       aCGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCttgt                                                                           >  1:752776/1‑151 (MQ=255)                                                                              cATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCtttt                                                                    >  1:251750/1‑151 (MQ=255)                                                                                   cATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTac                                                                >  1:574578/1‑150 (MQ=255)                                                                                           aGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGtt                                                       >  1:446179/1‑151 (MQ=255)                                                                                            gATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGttt                                                      <  1:69516/151‑1 (MQ=255)                                                                                             aTCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGtttt                                                     >  2:150517/1‑151 (MQ=255)                                                                                               cacaTCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTtgt                                                   >  1:87464/1‑151 (MQ=255)                                                                                               cacaTCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTtgt                                                   >  2:334717/1‑151 (MQ=255)                                                                                                    cAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGtgatg                                               >  2:636396/1‑150 (MQ=255)                                                                                                     aGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGtgatga                                              >  1:877967/1‑150 (MQ=255)                                                                                                       gTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAacg                                           <  2:37578/151‑1 (MQ=255)                                                                                                         gAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAacgtc                                         >  2:573980/1‑151 (MQ=255)                                                                                                          aaCGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAacgtc                                         >  2:866436/1‑150 (MQ=255)                                                                                                           aCGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAacgtcac                                       <  2:406369/151‑1 (MQ=255)                                                                                                              ccGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAacgtcacgtc                                    <  1:410460/151‑1 (MQ=255)                                                                                                              ccGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAacgtcacgtc                                    <  1:789240/151‑1 (MQ=255)                                                                                                                gTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCaa                                  <  1:852827/151‑1 (MQ=255)                                                                                                                  aaGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAAtt                                >  1:535953/1‑151 (MQ=255)                                                                                                                  aaGAAAATATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAAt                                 <  1:299712/150‑1 (MQ=255)                                                                                                                     aaaaTATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTAcc                             >  2:553235/1‑151 (MQ=255)                                                                                                                     aaaaTATCTTGTGATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTAc                              <  2:638261/150‑1 (MQ=255)                                                                                                                              tgtgATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACctctct                        <  2:689246/147‑1 (MQ=255)                                                                                                                              tgtgATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACctctct                        <  2:196816/147‑1 (MQ=255)                                                                                                                              tgtgATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTAcc                     <  1:194985/150‑1 (MQ=255)                                                                                                                              tgtgATTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTAcc                     <  2:179800/150‑1 (MQ=255)                                                                                                                                  aTTCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCta                <  2:633346/151‑1 (MQ=255)                                                                                                                                    tCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCtat               >  1:37388/1‑150 (MQ=255)                                                                                                                                    tCAGATCACAAAGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCtat               >  2:475117/1‑150 (MQ=255)                                                                                                                                              aaGATTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGtgat    <  1:616931/151‑1 (MQ=255)                                                                                                                                                 aTTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGtgattg  >  1:799766/1‑150 (MQ=255)                                                                                                                                                 aTTCCACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGtgattg  >  2:334921/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTG  >  NC_000913/3771836‑3772129

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