Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,530,925	T→C	Q103Q (CAA→CAG)	zipA ←	FtsZ stabilizer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,530,925	0	T	C	97.7%	303.0 / ‑2.2	86	Q103Q (CAA→CAG)	zipA	FtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (46/38);  total (47/39)
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.68e-01

CTGTTCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATGAT  >  NC_000913/2530826‑2531019                                                                                                    |                                                                                               cTGTTCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGt                                                                                                <  1:1360600/100‑1 (MQ=255)    ttCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGa                                                                                             <  1:1019731/100‑1 (MQ=255)     tCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGAc                                                                                            >  1:346874/1‑100 (MQ=255)      cAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGa                                                                                             <  1:859162/98‑1 (MQ=255)      cAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGAcg                                                                                           >  1:719865/1‑100 (MQ=255)       aGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGAcgc                                                                                          >  1:1878803/1‑100 (MQ=255)        ggCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGAcgcg                                                                                         >  1:1388781/1‑100 (MQ=255)          cTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGc                                                                                       <  1:1158303/100‑1 (MQ=255)           tGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCt                                                                                      >  2:1218044/1‑100 (MQ=255)            gCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTg                                                                                     <  2:54053/100‑1 (MQ=255)              ggCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCt                                                                                      >  1:1034773/1‑97 (MQ=255)               gCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTgcgc                                                                                  <  2:376002/100‑1 (MQ=255)                  gATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGa                                                                               <  2:586654/100‑1 (MQ=255)                   aTAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGAc                                                                              <  2:1106893/100‑1 (MQ=255)                    tAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACg                                                                             <  1:413726/100‑1 (MQ=255)                    tAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACg                                                                             >  1:312263/1‑100 (MQ=255)                     aggcaggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGc                                                                            >  2:1747385/1‑100 (MQ=255)                      ggcaggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCa                                                                           >  1:1328541/1‑100 (MQ=255)                      ggcaggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCa                                                                           >  2:60664/1‑100 (MQ=255)                      ggcaggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTACGCAGACGCa                                                                           >  1:1280023/1‑100 (MQ=255)                        caggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATa                                                                         <  1:1478006/100‑1 (MQ=255)                          ggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAg                                                                       <  2:638096/100‑1 (MQ=255)                          ggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAg                                                                       <  2:1275809/100‑1 (MQ=255)                          ggcTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAg                                                                       >  1:1399787/1‑100 (MQ=255)                              gctgcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCgg                                                                   <  2:978118/100‑1 (MQ=255)                              gctgcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCg                                                                    >  1:459665/1‑99 (MQ=255)                                tgcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGtt                                                                 >  2:907457/1‑100 (MQ=255)                                 gcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTg                                                                >  1:414377/1‑100 (MQ=255)                                 gcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTg                                                                >  2:1792418/1‑100 (MQ=255)                                 gcACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTg                                                                >  2:1791362/1‑100 (MQ=255)                                    ccGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTa                                                             >  1:1185982/1‑100 (MQ=255)                                     cGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTAc                                                            <  2:1884293/100‑1 (MQ=255)                                      ggCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACt                                                           <  2:1878815/100‑1 (MQ=255)                                      ggCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACt                                                           >  2:1851865/1‑100 (MQ=255)                                        cTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTgg                                                         >  2:1506611/1‑100 (MQ=255)                                          gcgcTTGATGCGGAGCATGCTGCGGCGGTCCCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCCGACGCATAAGGCGGTTGGTACTGgtg                                                       >  1:395169/1‑100 (MQ=255)                                          gcgcTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGgtg                                                       <  1:506785/100‑1 (MQ=255)                                              tGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgc                                                   >  1:1190345/1‑100 (MQ=255)                                               ggATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcg                                                  >  1:262931/1‑100 (MQ=255)                                               ggATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcg                                                  >  2:137908/1‑100 (MQ=255)                                               ggATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcg                                                  >  2:1230892/1‑100 (MQ=255)                                                 aTGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcggc                                                <  2:1185984/100‑1 (MQ=255)                                                  tGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcggcg                                               >  2:931507/1‑100 (MQ=255)                                                  tGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcggcg                                               <  1:127090/100‑1 (MQ=255)                                                  tGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTgcggcg                                               <  2:1328549/100‑1 (MQ=255)                                                     ggAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgacg                                            >  1:804930/1‑100 (MQ=255)                                                      gAGCATGTTGCGGCGGTACCTGCGCTTGAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgacgg                                           <  1:697118/100‑1 (MQ=255)                                                        gCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgacggac                                         >  2:260216/1‑100 (MQ=255)                                                        gCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgacgga                                          <  1:1816779/99‑1 (MQ=255)                                                         cATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgacggacg                                        <  2:1376610/100‑1 (MQ=255)                                                         cATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCgacggacg                                        <  2:1631581/100‑1 (MQ=255)                                                          atcttGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGa                                       <  1:576034/97‑1 (MQ=255)                                                           tGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGag                                      >  2:326187/1‑100 (MQ=255)                                                           tGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGag                                      <  2:915829/100‑1 (MQ=255)                                                            gTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagc                                     <  2:879896/100‑1 (MQ=255)                                                            gTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagc                                     >  2:679849/1‑100 (MQ=255)                                                            gTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagc                                     >  1:1899196/1‑100 (MQ=255)                                                              tGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagcag                                   >  2:804089/1‑100 (MQ=255)                                                              tGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGaagat                                   >  2:576036/1‑96 (MQ=255)                                                               gcggcgGTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGagcagc                                  >  1:32509/1‑100 (MQ=255)                                                                cggcggTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGcc                                 <  2:414379/100‑1 (MQ=255)                                                                cggcggTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGcc                                 >  2:760193/1‑100 (MQ=255)                                                                   cggTACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCa                              <  2:1280030/100‑1 (MQ=255)                                                                     gTACCTGCGCTTCAGGCGGCTGCTGGACCCGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGTACGAGCAGCCTCATg                            <  1:1603071/100‑1 (MQ=255)                                                                      tACCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGc                           <  2:312265/100‑1 (MQ=255)                                                                       aCCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCt                          >  1:239581/1‑100 (MQ=255)                                                                       aCCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCt                          >  1:1589674/1‑100 (MQ=255)                                                                       aCCTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCt                          >  2:407980/1‑100 (MQ=255)                                                                         cTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTcc                        <  1:1747374/100‑1 (MQ=255)                                                                         cTGCGCTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTcc                        >  2:1443949/1‑100 (MQ=255)                                                                           gcgcTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTg                      >  2:1294206/1‑100 (MQ=255)                                                                            cgcTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGa                     <  1:684461/100‑1 (MQ=255)                                                                             gcTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAg                    <  1:1083908/100‑1 (MQ=255)                                                                              cTTCAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGc                   >  2:1297560/1‑100 (MQ=255)                                                                              cTTCAGGAGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGc                   >  2:1027133/1‑100 (MQ=255)                                                                                 cAGGCGGCTGCTGGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTGGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGtt                >  1:1209743/1‑100 (MQ=255)                                                                                        ctgctgGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCggg         >  1:1597058/1‑100 (MQ=255)                                                                                        ctgctgGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCggg         <  1:182675/100‑1 (MQ=255)                                                                                        ctgctgGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCggg         <  2:32509/100‑1 (MQ=255)                                                                                          gctgGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGc       <  2:1899209/100‑1 (MQ=255)                                                                                          gctgGACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGc       <  2:1541376/100‑1 (MQ=255)                                                                                             ggACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATg    >  2:759264/1‑100 (MQ=255)                                                                                             ggACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATg    <  2:346876/100‑1 (MQ=255)                                                                                             ggACCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATg    <  2:223191/100‑1 (MQ=255)                                                                                               aCCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATGAt  <  2:759235/100‑1 (MQ=255)                                                                                               aCCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATGAt  >  2:1512861/1‑100 (MQ=255)                                                                                               aCCGGCTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATGAt  >  2:462361/1‑100 (MQ=255)                                                                                                    |                                                                                               CTGTTCAGGCTGCGGCTGATAGGCAGGCTGCTGCACCGGCTGCGCTGGATGCGGAGCATGTTGCGGCGGTACCTGCGCTTCAGGCGGCTGCTGGACCGGTTGACGCGGCTGCGCAGACGCATAAGGCGGTTGGTACTGGTGTTGCGGCGACGGACGAGCAGCCTCATGCTCCTGAGCGTTAGCCGGGGCATGAT  >  NC_000913/2530826‑2531019

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