Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,543,229	G→T	5.8%	L278M (CTG→ATG)	narU ←	nitrate/nitrite transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,543,229	0	G	T	5.8%	244.7 / 4.5	69	L278M (CTG→ATG)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base T (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‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACC  >  NC_000913/1543083‑1543375                                                                                                                                                    |                                                                                                                                                   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                                                                                                                                               <  3:546621/140‑1 (MQ=255)             gTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATgg                                                                                                                                               >  4:546621/1‑140 (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                                                                                                                                     >  2:51243/1‑135 (MQ=255)                            ccGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGc                                                                                                                                     <  1:51243/135‑1 (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                                                                                                                                        >  2:168158/1‑96 (MQ=255)                                                             accaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGG‑TTTTGCCA‑TGCTGGCaaaaac                                                                                                                                        <  1:168158/98‑3 (MQ=255)                                                                   cgaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGcc                                                                                <  1:302649/149‑1 (MQ=255)                                                                    gaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTaaa                                                                                               >  3:468879/1‑133 (MQ=255)                                                                    gaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTaaa                                                                                               <  4:468879/133‑1 (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                                                                                                                                         >  4:462678/1‑70 (MQ=38)                                                                                        aaaTGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGC‑TGGCaaaaa                                                                                                                                         <  3:462678/71‑2 (MQ=38)                                                                                          aTGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACa                                                         <  2:33465/149‑1 (MQ=255)                                                                                            gggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                       <  2:95588/149‑1 (MQ=255)                                                                                            gggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                       <  3:268202/149‑1 (MQ=255)                                                                                              gCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGa                                                     >  2:224121/1‑149 (MQ=255)                                                                                               ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCACAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAg                                                    >  3:406502/1‑149 (MQ=255)                                                                                               ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGt                                                               <  1:68747/138‑1 (MQ=255)                                                                                               ccAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGt                                                               >  2:68747/1‑138 (MQ=255)                                                                                                cAAAGAACGCCAGGCGCAGAATATTGACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGGAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                             <  4:198802/139‑1 (MQ=255)                                                                                                cAAAGAACGCCAGGCGCAGAATATTGACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGGAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                             >  3:198802/1‑139 (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:155807/149‑1 (MQ=255)                                                                                                                   aaTATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgcgc                                <  1:340890/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              <  1:283965/148‑1 (MQ=255)                                                                                                                                      tGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcat              >  2:283965/1‑148 (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)                                                                                                                                                    |                                                                                                                                                   CATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCA‑GCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACC  >  NC_000913/1543083‑1543375

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