Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,386,922	Δ1 bp	coding (2449/3324 nt)	mscM ←	miniconductance mechanosensitive channel MscM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,386,922	0	T	.	100.0%	370.2 / NA	85	coding (2449/3324 nt)	mscM	miniconductance mechanosensitive channel MscM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (46/39);  total (46/39)

CTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCAT  >  NC_000913/4386788‑4387056                                                                                                                                       |                                                                                                                                       ctgcAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATc                                                                                                                                  >  2:555306/1‑140 (MQ=255)   gcAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATccc                                                                                                                                >  2:680143/1‑140 (MQ=255)   gcAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATccc                                                                                                                                >  2:645029/1‑140 (MQ=255)    cAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCca                                                                                                                               >  2:923500/1‑140 (MQ=255)    cAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCca                                                                                                                               <  2:829160/140‑1 (MQ=255)     aaaaTCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCcac                                                                                                                              <  1:49398/140‑1 (MQ=255)         tCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCg                                                                                                                          >  2:137146/1‑140 (MQ=255)         tCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCg                                                                                                                          >  1:1060358/1‑140 (MQ=255)            ccAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGaaa                                                                                                                       <  2:383814/140‑1 (MQ=255)              aGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAata                                                                                                                     >  1:717226/1‑140 (MQ=255)                         cgcgGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGaa                                                                                                          >  1:490929/1‑140 (MQ=255)                              gCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGa                                                                                                     <  1:555306/140‑1 (MQ=255)                                 aGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAg                                                                                                  <  2:660605/140‑1 (MQ=255)                                   ttGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCg                                                                                                <  2:108398/140‑1 (MQ=255)                                        cACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATg                                                                                           <  1:763157/140‑1 (MQ=255)                                        cACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATg                                                                                           >  2:163299/1‑140 (MQ=255)                                        cACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATg                                                                                           >  1:361234/1‑140 (MQ=255)                                        cACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATg                                                                                           >  2:901565/1‑140 (MQ=255)                                        cACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATg                                                                                           >  1:114991/1‑140 (MQ=255)                                        cACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATg                                                                                           <  2:990556/140‑1 (MQ=255)                                        cACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATg                                                                                           <  1:104833/140‑1 (MQ=255)                                        cACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GAC‑TCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATgg                                                                                          >  2:46920/1‑140 (MQ=255)                                          ccAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGa                                                                                         <  1:901565/140‑1 (MQ=255)                                             gCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAACGCGGAATGGAttt                                                                                      >  2:805056/1‑140 (MQ=255)                                             gCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGAttt                                                                                      >  2:1029081/1‑140 (MQ=255)                                               tGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCt                                                                                    <  2:114991/140‑1 (MQ=255)                                                   tggtgATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAcc                                                                                >  1:762271/1‑140 (MQ=255)                                                   tggtgATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAcc                                                                                <  2:1235757/140‑1 (MQ=255)                                                       gatgatAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAg                                                                            <  1:1029081/140‑1 (MQ=255)                                                        atgatAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGc                                                                           >  1:722656/1‑140 (MQ=255)                                                         tgatAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGca                                                                          >  2:486924/1‑140 (MQ=255)                                                         tgatAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GAC‑TCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGcac                                                                         <  1:46920/140‑1 (MQ=255)                                                          gatAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGcac                                                                         <  2:717226/140‑1 (MQ=255)                                                           atAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGcaca                                                                        >  1:1040587/1‑140 (MQ=255)                                                           atAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGcaca                                                                        >  1:426914/1‑140 (MQ=255)                                                            tAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACaa                                                                       <  1:163299/140‑1 (MQ=255)                                                               acacCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGa                                                                    <  1:680143/140‑1 (MQ=255)                                                                cacCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGAc                                                                   >  1:732461/1‑140 (MQ=255)                                                                 acCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACa                                                                  <  1:672032/140‑1 (MQ=255)                                                                  ccAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAc                                                                                 >  1:596026/1‑124 (MQ=255)                                                                  ccAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGAc                                                                                 <  2:596026/124‑1 (MQ=255)                                                                  ccAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACaa                                                                       >  2:934231/1‑134 (MQ=255)                                                                  ccAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACaa                                                                       <  1:934231/134‑1 (MQ=255)                                                                  ccAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAg                                                                 >  2:873954/1‑140 (MQ=255)                                                                   cAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGa                                                                <  2:762271/140‑1 (MQ=255)                                                                      aaTGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAg                                                             >  2:536837/1‑140 (MQ=255)                                                                         ggCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACaa                                                                       >  1:174553/1‑127 (MQ=255)                                                                         ggCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACaa                                                                       <  2:174553/127‑1 (MQ=255)                                                                           cAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAgag                                                        >  1:285124/1‑138 (MQ=255)                                                                             aTCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAggg                                                        <  1:427206/138‑1 (MQ=255)                                                                             aTCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAggg                                                        >  2:427206/1‑138 (MQ=255)                                                                                aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATc                                                   <  2:722656/140‑1 (MQ=255)                                                                                aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATc                                                   <  2:843530/140‑1 (MQ=255)                                                                                aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATc                                                   >  1:967803/1‑140 (MQ=255)                                                                                  caccgcaccgAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCaa                                                 >  1:531424/1‑140 (MQ=255)                                                                                       caccgAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCTCAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCa                                            >  2:495710/1‑140 (MQ=255)                                                                                       caccgAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCa                                            >  1:1127514/1‑140 (MQ=255)                                                                                       caccgAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCa                                            >  1:634898/1‑140 (MQ=255)                                                                                          cgAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCaaaa                                         <  2:426914/140‑1 (MQ=255)                                                                                          cgAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCaaaa                                         <  2:1040587/140‑1 (MQ=255)                                                                                            aGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAAtt                                       >  2:363069/1‑140 (MQ=255)                                                                                            aGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAAtt                                       >  2:706330/1‑140 (MQ=255)                                                                                               gTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCTCAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAg                                    <  1:495710/140‑1 (MQ=255)                                                                                                  aTTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAgcgc                                 >  1:1142894/1‑140 (MQ=255)                                                                                                    tGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCAc                               <  1:873954/140‑1 (MQ=255)                                                                                                      gCTCCAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCa                             <  1:265228/140‑1 (MQ=255)                                                                                                         ccAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGcc                          >  2:72419/1‑140 (MQ=255)                                                                                                         ccAGACTTTCTAAGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAg                                    >  2:171460/1‑130 (MQ=255)                                                                                                         ccAGACTTTCTAAGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAg                                    <  1:171460/130‑1 (MQ=255)                                                                                                          cAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCg                         >  2:70510/1‑140 (MQ=255)                                                                                                          cAGACTTTCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCg                         >  2:712333/1‑140 (MQ=255)                                                                                                                ttCTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGa                   >  2:302793/1‑140 (MQ=255)                                                                                                                  cTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCa                             <  1:234263/128‑1 (MQ=255)                                                                                                                  cTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCa                             >  2:234263/1‑128 (MQ=255)                                                                                                                  cTACGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGAtt                 <  1:678485/140‑1 (MQ=255)                                                                                                                     cGCCCTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTgcg              <  1:923500/140‑1 (MQ=255)                                                                                                                        ccTGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCAc           <  1:712333/140‑1 (MQ=255)                                                                                                                          tGTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTg         <  2:1142894/140‑1 (MQ=255)                                                                                                                           gTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGa                   <  1:1238425/129‑1 (MQ=255)                                                                                                                           gTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGa                   >  2:1238425/1‑129 (MQ=255)                                                                                                                           gTACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGa        >  2:1233011/1‑140 (MQ=255)                                                                                                                            tACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGAt       >  2:581388/1‑140 (MQ=255)                                                                                                                            tACCGTGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGAt       <  2:1165193/140‑1 (MQ=255)                                                                                                                             aCCGTGGAGG‑GACATCCCCCAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATg      <  2:285124/140‑1 (MQ=255)                                                                                                                                 tGGAGG‑GACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCAt  >  2:952039/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                       CTGCAAAATCGCCAGTTCCAGCAGCGCGGGCAAGTTGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGAAGCCGAAAGCGGAATGGATTTCTGACCACAGCACAATGACAGAAAGCAGGGCGATCAACATCAAAATTGAGCGCACCAGCCGCAAGGATTGCGCACTGATGGCAT  >  NC_000913/4386788‑4387056

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: ‑