Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,453,508	G→A	pseudogene (1494/2646 nt)	yfcU ←	putative outer membrane usher protein YfcU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,453,508	0	G	A	100.0%	232.1 / NA	67	pseudogene (1494/2646 nt)	yfcU	putative outer membrane usher protein YfcU
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/39);  total (28/39)

GTAGTTTGTCTGCTCCTCACGATCCCAGTAGGTATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCG  >  NC_000913/2453373‑2453637                                                                                                                                        |                                                                                                                                  gtAGTTTGTCTGCTCCTCACGATCCCAGTAGGTATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGt                                                                                                                               <  2:702662/140‑1 (MQ=255)       tGTCTGCTCCTCACGATCCCAGTAGGTATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTgg                                                                                                                         <  1:61013/140‑1 (MQ=255)       tGTCTGCTCCTCACGATCCCAGTAGGTATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTgg                                                                                                                         <  1:717835/140‑1 (MQ=255)       tGTCTGCTCCTCACGATCCCAGTAGGTATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTgg                                                                                                                         <  1:134653/140‑1 (MQ=255)               cctcACGATCCCAGTAGGTATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGa                                                                                                                 <  2:522520/140‑1 (MQ=255)                  cACGATCCCAGTAGGTATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATAc                                                                                                              <  1:1145397/140‑1 (MQ=255)                         ccAGTAGGTATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCa                                                                                                       <  1:1194618/140‑1 (MQ=255)                            gTAGGTATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATgg                                                                                                    <  1:287675/140‑1 (MQ=255)                              aGGTATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTc                                                                                                  <  1:397921/140‑1 (MQ=255)                                 tATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATa                                                                                               <  1:891476/140‑1 (MQ=255)                                      cGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGtt                                                                                          <  1:324019/140‑1 (MQ=255)                                         gTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTtctc                                                                                       <  1:370227/140‑1 (MQ=255)                                               ttGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGAGTACATCTCTTTGTCGTTGCCCGAACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCga                                                                                 >  2:484877/1‑140 (MQ=255)                                                    aTAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGc                                                                            >  1:1197882/1‑140 (MQ=255)                                                    aTAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGc                                                                            >  1:202398/1‑140 (MQ=255)                                                    aTAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGc                                                                            >  2:728028/1‑140 (MQ=255)                                                        aCCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGata                                                                        >  2:1126677/1‑140 (MQ=255)                                                        aCCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGata                                                                        >  1:703100/1‑140 (MQ=255)                                                         ccGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCgag                                                                                <  1:45289/131‑1 (MQ=255)                                                         ccGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCgag                                                                                >  2:45289/1‑131 (MQ=255)                                                         ccGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGatat                                                                       >  2:411696/1‑140 (MQ=255)                                                                 cccGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGaa                                                               <  1:1170777/140‑1 (MQ=255)                                                                      aTCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGa                                                          <  2:202398/140‑1 (MQ=255)                                                                      aTCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGa                                                          <  1:411696/140‑1 (MQ=255)                                                                         gcgGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGAcgc                                                       <  2:1197882/140‑1 (MQ=255)                                                                          cgGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGAcgcg                                                      >  1:149554/1‑140 (MQ=255)                                                                           ggAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCgg                                                     <  1:1072970/140‑1 (MQ=255)                                                                           ggAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGCCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCgg                                                     <  1:261456/140‑1 (MQ=255)                                                                            gAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGc                                                    <  1:1379011/140‑1 (MQ=255)                                                                                  cTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCACTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGa                                              >  1:1181424/1‑140 (MQ=255)                                                                                   tGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAg                                             <  2:1299269/140‑1 (MQ=255)                                                                                    ggTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGc                                            <  1:704783/140‑1 (MQ=255)                                                                                          aaGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCCTAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCa                                      >  2:375156/1‑140 (MQ=255)                                                                                          aaGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGC‑GCTGTTGAGCTGGTCaa                                     >  1:315289/1‑140 (MQ=255)                                                                                               gCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCaaa                                    >  1:75809/1‑137 (MQ=255)                                                                                               gCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCaaa                                    <  2:75809/137‑1 (MQ=255)                                                                                               gCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTc                                 >  2:618567/1‑140 (MQ=255)                                                                                                cGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGc                                            <  2:793518/128‑1 (MQ=255)                                                                                                cGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGc                                            >  1:793518/1‑128 (MQ=255)                                                                                                cGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCt                                >  2:37207/1‑140 (MQ=255)                                                                                                 ggTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCtt                               <  2:1181424/140‑1 (MQ=255)                                                                                                 ggTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCtt                               <  2:149554/140‑1 (MQ=255)                                                                                                   tgtACATCTCTTTGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTg                             >  1:408577/1‑140 (MQ=255)                                                                                                             tttGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCtt                               <  1:1069667/128‑1 (MQ=255)                                                                                                             tttGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCtt                               >  2:1069667/1‑128 (MQ=255)                                                                                                             tttGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACt                     >  1:1478432/1‑138 (MQ=255)                                                                                                             tttGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACt                     <  2:1478432/138‑1 (MQ=255)                                                                                                             tttGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTca                   <  2:408577/140‑1 (MQ=255)                                                                                                             tttGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTca                   >  2:943476/1‑140 (MQ=255)                                                                                                             tttGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTca                   <  2:29642/140‑1 (MQ=255)                                                                                                             tttGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTca                   >  1:1254814/1‑140 (MQ=255)                                                                                                             tttGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTca                   <  1:24991/140‑1 (MQ=255)                                                                                                             tttGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTca                   <  2:1175661/140‑1 (MQ=255)                                                                                                             tttGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTca                   <  2:1078326/140‑1 (MQ=255)                                                                                                               tGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTcaca                 >  1:251128/1‑140 (MQ=255)                                                                                                               tGTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGC‑GCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCacac                <  2:315289/140‑1 (MQ=255)                                                                                                                gTCGTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCacac                <  2:310235/140‑1 (MQ=255)                                                                                                                   gTTGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGa             >  2:1451141/1‑140 (MQ=255)                                                                                                                    ttGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGaa            >  1:40289/1‑140 (MQ=255)                                                                                                                    ttGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGaa            >  1:913514/1‑140 (MQ=255)                                                                                                                     tGCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAg           <  1:943476/140‑1 (MQ=255)                                                                                                                      gCCCGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAgg          <  2:913514/140‑1 (MQ=255)                                                                                                                       cccGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATaa                       >  2:122851/1‑126 (MQ=255)                                                                                                                       cccGTACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATaa                       <  1:122851/126‑1 (MQ=255)                                                                                                                           tACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTca                   >  1:171771/1‑126 (MQ=255)                                                                                                                           tACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTca                   <  2:171771/126‑1 (MQ=255)                                                                                                                           tACGGACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTa     >  2:441310/1‑140 (MQ=255)                                                                                                                              ggACCATTTCACTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCg  <  2:251128/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  GTAGTTTGTCTGCTCCTCACGATCCCAGTAGGTATGGCGGGTATAGTTGAGATAAACCGAAACACCCGCATCGCGGAAGTTCTGGTTATAAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCG  >  NC_000913/2453373‑2453637

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