Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,543,230	C→G	5.8%	M277I (ATG→ATC) ‡	narU ←	nitrate/nitrite transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,543,230	0	C	G	5.8%	240.6 / 5.5	69	M277I (ATG→ATC) ‡	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base G (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATG  >  NC_000913/1543083‑1543378                                                                                                                                                     |                                                                                                                                                     cATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCa                                                                                                                                                     <  2:655083/149‑1 (MQ=255)  aTAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCAt                                                                                                                                                    >  2:272949/1‑149 (MQ=255)   tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATg                                                                                                                                                   >  2:340890/1‑149 (MQ=255)   tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATg                                                                                                                                                   <  3:607460/149‑1 (MQ=255)         tAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCaaaa                                                                                                                                             >  4:268202/1‑149 (MQ=255)          aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGC‑TGGCaaaaa                                                                                                                                            <  1:128206/149‑2 (MQ=255)             gTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATgg                                                                                                                                                  >  4:546621/1‑140 (MQ=255)             gTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATgg                                                                                                                                                  <  3:546621/140‑1 (MQ=255)                gATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAg                                                                                                                                      >  4:622709/1‑149 (MQ=255)                 aTCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGa                                                                                                                                     >  2:168300/1‑149 (MQ=255)                       gTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAAcc                                                                                                                               >  2:102167/1‑149 (MQ=255)                            ccGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGc                                                                                                                                        <  1:51243/135‑1 (MQ=255)                            ccGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGc                                                                                                                                        >  2:51243/1‑135 (MQ=255)                            ccGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATaa                                                                                                                          >  3:630842/1‑149 (MQ=255)                              acACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGAACGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGTAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAAc                                                                                                                        >  4:225135/2‑149 (MQ=255)                                      ccgAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGaa                                                                                                                >  1:100360/1‑149 (MQ=255)                                                  gAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATaaa                                                                                                                         >  1:176621/1‑128 (MQ=255)                                                  gAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATaaa                                                                                                                         <  2:176621/128‑1 (MQ=255)                                                  gAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                    >  1:33465/1‑149 (MQ=255)                                                     aTAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAg                                                                                                 >  1:95588/1‑149 (MQ=255)                                                      tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGc                                                                                                <  3:115344/149‑1 (MQ=255)                                                        gCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                              <  3:135499/149‑1 (MQ=255)                                                         caccacCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAgg                                                                                             <  1:144753/149‑1 (MQ=255)                                                           ccaccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCt                                                                                           >  2:155807/1‑149 (MQ=255)                                                            caccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTc                                                                                          >  4:641357/1‑149 (MQ=255)                                                             accaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGG‑TTTTGCCA‑TGCTGGCaaaaac                                                                                                                                           <  1:168158/98‑3 (MQ=255)                                                             accaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGG‑TTTTGCCA‑TGCTGGCaaaaac                                                                                                                                           >  2:168158/1‑96 (MQ=255)                                                                   cgaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGcc                                                                                   <  1:302649/149‑1 (MQ=255)                                                                    gaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTaaa                                                                                                  <  4:468879/133‑1 (MQ=255)                                                                    gaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTaaa                                                                                                  >  3:468879/1‑133 (MQ=255)                                                                    gaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCa                                                                                  >  2:469980/1‑149 (MQ=255)                                                                            gATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAgaga                                                                              >  2:355465/1‑145 (MQ=255)                                                                            gATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAgaga                                                                              <  1:355465/145‑1 (MQ=255)                                                                              tGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAgg                                                                        <  1:168300/149‑1 (MQ=255)                                                                                  aCCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGtt                                                                    <  2:100360/149‑1 (MQ=255)                                                                                        aaaTGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGC‑TGGCaaaaa                                                                                                                                            <  3:462678/71‑2 (MQ=38)                                                                                        aaaTGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGC‑TGGCaaaaa                                                                                                                                            >  4:462678/1‑70 (MQ=38)                                                                                          aTGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACa                                                            <  2:33465/149‑1 (MQ=255)                                                                                            gggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                          <  3:268202/149‑1 (MQ=255)                                                                                            gggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                          <  2:95588/149‑1 (MQ=255)                                                                                              gCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGa                                                        >  2:224121/1‑149 (MQ=255)                                                                                               ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCACAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAg                                                       >  3:406502/1‑149 (MQ=255)                                                                                               ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGt                                                                  >  2:68747/1‑138 (MQ=255)                                                                                               ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGt                                                                  <  1:68747/138‑1 (MQ=255)                                                                                                cAAAGAACGCCAGGCGCAGAATATTGACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGGAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                >  3:198802/1‑139 (MQ=255)                                                                                                cAAAGAACGCCAGGCGCAGAATATTGACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGGAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                                <  4:198802/139‑1 (MQ=255)                                                                                                 aaaGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCt                                                     >  1:5460/1‑149 (MQ=255)                                                                                                  aaGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGAc                                                             >  4:585675/1‑140 (MQ=255)                                                                                                  aaGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGAc                                                             <  3:585675/140‑1 (MQ=255)                                                                                                            ggCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAAtt                                          <  1:224121/149‑1 (MQ=255)                                                                                                                 aGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgc                                     >  3:628826/1‑149 (MQ=255)                                                                                                                   aaTATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgcgc                                   <  1:340890/149‑1 (MQ=255)                                                                                                                   aaTATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgcgc                                   <  1:155807/149‑1 (MQ=255)                                                                                                                     tatTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGt                                 <  1:469980/149‑1 (MQ=255)                                                                                                                     tatTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGt                                 >  4:268901/1‑149 (MQ=255)                                                                                                                          acaTCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACt                            >  3:352223/1‑149 (MQ=255)                                                                                                                             tCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGc                         <  2:5460/149‑1 (MQ=255)                                                                                                                               cGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGa                       <  1:102167/149‑1 (MQ=255)                                                                                                                               cGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGa                       <  3:622709/149‑1 (MQ=255)                                                                                                                                  gAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGatat                    <  4:628826/149‑1 (MQ=255)                                                                                                                                      tGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcatt                >  1:256867/1‑149 (MQ=255)                                                                                                                                      tGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcatt                >  1:552802/1‑149 (MQ=255)                                                                                                                                      tGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcat                 >  2:283965/1‑148 (MQ=255)                                                                                                                                      tGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcat                 <  1:283965/148‑1 (MQ=255)                                                                                                                                        gggttttGCCA‑TGCTGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGtt                                                                    <  4:122110/92‑1 (MQ=255)                                                                                                                                        gggttttGCCA‑TGCTGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGtt                                                                    >  3:122110/4‑95 (MQ=255)                                                                                                                                         gTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat             <  4:630842/149‑1 (MQ=255)                                                                                                                                          tttttGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc            <  1:272949/149‑1 (MQ=255)                                                                                                                                                 cA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAcc     >  2:82549/1‑149 (MQ=255)                                                                                                                                                     cATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATg  <  3:641357/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATG  >  NC_000913/1543083‑1543378

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