Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,543,230	C→G	11.3%	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	11.3%	167.6 / 14.9	62	M277I (ATG→ATC) ‡	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (31/24);  new base G (4/3);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.14e-01

AAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGC  >  NC_000913/1543092‑1543379                                                                                                                                           |                                                                                                                                                       aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGC‑aaaaac                                                                                                                                            <  1:490750/150‑3 (MQ=255) aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGAGATGGCACCGATAAATGGGCAAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGC‑aaaaac                                                                                                                                            <  2:823939/150‑3 (MQ=255)  aaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGC‑AAAAacc                                                                                                                                           <  1:1082674/150‑1 (MQ=255)             cGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAAcc                                                                                                                                <  1:781309/150‑1 (MQ=255)             cGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGCGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAAcc                                                                                                                                <  1:848774/150‑1 (MQ=255)                  cccGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCCCCGATAAATGGGCCAAAGACCGCCAGGCGCAGAATAGTCCCATCCGGGAACTGGGTTTTTGCCAGCATGGC‑‑ACCCCCGCACAAAAACCGATaaa                                                                                                                          >  1:134841/1‑150 (MQ=255)                     gCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACg                                                                                                                        <  1:1314671/150‑1 (MQ=255)                           cgccgAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑ACAACCCGCAGAAAAACCGATAACCGAGCCGa                                                                                                                  >  1:365085/1‑150 (MQ=255)                            gccgAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGaa                                                                                                                 >  1:142352/1‑150 (MQ=255)                                  cTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGc                                                                                                           >  2:347327/1‑150 (MQ=255)                                       cGGAATTAGCACCCCCAACCGACCGCGCGATGGCACCGATTAGTGGGCCATTGAACGGCAGGCCCAGACTCTTCACATCCGGGAACTGGGTTTTTCCCTGCAGGGC‑AATCCCCTCATAAAATCCGATACACGAACCGAAGGAGGCAAggg                                                                                                      >  2:851621/1‑149 (MQ=255)                                       cGGAAATAGCACCACCAACCGACCGCGCTAAGGCAACGATAAATTGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGGTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGt                                                                                                      <  1:243452/150‑1 (MQ=255)                                        ggAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGACCGCCAGGCGCAGAATATTCACATCCGGGATCTGGGTTTTTGCCAGCATGGC‑AAAAGCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                     >  1:83854/1‑150 (MQ=255)                                        ggAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGtt                                                                                                     >  2:342691/1‑149 (MQ=255)                                          aaaTAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGAAAAACAAACCGAAGGTGGCAAGGTaaa                                                                                                   <  2:393903/150‑1 (MQ=255)                                             tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATACACGAACCGAAGGTGGCAAGGTAAAGcc                                                                                                >  1:1146159/1‑149 (MQ=255)                                             tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCa                                                                                                >  1:717194/1‑150 (MQ=255)                                              aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AACACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                               >  1:1267119/1‑150 (MQ=255)                                              aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                               <  2:690124/150‑1 (MQ=255)                                              aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                               <  1:1299836/150‑1 (MQ=255)                                                    accaACCGACAGCGCGATGGCACCGATAAATGGGCCAAACACCGCCCGGCGCAGAATAGTCACATCCGGGAACTGGGTTTATGCCAGCATGGA‑AAAACACGCAGAAAAACCGAAAAACGCACCCAAGGTGGCCAGGTCAAGCAGGcgccg                                                                                         >  2:810583/1‑146 (MQ=255)                                                      caACCGACCGCGCGATGGCACCGATAAATGGGCCACAGAACGCCAGGCGCAGAATATTCACATCCGGGAACGGGCTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAATGTGGCAAGGTAAAGCAGGCTcagca                                                                                       <  1:794201/150‑1 (MQ=255)                                                      caACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGACCTGGGTTTATGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTccgcc                                                                                       >  2:811115/1‑149 (MQ=255)                                                             ccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCagag                                                                                <  1:904631/150‑1 (MQ=255)                                                                  cgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCa                                                                           >  2:548894/1‑150 (MQ=255)                                                                   gATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAg                                                                          >  2:1036242/1‑150 (MQ=255)                                                                    aTGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAgg                                                                         >  1:513020/1‑150 (MQ=255)                                                                        cACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGAAAGGTAAAGCAGGATCAGAAGCCAGAGATGCAGGCGtt                                                                     >  2:477482/1‑150 (MQ=255)                                                                                 aTGGGCCAAAGACCGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGA‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAg                                                            <  2:1135611/150‑1 (MQ=255)                                                                                 aTGGGCCAAAGAACGCCAGTCGCAGATTATTCACATCCGGGAACTGGGTTTTTGCCTGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAg                                                            <  1:995087/150‑1 (MQ=255)                                                                                  tGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                           <  2:1003134/150‑1 (MQ=255)                                                                                    ggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGa                                                         >  2:1246214/1‑150 (MQ=255)                                                                                       cAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGGTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCt                                                      <  2:1076878/150‑1 (MQ=255)                                                                                           gAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTc                                                  >  1:1091023/1‑150 (MQ=255)                                                                                             aCGCCAGGCGCAGAATATTCACATCCGGGAACTGTGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCgg                                                <  1:848695/150‑1 (MQ=255)                                                                                                ccAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCCGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCaa                                             <  2:807116/150‑1 (MQ=255)                                                                                                 cAGGCGCAGAATATTCACATCCTGGATCTTGGTTTTTGCCAGCATTGC‑AAAACCCGCAGAAAAACCTATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCATCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAAt                                            <  1:1003710/150‑1 (MQ=255)                                                                                                 cAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACGGGGAGGTGGTCGGCatt                                            >  2:47720/1‑148 (MQ=255)                                                                                                              ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                               >  1:991960/1‑150 (MQ=255)                                                                                                              ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                               >  1:986042/1‑150 (MQ=255)                                                                                                              ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                               >  2:174882/1‑150 (MQ=255)                                                                                                               tCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAAc                              <  1:365075/150‑1 (MQ=255)                                                                                                               tCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAAc                              <  1:347327/150‑1 (MQ=255)                                                                                                                cacaTCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACt                             >  1:365902/1‑150 (MQ=255)                                                                                                                   ctCCGGGAACTGGGTTTTTGCCAGCCTGGC‑AAACCCCGCAGAAAAACCGATAACCGAACCGAAGGGGGCAAGGGCAAGCAGGCACAGCAGCCAGAGATGCAGGCGTTGTAAGACAGCGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGc                          >  2:930396/2‑150 (MQ=255)                                                                                                                   aTCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGc                          <  1:586029/150‑1 (MQ=255)                                                                                                                      cGGCAACGGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAACCGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGat                       <  1:954228/150‑1 (MQ=255)                                                                                                                       gggAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGata                      <  1:1193973/150‑1 (MQ=255)                                                                                                                            tgcgggttttGCCATGCTGGCAAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcat                  >  2:621510/7‑150 (MQ=255)                                                                                                                             tGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattc                <  1:477503/150‑1 (MQ=255)                                                                                                                               gggttttGCCATGCTGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat              >  2:1118953/4‑150 (MQ=255)                                                                                                                               ggTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat              >  1:671490/1‑150 (MQ=255)                                                                                                                                ggttttGCCATGCTGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc             >  1:636112/3‑150 (MQ=255)                                                                                                                                 tttttGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGGGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATcc            <  2:642606/150‑1 (MQ=255)                                                                                                                                 gttttGCCATGCTGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATcc            >  2:1116637/2‑150 (MQ=255)                                                                                                                                      gCCAGCATGGC‑AAAACCCGCAGAACAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAc       >  1:749867/1‑150 (MQ=255)                                                                                                                                       ccAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGCACCGAAGGGGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGAGTTGTATGACAGGGAGCTGGGCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAca      >  2:1218421/1‑149 (MQ=255)                                                                                                                                       ccAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCCGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGTTATCATTCCTCCCTGaac      >  2:200322/1‑148 (MQ=255)                                                                                                                                       ccAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTATGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAcc      >  2:1213583/1‑150 (MQ=255)                                                                                                                                       ccAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAcc      <  1:716108/150‑1 (MQ=255)                                                                                                                                        cAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAACGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTCGTCGGCAATTGCGGCGCGTGAACTGGCGAAATCATTCATCCCAGACCa     >  1:217151/1‑150 (MQ=255)                                                                                                                                          gCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATg   >  2:808100/1‑150 (MQ=255)                                                                                                                                           cATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGCGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATgc  >  1:362345/1‑150 (MQ=255)                                                                                                                                           |                                                                                                                                                       AAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGC  >  NC_000913/1543092‑1543379

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