Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,531,596	(G)6→5	coding (134/762 nt)	cysZ →	sulfate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,531,591	0	G	.	100.0%	264.3 / NA	63	coding (129/762 nt)	cysZ	sulfate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (26/37);  total (26/37)

CATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTT  >  NC_000913/2531453‑2531726                                                                                                                                           |                                                                                                                                        cATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgg                                                                                                                                        >  2:435659/1‑140 (MQ=255)    ccTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATggggg                                                                                                                                     >  1:856537/1‑140 (MQ=255)    ccTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATggggg                                                                                                                                     <  1:17458/140‑1 (MQ=255)      tGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggggcg                                                                                                                                   <  1:141661/140‑3 (MQ=255)        aTTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggggcgca                                                                                                                                 <  1:435659/140‑5 (MQ=255)          tATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCAtt                                                                                                                              <  1:996841/140‑1 (MQ=255)              gTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtgg                                                                                                                          >  1:15498/1‑140 (MQ=255)              gTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtgg                                                                                                                          >  1:788802/1‑140 (MQ=255)                ttCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtg                                                                                                                        <  2:788802/140‑1 (MQ=255)                ttCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtg                                                                                                                        <  1:445932/140‑1 (MQ=255)                  catcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGc                                                                                                                      >  1:1065072/1‑140 (MQ=255)                  catcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGc                                                                                                                      >  1:275653/1‑140 (MQ=255)                  catcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGc                                                                                                                      >  2:328557/1‑140 (MQ=255)                  catcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGc                                                                                                                      >  1:384865/1‑140 (MQ=255)                  catcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTAGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGc                                                                                                                      >  2:800340/1‑140 (MQ=255)                   atcatTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGct                                                                                                                     <  2:1120352/140‑1 (MQ=255)                         cacaTCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTac                                                                                                               >  1:568455/1‑140 (MQ=255)                          acaTCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCCTTaca                                                                                                              <  2:253537/140‑1 (MQ=255)                           caTCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTacac                                                                                                             <  1:328557/140‑1 (MQ=255)                               tGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtg                                                                                                                        <  2:150742/125‑1 (MQ=255)                               tGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtg                                                                                                                        >  1:150742/1‑125 (MQ=255)                               tGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCt                                                                                                         >  2:1060238/1‑140 (MQ=255)                               tGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCt                                                                                                         <  2:856537/140‑1 (MQ=255)                                 ccccACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtgg                                                                                                                       >  1:59691/1‑124 (MQ=255)                                 ccccACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCtggtgg                                                                                                                       <  2:59691/124‑1 (MQ=255)                                 ccccACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCg                                                                                                       >  2:234579/1‑140 (MQ=255)                                      cGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCt                                                                                                         >  2:765870/1‑133 (MQ=255)                                      cGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCt                                                                                                         <  1:765870/133‑1 (MQ=255)                                      cGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGa                                                                                                      >  2:121753/1‑136 (MQ=255)                                      cGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGa                                                                                                      <  1:121753/136‑1 (MQ=255)                                        cAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                                                                <  2:249132/140‑1 (MQ=255)                                           cGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGAt                                                                                             <  2:384865/140‑1 (MQ=255)                                           cGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGAt                                                                                             >  1:135721/1‑140 (MQ=255)                                                tttACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGa                                                                                        <  2:135721/140‑1 (MQ=255)                                                tttACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGa                                                                                        <  2:755492/140‑1 (MQ=255)                                                       ttttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCAt                                                                                 >  2:574771/1‑140 (MQ=255)                                                          tGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAg                                                                              <  1:647108/140‑1 (MQ=255)                                                          tGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAg                                                                              >  2:207309/1‑140 (MQ=255)                                                          tGCGCAAGGCTGGAAGCTCGACTCGCAACCTGGGCTTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTATCATGAg                                                                              >  2:442613/1‑140 (MQ=255)                                                            cgcAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGtt                                                                            <  2:100023/140‑1 (MQ=255)                                                                    tGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCg                                                                    >  2:528513/1‑140 (MQ=255)                                                                     ggAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCgg                                                                   <  1:343678/140‑1 (MQ=255)                                                                      gAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGa                                                                  <  1:1060238/140‑1 (MQ=255)                                                                         gCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTg                                                               <  1:413748/140‑1 (MQ=255)                                                                              tctcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTAc                                                          <  2:938963/140‑1 (MQ=255)                                                                                     aCCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCt                                                   <  1:100639/140‑1 (MQ=255)                                                                                     aCCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCt                                                   <  2:659839/140‑1 (MQ=255)                                                                                     aCCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCt                                                   <  2:668866/140‑1 (MQ=255)                                                                                            aTTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTAt                                            <  2:867243/140‑1 (MQ=255)                                                                                            aTTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTAt                                            <  2:503719/140‑1 (MQ=255)                                                                                                ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCtgt                                        >  2:1053751/1‑140 (MQ=255)                                                                                                ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCtgt                                        >  1:194493/1‑140 (MQ=255)                                                                                                         ttATTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCctc                               <  1:574771/140‑1 (MQ=255)                                                                                                           aTTTTACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTg                             <  1:207309/140‑1 (MQ=255)                                                                                                               tACCGCTGCTGGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTg                         >  1:318194/1‑140 (MQ=255)                                                                                                                    ctgctgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCAtc                    <  2:194493/140‑1 (MQ=255)                                                                                                                      gctgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGt                        <  2:330671/134‑1 (MQ=255)                                                                                                                      gctgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGt                        >  1:330671/1‑134 (MQ=255)                                                                                                                       ctgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATctct                 >  2:962828/1‑140 (MQ=255)                                                                                                                       ctgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATctct                 >  2:736984/1‑140 (MQ=255)                                                                                                                        tgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCtg                <  1:736984/140‑1 (MQ=255)                                                                                                                        tgGTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCtg                <  1:962828/140‑1 (MQ=255)                                                                                                                          gTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCtgtg              >  2:886635/1‑140 (MQ=255)                                                                                                                          gTCAATATTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCtgtg              <  2:318194/140‑1 (MQ=255)                                                                                                                              atatTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGt          <  2:686417/140‑1 (MQ=255)                                                                                                                              atatTTTGTTGAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGt          <  2:568455/140‑1 (MQ=255)                                                                                                                                    tgttgAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAgtgt    >  2:927631/1‑140 (MQ=255)                                                                                                                                     gttgAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGtt   <  1:528513/140‑1 (MQ=255)                                                                                                                                      ttgAT‑GGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGttt  <  1:234579/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                        CATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTT  >  NC_000913/2531453‑2531726

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