Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,541,962	T→G	100%	I519S (ATC→AGC)	feoB →	ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,541,962	0	T	G	100.0%	276.6 / NA	79	I519S (ATC→AGC)	feoB	ferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (45/34);  total (45/34)

TGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCA  >  NC_000913/3541866‑3542057                                                                                                 |                                                                                                tGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagcat                                                                                              >  2:2456830/1‑100 (MQ=255) tGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagca                                                                                               <  2:246667/99‑1 (MQ=255)   gAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagcatca                                                                                            <  1:475509/100‑1 (MQ=255)   gAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagcatca                                                                                            <  2:1626266/100‑1 (MQ=255)    aGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagcatcat                                                                                           >  1:2550415/1‑100 (MQ=255)      ttGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCg                                                                                         <  1:821625/100‑1 (MQ=255)      ttGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCg                                                                                         <  1:765978/100‑1 (MQ=255)      ttGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCg                                                                                         <  1:930771/100‑1 (MQ=255)       tGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGt                                                                                        >  2:2494808/1‑100 (MQ=255)        gCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTc                                                                                       >  2:402274/1‑100 (MQ=255)        gCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTACTGGTAAAGTGAGCATCATCGTc                                                                                       >  1:854337/1‑100 (MQ=255)         ccGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCa                                                                                      >  2:1905574/1‑100 (MQ=255)           ggTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGc                                                                                    >  1:632274/1‑100 (MQ=255)           ggTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGc                                                                                    >  1:50238/1‑100 (MQ=255)              cTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCAtt                                                                                 >  1:2091842/1‑100 (MQ=255)                aTCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCAtttt                                                                               <  1:19870/100‑1 (MQ=255)                 tCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTc                                                                              >  2:1818631/1‑100 (MQ=255)                  cATGTACGACACGTTAAAAGCCGGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCGGGTAAAGTGAGCATCATCGTCAGCATTTTcc                                                                             >  2:527174/1‑100 (MQ=255)                  cATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTcc                                                                             <  1:1334452/100‑1 (MQ=255)                  cATGTACCACACGTTAAAAGCCTGATTACCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTcc                                                                             <  2:113665/100‑1 (MQ=255)                   aTGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCACCGTCTGACAGGCATCGTTCTGCGTACTGGTACAGTGAGCATCACCGGCAGCATTTTCCt                                                                            >  2:1022521/1‑100 (MQ=255)                    tGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTg                                                                           <  1:2110221/100‑1 (MQ=255)                    tGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTg                                                                           <  1:956702/100‑1 (MQ=255)                     gTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGa                                                                          <  2:2385970/100‑1 (MQ=255)                      tACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAg                                                                         <  1:387393/100‑1 (MQ=255)                       aCCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAgc                                                                        <  1:706254/100‑1 (MQ=255)                       aCCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAgc                                                                        >  1:1987155/1‑100 (MQ=255)                         cacaCGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATTGTCAGCATTTTCCTGAgcgc                                                                      <  2:1568083/100‑1 (MQ=255)                         cacaCGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAgcgc                                                                      >  2:1312917/1‑100 (MQ=255)                         cacaCGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAgcgc                                                                      >  2:1321364/1‑100 (MQ=255)                         cacaCGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAgcgc                                                                      <  1:185889/100‑1 (MQ=255)                          acacGTTAAGAGCTTGATTAGCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTTCGTGCTGGTAAAGGTAGCATCATCGTCAGCATTTTCTTGAgcg                                                                       >  2:756343/1‑98 (MQ=255)                          acacGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCt                                                                     >  2:2047681/1‑100 (MQ=255)                          acacGTTAAAAGCATGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATAATCGTCAGCATTTTCCTGAGCGCt                                                                     >  2:867176/1‑100 (MQ=255)                           cacGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCt                                                                     >  1:2264291/1‑99 (MQ=255)                             cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTc                                                                  <  2:632267/100‑1 (MQ=255)                             cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTc                                                                  >  2:68023/1‑100 (MQ=255)                             cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTc                                                                  <  1:1011617/100‑1 (MQ=255)                             cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTGCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTc                                                                  <  2:656494/100‑1 (MQ=255)                              gTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCa                                                                 >  1:1599588/1‑100 (MQ=255)                               ttAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCaa                                                                >  2:992790/1‑100 (MQ=255)                                tAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAAc                                                               <  1:1778008/100‑1 (MQ=255)                                 aaaaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACa                                                              >  2:1573367/1‑100 (MQ=255)                                   aaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAg                                                             >  1:224640/1‑99 (MQ=255)                                   aaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGc                                                            <  1:939535/100‑1 (MQ=255)                                    aGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCt                                                           >  1:688970/1‑100 (MQ=255)                                    aGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCt                                                           >  2:944415/1‑100 (MQ=255)                                        tGATTATCCAGACCTGGCAGCGTCTGAAAGGCCTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTtctc                                                       <  2:681016/100‑1 (MQ=255)                                        tGATTAGCCAGCCCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGc                                                            <  2:2417688/95‑1 (MQ=255)                                         gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGCAAAGTGAGCATCATGGTCAGCATTGTCCTGAGCGCTTTCAACAGCTtct                                                        >  2:149185/1‑98 (MQ=255)                                          aTTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGc                                                     >  2:730841/1‑100 (MQ=255)                                             aTCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCGGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                  >  1:1797977/1‑100 (MQ=255)                                             aTCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                  >  1:1955797/1‑100 (MQ=255)                                             aTCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGa                                                  >  2:1234261/1‑100 (MQ=255)                                               ccAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGc                                                >  2:1821550/1‑100 (MQ=255)                                               ccAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGc                                                >  1:1391579/1‑100 (MQ=255)                                                cAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACATCTTCTCGCTGAGCg                                               >  1:889094/1‑100 (MQ=255)                                                 aGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCgg                                              >  1:423215/1‑100 (MQ=255)                                                    ccTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCACCCGCTTCTCGCTGAGCGGGaa                                           >  2:2031227/1‑100 (MQ=255)                                                    ccTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGaa                                           <  2:167711/100‑1 (MQ=255)                                                    ccTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGaa                                           >  2:77415/1‑100 (MQ=255)                                                       ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAt                                        >  2:1886758/1‑100 (MQ=255)                                                       ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCCCGCTCAGCGAGaagc                                         <  1:1139046/99‑3 (MQ=255)                                                          aGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCCCGCTCAGCGAGaagctgtt                                     <  2:7595/100‑7 (MQ=255)                                                           gCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcgtc                                    <  1:1904503/100‑1 (MQ=255)                                                           gCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcgtc                                    <  2:1022728/100‑1 (MQ=255)                                                           gCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcgtc                                    <  1:2255555/100‑1 (MQ=255)                                                            cGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtcgtcg                                   >  1:1853567/1‑100 (MQ=255)                                                               cTGAACGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCCTCATCGTCAGCATTTTCCTGAGCGCTTGCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATa                                <  2:851209/100‑1 (MQ=255)                                                                tGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATaa                               >  2:1504670/1‑100 (MQ=255)                                                                tGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATaa                               >  1:361043/1‑100 (MQ=255)                                                                 gAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAAc                              <  2:2206566/100‑1 (MQ=255)                                                                       cTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAAc                        <  1:2093312/100‑1 (MQ=255)                                                                            ttCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTc                   <  2:589187/100‑1 (MQ=255)                                                                               tGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGgc                >  2:425835/1‑100 (MQ=255)                                                                                  gTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCt             <  1:2123489/100‑1 (MQ=255)                                                                                   tGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTg            >  1:744747/1‑100 (MQ=255)                                                                                    gCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTgg           <  2:2228176/100‑1 (MQ=255)                                                                                     cTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc          >  1:709436/1‑100 (MQ=255)                                                                                      tGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg         >  1:460039/1‑100 (MQ=255)                                                                                         tAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtcc      <  1:930904/100‑1 (MQ=255)                                                                                         tAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtcc      >  2:384156/1‑100 (MQ=255)                                                                                            aGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgtc   <  1:881882/100‑1 (MQ=255)                                                                                             gTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCCCTGGCGTCCGTCa  >  2:368370/1‑100 (MQ=255)                                                                                                 |                                                                                                TGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCA  >  NC_000913/3541866‑3542057

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