Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,659,608	C→T	5.9%	P93P (CCG→CCA)	iscA ←	FeS cluster assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,659,608	0	C	T	5.9%	223.1 / 2.5	68	P93P (CCG→CCA)	iscA	FeS cluster assembly protein
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base T (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAACACGATGTCTTCCGGCGTCGG  >  NC_000913/2659468‑2659745                                                                                                                                             |                                                                                                                                          aaCATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTgg                                                                                                                                    >  4:1652/1‑148 (MQ=255) aaCATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTgg                                                                                                                                    >  2:100633/1‑148 (MQ=255)  aCATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGt                                                                                                                                   >  6:85305/1‑148 (MQ=255)   cATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGt                                                                                                                                   >  4:12261/1‑147 (MQ=255)     ttCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGaa                                                                                                                                >  2:195720/1‑148 (MQ=255)          aGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAAttt                                                                                                                             <  1:324101/146‑1 (MQ=255)          aGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTa                                                                                                                            <  4:132556/147‑1 (MQ=255)           gggCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTTGGGTTGGTGAATTTaaa                                                                                                                          >  8:89481/1‑148 (MQ=255)                gCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCtt                                                                                                                     <  7:29760/148‑1 (MQ=255)                gCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCtt                                                                                                                     >  5:385804/1‑148 (MQ=255)                    tCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGtt                                                                                                                 >  4:159290/1‑148 (MQ=255)                      cgccgcGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTTGGGTTGGTGAATTTAAACCCTTCGTTc                                                                                                                <  7:89481/147‑1 (MQ=255)                      cgccgcGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTc                                                                                                                <  1:195720/147‑1 (MQ=255)                      cgccgcGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTc                                                                                                                <  1:203299/147‑1 (MQ=255)                      cgccgcGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTc                                                                                                                <  3:1652/147‑1 (MQ=255)                      cgccgcGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTc                                                                                                                <  3:12261/147‑1 (MQ=255)                               aaCGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAAttt                                                                                                                             <  1:49121/125‑1 (MQ=255)                               aaCGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAAttt                                                                                                                             >  2:49121/1‑125 (MQ=255)                                 cGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCtttt                                                                                                    <  3:159290/148‑1 (MQ=255)                                 cGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCttt                                                                                                     >  7:230789/1‑147 (MQ=255)                                     gAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGa                                                                                                <  1:100633/148‑1 (MQ=255)                                           ggTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTTGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCa                                                                                          <  8:230789/148‑1 (MQ=255)                                             tAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGc                                                                                        >  5:367220/1‑148 (MQ=255)                                              aaaaCAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTaa                                                                                                                           >  3:118954/1‑112 (MQ=255)                                              aaaaCAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTaa                                                                                                                           <  4:118954/112‑1 (MQ=255)                                               aaacaaacCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTg                                                                                      <  6:385804/148‑1 (MQ=255)                                                aacaaacCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCtttt                                                                                                    >  2:140436/1‑133 (MQ=255)                                                aacaaacCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCtttt                                                                                                    <  1:140436/133‑1 (MQ=255)                                                 acaaacCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCg                                                                                    >  7:388810/1‑148 (MQ=255)                                                   aaacCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTa                                                                                  <  5:85305/148‑1 (MQ=255)                                                         cACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGa                                                                                                <  6:179395/128‑1 (MQ=255)                                                         cACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGa                                                                                                >  5:179395/1‑128 (MQ=255)                                                           cgcgCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTc                                                                             >  4:151195/1‑145 (MQ=255)                                                            gcgcAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAg                                                                          >  2:158916/1‑147 (MQ=255)                                                                  gCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGc                                                                   <  3:151195/148‑1 (MQ=255)                                                                                        tATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTg                                                                    >  1:349639/1‑125 (MQ=255)                                                                                        tATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTg                                                                    <  2:349639/125‑1 (MQ=255)                                                                                        tATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGa                                             >  7:346313/1‑148 (MQ=255)                                                                                           gcgcATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTTGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACtt                                          >  2:180462/1‑148 (MQ=255)                                                                                             gcATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTc                                        <  8:388810/148‑1 (MQ=255)                                                                                              cATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCa                                       >  2:99025/1‑148 (MQ=255)                                                                                              cATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCa                                       >  6:2376/1‑148 (MQ=255)                                                                                                tCAAAAGTGGAAGCTTTCTCCGCAACCACACTCATCTTTGACGTTCGGTTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGTCTCTTGCCATCGACCACGACGTTc                                        <  1:101286/145‑1 (MQ=255)                                                                                                 cAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGaa                                                                        >  4:38169/1‑112 (MQ=255)                                                                                                 cAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGaa                                                                        <  3:38169/112‑1 (MQ=255)                                                                                                 cAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGc                                    >  4:77171/1‑148 (MQ=255)                                                                                                 cAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGc                                    >  8:224456/1‑148 (MQ=255)                                                                                                 cAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTCACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGATCTTGCCATCTACCACGACTTTCa                                       >  2:101286/1‑145 (MQ=255)                                                                                                   aaCGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCt                                  >  7:381331/1‑148 (MQ=255)                                                                                                    aCGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCtt                                 <  3:77171/148‑1 (MQ=255)                                                                                                    aCGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCta                                 >  5:357194/1‑147 (MQ=255)                                                                                                      gTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTTGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTg                               <  1:180462/148‑1 (MQ=255)                                                                                                      gTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCGTCGACCACGACTTTCACGCCTTTg                               <  6:367220/148‑1 (MQ=255)                                                                                                                tCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAaca                     >  5:193335/1‑148 (MQ=255)                                                                                                                  gccgcAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAACACg                   <  7:224456/148‑1 (MQ=255)                                                                                                                  gccgcAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAACACg                   <  5:248211/148‑1 (MQ=255)                                                                                                                  gccgcAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTATGTCTTCAAACACg                   <  6:357194/148‑1 (MQ=255)                                                                                                                   ccgcAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAACACGa                  <  8:346313/148‑1 (MQ=255)                                                                                                                            cacTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAACACGATGTCTTCCg         >  2:385250/1‑148 (MQ=255)                                                                                                                             acTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAACACGATGTCTTCCgg        >  6:370855/1‑148 (MQ=255)                                                                                                                              cTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTg                               >  3:127201/1‑124 (MQ=255)                                                                                                                              cTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTg                               <  4:127201/124‑1 (MQ=255)                                                                                                                                 aTCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAACACGATGTCTTCCGGCGTc    <  8:381331/148‑1 (MQ=255)                                                                                                                                   cTTTGAAGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACAGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAACACGATGTCTTCCGGCGTCgg  >  2:101725/1‑148 (MQ=255)                                                                                                                                     ttGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTc                                        >  4:192479/1‑108 (MQ=255)                                                                                                                                     ttGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCGCTTGCCATCGACCACGACTTTc                                        <  3:192479/108‑1 (MQ=255)                                                                                                                                            tCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAACACGATGTCTTCCGGCGTCg   <  1:347812/138‑1 (MQ=255)                                                                                                                                            tCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAACACGATGTCTTCCGGCGTCg   >  2:347812/1‑138 (MQ=255)                                                                                                                                             |                                                                                                                                          AACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTTATCGGTATGCGCATCAAACGTGGAAGCTTTCGCCGCAACCACACTCATCTTTGACGTTCGGGTTGGTGAATTTAAACCCTTCGTTCAGGCCTTCTTTTACGAAGTCCAGCTGCGTACCGTCCAGAAATTGCAGGCTCTTGCCATCGACCACGACTTTCACGCCTTTGTCTTCAAACACGATGTCTTCCGGCGTCGG  >  NC_000913/2659468‑2659745

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