Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,850,155	T→G	5.2%	Q94P (CAG→CCG)	hycA ←	regulator of the transcriptional regulator FhlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,850,155	0	T	G	5.2%	276.7 / 2.7	77	Q94P (CAG→CCG)	hycA	regulator of the transcriptional regulator FhlA
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base G (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCGGCGCGCAGCTCATGG  >  NC_000913/2850008‑2850302                                                                                                                                                    |                                                                                                                                                    cTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTg                                                                                                                                                    >  4:451154/1‑149 (MQ=255)  ttAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCttc                                                                                                                                                   >  1:426729/1‑147 (MQ=255)  ttAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTgc                                                                                                                                                   >  2:286174/1‑149 (MQ=255)       ttCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCGGTTGCGCAGGCGCagc                                                                                                                                              <  4:126478/149‑3 (MQ=255)       ttCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAAt                                                                                                                                              <  4:123694/149‑1 (MQ=255)        tCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATc                                                                                                                                             >  4:30997/1‑149 (MQ=255)            cGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCt                                                                                                                                         >  1:294965/1‑149 (MQ=255)              gCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCt                                                                                                                                                     <  3:354385/135‑1 (MQ=255)              gCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCt                                                                                                                                                     >  4:354385/1‑135 (MQ=255)                acaTTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGttt                                                                                                                                     <  1:140130/149‑1 (MQ=255)                acaTTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGttt                                                                                                                                     >  4:431703/1‑149 (MQ=255)                 caTTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATcc                                                                                                                                            <  1:25452/141‑1 (MQ=255)                 caTTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATcc                                                                                                                                            >  2:25452/1‑141 (MQ=255)                      ttttCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTc                                                                                                                               >  4:38515/1‑149 (MQ=255)                           aaTAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCAt                                                                                                                          >  1:366824/1‑149 (MQ=255)                                cGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCtttt                                                                                                                     <  4:367904/149‑1 (MQ=255)                                     tataAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTc                                                                                                                >  2:267732/1‑149 (MQ=255)                                        tcgCTGGCGTTTACGACAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTc                                                                                                                >  1:411079/3‑146 (MQ=255)                                        aaGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGAc                                                                                                             >  3:343022/1‑149 (MQ=255)                                              tCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATa                                                                                                       <  2:444769/149‑1 (MQ=255)                                                 tAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTc                                                                                                    >  4:315974/1‑149 (MQ=255)                                                   aaCGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGa                                                                                                  <  4:460039/149‑1 (MQ=255)                                                     cGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCAGGCGCaactga                                                                                                                                           >  2:142110/1‑102 (MQ=255)                                                     cGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCAGGCGCaactga                                                                                                                                           <  1:142110/106‑5 (MQ=255)                                                       ccAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGgt                                                                                              <  3:256445/149‑1 (MQ=255)                                                        cAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGgtg                                                                                             >  2:205022/1‑149 (MQ=255)                                                         aGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGgtgt                                                                                            <  4:146708/149‑1 (MQ=255)                                                         aGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGgtgt                                                                                            >  2:74549/1‑149 (MQ=255)                                                             aTTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTgcgc                                                                                                                                                 >  3:322277/1‑92 (MQ=255)                                                             aTTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTgcgc                                                                                                                                                 <  4:322277/92‑1 (MQ=255)                                                             aTTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTTCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGTTGTGGCt                                                                                        <  2:426729/149‑1 (MQ=255)                                                             aTTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCt                                                                                        <  3:376937/149‑1 (MQ=255)                                                             aTTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCt                                                                                        >  2:166692/1‑149 (MQ=255)                                                                  ttCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTaa                                                                                   <  3:119508/149‑1 (MQ=255)                                                                    cGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAg                                                                                 <  4:116798/149‑1 (MQ=255)                                                                        aTAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCAt                                                                             >  2:431194/1‑149 (MQ=255)                                                                            tgtTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGc                                                                         >  3:97156/1‑149 (MQ=255)                                                                                   gAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTg                                                                  <  1:337580/149‑1 (MQ=255)                                                                                     aCCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGAc                                                                <  3:94906/149‑1 (MQ=255)                                                                                     aCCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGAc                                                                >  4:302430/1‑149 (MQ=255)                                                                                       cTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGt                                                              >  1:94225/1‑149 (MQ=255)                                                                                        tGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTa                                                             >  4:107163/1‑149 (MQ=255)                                                                                        tGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTa                                                             >  3:471940/1‑149 (MQ=255)                                                                                          ttCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTaaa                                                           >  1:97016/1‑149 (MQ=255)                                                                                          ttCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTaaa                                                           <  2:386934/149‑1 (MQ=255)                                                                                          ttCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTaaa                                                           >  1:203074/1‑149 (MQ=255)                                                                                               gATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGa                                                      >  4:25265/1‑149 (MQ=255)                                                                                                aTCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGaa                                                     >  4:260541/1‑149 (MQ=255)                                                                                                            gACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCaaaaa                                         <  3:431703/149‑1 (MQ=255)                                                                                                            gACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCaaaaa                                         <  2:97016/149‑1 (MQ=255)                                                                                                             aCCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAg                                        <  1:74549/149‑1 (MQ=255)                                                                                                             aCCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAg                                        <  1:267732/149‑1 (MQ=255)                                                                                                             aCCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAg                                        >  1:322532/1‑149 (MQ=255)                                                                                                             aCCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTATGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAg                                        <  4:343022/149‑1 (MQ=255)                                                                                                              ccTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAgg                                       <  2:322532/149‑1 (MQ=255)                                                                                                               cTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGa                                      >  2:158126/1‑149 (MQ=255)                                                                                                                ttGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGAc                                     >  3:274225/1‑149 (MQ=255)                                                                                                                     atcaatcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTc                                              <  3:10897/135‑1 (MQ=255)                                                                                                                     atcaatcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTc                                              >  4:10897/1‑135 (MQ=255)                                                                                                                     atcaatcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAAc                                >  3:236848/1‑149 (MQ=255)                                                                                                                     atcaatcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAAc                                >  4:173782/1‑149 (MQ=255)                                                                                                                     atcaatcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAAc                                >  3:18454/1‑149 (MQ=255)                                                                                                                          tcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGt                           <  3:260541/149‑1 (MQ=255)                                                                                                                          tcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCAGTTTGTCTCCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGt                           <  3:451154/149‑1 (MQ=255)                                                                                                                            aTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATc                                                                            <  1:58098/98‑1 (MQ=255)                                                                                                                            aTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATc                                                                            >  2:58098/1‑98 (MQ=255)                                                                                                                             tGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCt                        <  3:38515/149‑1 (MQ=255)                                                                                                                             tGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCt                        <  3:25265/149‑1 (MQ=255)                                                                                                                              gCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCtt                       <  2:294965/149‑1 (MQ=255)                                                                                                                                 gTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGt                    >  4:466777/1‑149 (MQ=255)                                                                                                                                 gTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATACGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGt                    <  1:286174/149‑1 (MQ=255)                                                                                                                                  tCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTc                   <  3:107163/149‑1 (MQ=255)                                                                                                                                    aaTGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCg                 <  3:466777/149‑1 (MQ=255)                                                                                                                                       gCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCGgcg              >  1:115759/1‑149 (MQ=255)                                                                                                                                           gattGCGCAGGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCGGCGCGCa          >  3:460301/3‑149 (MQ=255)                                                                                                                                               gcgcCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCGGCGCGCAGCTc      <  2:94225/149‑1 (MQ=255)                                                                                                                                                 gcCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGc                                                                                >  2:208359/1‑73 (MQ=255)                                                                                                                                                 gcCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGc                                                                                <  1:208359/73‑1 (MQ=255)                                                                                                                                                  ccTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCGGCGCGCAGCTCATg   <  4:97156/149‑1 (MQ=255)                                                                                                                                                   cTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCGGCGCGCAGCTCATgg  >  4:13810/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCGGCGCGCAGCTCATGG  >  NC_000913/2850008‑2850302

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