Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,348,462	T→G	L4* (TTA→TGA)	yrbL →	Mg(2+)‑starvation‑stimulated protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,348,462	0	T	G	100.0%	243.3 / NA	63	L4* (TTA→TGA)	yrbL	Mg(2+)‑starvation‑stimulated protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (35/28);  total (35/28)

AATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTTATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATAAAGAGATCCGCCGTGAGTTAAAGTACTACGCGC  >  NC_000913/3348314‑3348605                                                                                                                                                     |                                                                                                                                                aatcaaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGAt                                                                                                                                               >  2:1134555/1‑151 (MQ=255) aatcaaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGAt                                                                                                                                               >  1:639255/1‑151 (MQ=255)    caaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTg                                                                                                                                            >  1:818664/1‑151 (MQ=255)     aaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTg                                                                                                                                            >  2:1074988/1‑150 (MQ=255)      aCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGaa                                                                                                                                          <  1:1046679/151‑1 (MQ=255)      aCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGaa                                                                                                                                          <  2:145188/151‑1 (MQ=255)       cGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAAc                                                                                                                                         >  2:508848/1‑151 (MQ=255)       cGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAAc                                                                                                                                         >  2:366856/1‑151 (MQ=255)         tCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGCCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACa                                                                                                                                        >  2:644226/1‑150 (MQ=255)                 aCATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTcccc                                                                                                                                <  1:937792/150‑1 (MQ=255)                  cATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTg                                                                                                                              >  2:983554/1‑151 (MQ=255)                   aTTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTgg                                                                                                                             <  1:116258/151‑1 (MQ=255)                       tgtgCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCAc                                                                                                                         >  2:414675/1‑151 (MQ=255)                          gcgcTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCgg                                                                                                                      <  1:568068/151‑1 (MQ=255)                               tttCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGa                                                                                                                  <  2:63449/150‑1 (MQ=255)                                 tCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACAt                                                                                                               >  2:73862/1‑151 (MQ=255)                                  cACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATc                                                                                                              <  1:74894/151‑1 (MQ=255)                                   aCTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCg                                                                                                             >  1:565828/1‑151 (MQ=255)                                        cAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAAt                                                                                                        >  2:218322/1‑151 (MQ=255)                                          gaagaaCCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGt                                                                                                      >  2:176707/1‑151 (MQ=255)                                              aaCCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATg                                                                                                  <  1:113291/151‑1 (MQ=255)                                                 cTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATgcgc                                                                                               >  1:1170653/1‑151 (MQ=255)                                                 cTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATgcgc                                                                                               >  2:326722/1‑151 (MQ=255)                                                    aaGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCAt                                                                                             <  2:490486/150‑1 (MQ=255)                                                    aaGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATc                                                                                            <  2:580536/151‑1 (MQ=255)                                                    aaGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATc                                                                                            >  1:53654/1‑151 (MQ=255)                                                      gAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCg                                                                                          <  1:905114/151‑1 (MQ=255)                                                      gAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCg                                                                                          <  2:685069/151‑1 (MQ=255)                                                            cTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGa                                                                                        >  2:1043488/1‑147 (MQ=255)                                                              tAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGAt                                                                                    <  2:1097908/149‑1 (MQ=255)                                                              tAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGc                                                                                  >  1:720248/1‑151 (MQ=255)                                                                 gagGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGccc                                                                                >  2:643284/1‑150 (MQ=255)                                                                 gagGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGcc                                                                                 >  2:1146752/1‑149 (MQ=255)                                                                 gagGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCa                                                                               <  1:610798/151‑1 (MQ=255)                                                                  agGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCa                                                                               >  2:736694/1‑150 (MQ=255)                                                                  agGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCa                                                                               >  2:462476/1‑150 (MQ=255)                                                                   ggCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAAc                                                                             <  1:708132/151‑1 (MQ=255)                                                                    gCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACg                                                                            <  2:83031/151‑1 (MQ=255)                                                                           ttAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATc                                                                     >  1:897023/1‑151 (MQ=255)                                                                            tAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATc                                                                     <  1:604470/150‑1 (MQ=255)                                                                                   tGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGt                                                             >  2:94537/1‑151 (MQ=255)                                                                                     tttAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCt                                                           >  2:327094/1‑151 (MQ=255)                                                                                         aGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCa                                                       <  2:208740/151‑1 (MQ=255)                                                                                          gTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCAt                                                      >  2:1023214/1‑151 (MQ=255)                                                                                            aaaTCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATc                                                     <  1:754450/149‑1 (MQ=255)                                                                                             aaTCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGt                                                   >  2:831749/1‑151 (MQ=255)                                                                                                    cATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGAtg                                            >  2:502691/1‑151 (MQ=255)                                                                                                    cATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGAtg                                            >  2:853247/1‑151 (MQ=255)                                                                                                       aCTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGAtggtg                                         >  1:874982/1‑151 (MQ=255)                                                                                                         tGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGCCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGc                                       >  1:400427/1‑151 (MQ=255)                                                                                                          ggAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGc                                       <  1:1149663/150‑1 (MQ=255)                                                                                                           gagaTCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGa                                     <  1:998889/151‑1 (MQ=255)                                                                                                             gaTCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATa                                   <  2:545158/151‑1 (MQ=255)                                                                                                                cGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATAAag                                <  2:874570/151‑1 (MQ=255)                                                                                                                 gTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATaa                                  >  2:214574/1‑148 (MQ=255)                                                                                                                 gTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATAAaga                               <  2:565802/151‑1 (MQ=255)                                                                                                                  tCAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATAAag                                <  2:258879/149‑1 (MQ=255)                                                                                                                   cAGAAAATATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATAAag                                >  2:577344/1‑148 (MQ=255)                                                                                                                      aaaaTATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATAAAGAGATcc                          <  2:302153/151‑1 (MQ=255)                                                                                                                       aaaTATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATAAAGAGATccg                         <  1:524775/151‑1 (MQ=255)                                                                                                                       aaaTATTTCCAGGAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATAAAGAGATccg                         <  2:1106319/151‑1 (MQ=255)                                                                                                                                  ggAGATGGCATGATTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATAAAGAGATCCGCCGTGAGTTa               >  1:622241/1‑150 (MQ=255)                                                                                                                                              aTTCGCTGATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATAAAGAGATCCGCCGTGAGTTAAAGTACTACgcgc  >  1:205203/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                AATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTTTAAGTTAATCGACCATACTGGAGATCGTCAGAAAATATTTCCAGGAGATGGCATGATTCGCTTATCTGAACAAAGTCCCCTGGGCACCGGAAGACATCGGAAATGTTATGCGCATCCGGAAGATGCCCAACGCTGTATCAAGATTGTCTACCATCGTGGCGATGGTGGCGATAAAGAGATCCGCCGTGAGTTAAAGTACTACGCGC  >  NC_000913/3348314‑3348605

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