Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,245,151	T→A	7.7%	N540I (AAC→ATC)	cirA ←	colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,245,151	0	T	A	7.7%	176.8 / 5.9	65	N540I (AAC→ATC)	cirA	colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter
Reads supporting (aligned to +/- strand):  ref base T (32/27);  new base A (3/2);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

ACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGT  >  NC_000913/2245008‑2245299                                                                                                                                                |                                                                                                                                                     aCCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAGCGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAAc                                                                                                                                                 <  2:1226293/149‑1 (MQ=255)  ccAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGCCCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACg                                                                                                                                                >  1:467928/1‑149 (MQ=255)     gTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGaa                                                                                                                                             <  2:556528/149‑1 (MQ=255)      tCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAAt                                                                                                                                            >  1:1209158/1‑149 (MQ=255)        ttCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAAttt                                                                                                                                          >  1:495305/1‑149 (MQ=255)        ttCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAAttt                                                                                                                                          >  2:1545989/1‑149 (MQ=255)        ttCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAAttt                                                                                                                                          >  2:429493/1‑149 (MQ=255)          ccagcgccagcgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCGTTCAACGATGAATGGaaactgt                                                                                                                                        <  1:926539/149‑6 (MQ=255)               gccagcgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTc                                                                                                                                   <  2:549699/149‑1 (MQ=255)                ccagcgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCAGTCATCGGTGAACGGAATTTTCAGTTCg                                                                                                                                  >  2:1260651/1‑149 (MQ=255)                  agcgGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGATTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGt                                                                                                                                <  1:1545522/149‑1 (MQ=255)                       tttCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCa                                                                                                                           >  1:1406517/1‑149 (MQ=255)                         tccagtccagCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACg                                                                                                                         >  1:1495187/1‑149 (MQ=255)                           cagtccagCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGcc                                                                                                                       >  1:12663/1‑149 (MQ=255)                                 agCGTTCCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCAACGTTGAACGGAATTTTCAGTTCGTTTTCCACGCCCTGatt                                                                                                                 >  1:515615/1‑147 (MQ=255)                                 agCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAAt                                                                                                                 <  1:185086/149‑1 (MQ=255)                                     tACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGCACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGa                                                                                                             >  2:1165033/1‑149 (MQ=255)                                     tACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGa                                                                                                             <  2:199113/149‑1 (MQ=255)                                          ttAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCttt                                                                                                        <  1:672097/149‑1 (MQ=255)                                              cAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTa                                                                                                    <  1:849516/149‑1 (MQ=255)                                                ctATGGAACGGCAGATCGGCTAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTTGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAAc                                                                                                  <  2:139759/148‑1 (MQ=255)                                                   tGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGtt                                                                                               <  1:664613/149‑1 (MQ=255)                                                    gtaaCGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                              <  2:129625/147‑1 (MQ=255)                                                    ggAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                              <  2:1495187/149‑1 (MQ=255)                                                    ggAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                              <  2:342204/149‑1 (MQ=255)                                                    ggAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTg                                                                                              <  2:1209158/149‑1 (MQ=255)                                                     gAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTgt                                                                                             >  1:770616/1‑149 (MQ=255)                                                              aTCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTaa                                                                                    >  2:1221424/1‑149 (MQ=255)                                                              aTCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTaa                                                                                    <  2:1406517/149‑1 (MQ=255)                                                              aTCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCAAGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTaa                                                                                    >  2:1109823/1‑149 (MQ=255)                                                                 ggATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATa                                                                                 >  2:1505750/1‑149 (MQ=255)                                                                   aTAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATAcc                                                                               <  2:1287956/149‑1 (MQ=255)                                                                    tAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCg                                                                              <  2:578960/149‑1 (MQ=255)                                                                     aGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGTTAAATACCgg                                                                             <  2:1358562/149‑1 (MQ=255)                                                                     aGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCgg                                                                             >  1:113690/1‑149 (MQ=255)                                                                              ttttCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGccg                                                                    >  1:491358/1‑149 (MQ=255)                                                                              ttttCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGccg                                                                    >  1:1521897/1‑149 (MQ=255)                                                                                    ccGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGtt                                                              >  2:588849/1‑149 (MQ=255)                                                                                          cTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGcc                                                        >  1:187826/1‑149 (MQ=255)                                                                                               aTCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGtct                                                   >  1:516023/1‑149 (MQ=255)                                                                                               aTCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGtct                                                   >  1:1512044/1‑149 (MQ=255)                                                                                                tCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGtctc                                                  >  1:346821/1‑149 (MQ=255)                                                                                                    gACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCaaaa                                              <  1:1191567/149‑1 (MQ=255)                                                                                                     aCCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAAc                                             <  2:541195/149‑1 (MQ=255)                                                                                                       cATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACca                                           >  2:1441602/1‑149 (MQ=255)                                                                                                         tCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACcaac                                         >  2:1121695/1‑149 (MQ=255)                                                                                                         tCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACcaac                                         >  1:1066842/1‑149 (MQ=255)                                                                                                               tAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGtt                                   >  1:815806/1‑149 (MQ=255)                                                                                                                aGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTGCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGttt                                  <  1:976589/149‑1 (MQ=255)                                                                                                                 ggTGTAGTTGATCGACAGTTTCCATTCATCGTTGAAAGGAATTTTCAGTTCGGTTTCCACGCGCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGTTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGtttt                                 <  2:467928/149‑1 (MQ=255)                                                                                                                            tCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGt                      <  1:429493/149‑1 (MQ=255)                                                                                                                              tgaaaattCCGTTCAACGATGAATGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGc                    >  2:733893/7‑149 (MQ=255)                                                                                                                              gACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGc                    >  2:1387369/1‑149 (MQ=255)                                                                                                                               gaaaattCCGTTCAACAATGAATGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCa                   >  1:227533/6‑149 (MQ=255)                                                                                                                               aCAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCa                   >  1:1557174/1‑149 (MQ=255)                                                                                                                               aCAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCa                   >  1:1470429/1‑149 (MQ=255)                                                                                                                                 aGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCACCAAAGTTTTGGTAGCCCGGTGCAGc                 <  2:770616/149‑1 (MQ=255)                                                                                                                                  aattCCGTTCAACGATGAATGGAATTTTCAGTTCGGTTTCCACGCCCTGTATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTgg                               >  2:123491/3‑134 (MQ=255)                                                                                                                                  aattCCGTTCAACGATGAATGGAATTTTCAGTTCGGTTTCCACGCCCTGTATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTgg                               <  1:123491/132‑1 (MQ=255)                                                                                                                                    ttCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGtt              <  2:1368960/149‑1 (MQ=255)                                                                                                                                    ttCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGtt              <  2:1512044/149‑1 (MQ=255)                                                                                                                                      ccATTCATCGTTGAACGGAATTTTCTGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGGGCAGCGTTGa            >  2:889509/1‑149 (MQ=255)                                                                                                                                         ttCATCGTTGTACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGt         <  1:1545989/149‑1 (MQ=255)                                                                                                                                            aTCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAg      >  1:707477/1‑149 (MQ=255)                                                                                                                                               gTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACg   >  2:621317/1‑149 (MQ=255)                                                                                                                                                ttGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGt  >  1:654924/1‑149 (MQ=255)                                                                                                                                                ttGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGt  >  2:753961/1‑149 (MQ=255)                                                                                                                                                ttGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGt  >  2:879014/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     ACCAGTCTTCCAGCGCCAGCGGTTTCCAGTCCAGCGTACCGTTAGCAGTATGGAACGGCAGATCGGATAGCGGTTTGTTTTCGCCGTTGCTGACATCACGACCATCGTTGTAGGTGTAGTTGATCGACAGTTTCCATTCATCGTTGAACGGAATTTTCAGTTCGGTTTCCACGCCCTGAATACGAGCTTTGTTAACGTTGTAGTAGCTAAATACCGGTATGCGCCGTCCGTTAGCGCCCGTCTCAAAACCAACAAAGTTTTGGTAGCCCGGTGCAGCGTTGACGTCAGACGT  >  NC_000913/2245008‑2245299

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