Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	523,311	T→C	G17G (GGT→GGC)	rhsD →	protein RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	523,311	0	T	C	100.0%	229.3 / NA	66	G17G (GGT→GGC)	rhsD	protein RhsD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (39/27);  total (39/27)

AATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGA  >  NC_000913/523176‑523444                                                                                                                                        |                                                                                                                                       aaTACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCa                                                                                                                                    <  2:903903/140‑1 (MQ=255) aaTACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCa                                                                                                                                    <  2:84913/140‑1 (MQ=255)       ttAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTcc                                                                                                                              >  1:370860/1‑140 (MQ=255)             ggTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGtt                                                                                                                        >  2:568315/1‑140 (MQ=255)              gTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTc                                                                                                                       >  1:248885/1‑140 (MQ=255)               tCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCg                                                                                                                      >  1:60103/1‑140 (MQ=255)                 ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGc                                                                                                                    <  1:568315/140‑1 (MQ=255)                 ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGc                                                                                                                    <  2:189231/140‑1 (MQ=255)                 ttGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGc                                                                                                                    <  2:60103/140‑1 (MQ=255)                    tgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCgg                                                                                                                     >  1:442224/1‑136 (MQ=255)                    tgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCgg                                                                                                                     <  2:442224/136‑1 (MQ=255)                    tgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGc                                                                                                                    <  1:102544/137‑1 (MQ=255)                    tgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGc                                                                                                                    >  2:102544/1‑137 (MQ=255)                    tgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAgg                                                                                                                 >  2:336461/1‑140 (MQ=255)                    tgtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAgg                                                                                                                 >  1:244377/1‑140 (MQ=255)                     gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCgg                                                                                                                     <  2:8047/135‑1 (MQ=255)                     gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCgg                                                                                                                     >  1:8047/1‑135 (MQ=255)                        tCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAGCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCATCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTa                                                                                                             >  1:617289/1‑140 (MQ=255)                              tattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAAtt                                                                                                       >  1:873080/1‑140 (MQ=255)                                                acacAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtg                                                                                     <  2:873080/140‑1 (MQ=255)                                                      ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtggcg                                                                                  <  2:1100356/137‑1 (MQ=255)                                                      ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtggcg                                                                                  >  1:1100356/1‑137 (MQ=255)                                                      ggATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGc                                                                               >  1:863608/1‑140 (MQ=255)                                                           tGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGgt                                                                          >  1:635517/1‑140 (MQ=255)                                                           tGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCggg                                                                          >  1:710202/1‑139 (MQ=255)                                                            gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCGGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCgg                                                                                         <  1:362564/124‑1 (MQ=255)                                                            gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCGGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCgg                                                                                         >  2:362564/1‑124 (MQ=255)                                                            gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGgtg                                                                         >  2:70595/1‑140 (MQ=255)                                                            gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCgggg                                                                         >  2:751261/1‑138 (MQ=255)                                                            gAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGgtg                                                                         >  2:463752/1‑140 (MQ=255)                                                             aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGgtgg                                                                        >  2:668186/1‑139 (MQ=255)                                                                aaTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTc                                                                     >  2:497703/1‑140 (MQ=255)                                                                  tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcg                                                                   >  1:78573/1‑140 (MQ=255)                                                                  tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcg                                                                   >  1:194823/1‑140 (MQ=255)                                                                    aCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcggg                                                                 <  1:70595/140‑1 (MQ=255)                                                                    aCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcggg                                                                 <  1:497703/140‑1 (MQ=255)                                                                     cTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCcgggg                                                                >  2:747272/1‑139 (MQ=255)                                                                          aGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGAt                                                          >  2:702371/1‑140 (MQ=255)                                                                          aGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGAt                                                          <  2:831420/140‑1 (MQ=255)                                                                            ggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGa                                                        >  1:866437/1‑140 (MQ=255)                                                                               aCTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTTGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACtt                                                     >  2:302500/1‑140 (MQ=255)                                                                                cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGG‑CGGGGTGACttt                                                    >  2:886256/1‑139 (MQ=255)                                                                                  ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCgg                                                  <  1:336461/140‑1 (MQ=255)                                                                                  ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCgg                                                  >  1:682929/1‑140 (MQ=255)                                                                                    atatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGc                                                <  1:581171/140‑1 (MQ=255)                                                                                    atatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGc                                                <  2:514302/140‑1 (MQ=255)                                                                                           cGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCaa                                              >  1:74507/1‑135 (MQ=255)                                                                                           cGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCaa                                              <  2:74507/135‑1 (MQ=255)                                                                                           cGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCgg                                         >  1:551132/1‑140 (MQ=255)                                                                                           cGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCgg                                         >  1:174675/1‑140 (MQ=255)                                                                                             gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCTCGGGGGGTCCGGG‑CGGGGTGACTTCGGGCAACCCggga                                       >  2:552437/1‑138 (MQ=255)                                                                                             gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCAATGCCCAGGGTTCGGCAGGCGGAAGAAATGGCGCGCCCCCCCGCGAGGCGTGCTCGGGGTGTGCGGGACGGGTGAATTTGGGGCACCCCgga                                        >  2:847313/1‑139 (MQ=255)                                                                                                   cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCg                                 <  2:533103/140‑1 (MQ=255)                                                                                                   cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCg                                 <  1:300231/140‑1 (MQ=255)                                                                                                      cggcgCGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCgctg                              <  2:244377/140‑1 (MQ=255)                                                                                                        gcgcGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCgctgct                            <  1:702371/140‑1 (MQ=255)                                                                                                        gcgcGTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCgctgct                            <  2:194823/140‑1 (MQ=255)                                                                                                           cgTCAGGGAGATATGACTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTggg                         <  2:551132/140‑1 (MQ=255)                                                                                                                      tatGACTCAGTATGGCGGCCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGG‑CGGGATGACCTCCGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgc              >  2:190239/1‑140 (MQ=255)                                                                                                                          aCTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGGTGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGGGCTGcc          >  2:359299/1‑140 (MQ=255)                                                                                                                           cTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGccc         >  1:998393/1‑140 (MQ=255)                                                                                                                           cTCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGccc         >  1:863811/1‑140 (MQ=255)                                                                                                                            tCAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCg        <  2:863608/140‑1 (MQ=255)                                                                                                                             cAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCgg       <  1:886256/140‑1 (MQ=255)                                                                                                                             cAGTATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCgg       >  1:443388/1‑140 (MQ=255)                                                                                                                                tATGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCga    >  2:76745/1‑140 (MQ=255)                                                                                                                                  tGGCGGCCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgaga  <  2:866437/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       AATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGG‑CGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGA  >  NC_000913/523176‑523444

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