Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,286,965	A→T	18.3%	intergenic (+439/+101)	narI → / ← rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,286,965	0	A	T	18.3%	140.2 / 2.5	61	intergenic (+439/+101)	narI/rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (14/34);  new base T (2/8);  total (18/43)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCG  >  NC_000913/1286818‑1287111                                                                                                                                                    |                                                                                                                                                   ccgccgtCCCCGCTTGGGCAGAATAGTTGAGTGGGGAGTCGTTTGAGCAGTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGAGTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGACg                                                                                                                                                  <  2:1181868/144‑1 (MQ=14) cATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACg                                                                                                                                                  <  2:886616/150‑1 (MQ=255) cATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACg                                                                                                                                                  >  2:1212940/1‑150 (MQ=255) cATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGAGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACg                                                                                                                                                  <  1:581875/150‑1 (MQ=255)  aTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcgc                                                                                                                                                 >  2:1377151/1‑150 (MQ=255)     cTTCCCCGCTTTGGCAGAATTTTTGATTGGGTCTTCGTTTGAGAATTCCGGGTCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCg                                                                                                                                              <  2:626540/150‑1 (MQ=14)      ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                             <  1:952649/150‑1 (MQ=255)      gtCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                             <  1:1444698/149‑1 (MQ=255)        ccccGCTTGGGCAGATTATTTGATTTCGGCTTCGTTTGAGTATTCGGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                           <  2:260784/150‑1 (MQ=14)        ccccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                           <  1:881278/150‑1 (MQ=255)        ccccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                           <  1:1147759/150‑1 (MQ=255)        ccccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                           <  2:1308833/150‑1 (MQ=255)        ccccGCGTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                           <  2:1109047/150‑1 (MQ=255)         cccGCTTGGGCAGATTTTTTGATTGCGGATTCGTTTGAGTATTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                                          <  2:361025/150‑1 (MQ=255)           cGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                                        <  1:640708/150‑1 (MQ=255)            gCTGGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGaa                                                                                                                                       <  1:864546/150‑1 (MQ=255)                           ttGATTGCGGATTCGTTTGAGAAATCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCAGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCACACTCGCAGACGCACCAGCCGCACCATAGTCGAAGcgt                                                                                                                        >  1:864625/1‑148 (MQ=14)                               ttGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                                                                                                    >  2:140754/1‑150 (MQ=255)                               ttGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                                                                                                    <  1:799545/150‑1 (MQ=255)                                              agaATTACGGGGCTTTTGAAAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCATCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                     <  1:1564836/150‑1 (MQ=255)                                               gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                    <  1:535620/150‑1 (MQ=255)                                                aaTTCCGGGGCTTTTGACAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                                   <  2:1072665/150‑1 (MQ=255)                                                 aTTCCTGGGCTTTTGATAGTGAGGGTTTTGGGGGTAGGATTACTCCGCGCTGCGCGCTTCTCGCTACGGGTCTTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                  <  2:1455638/150‑1 (MQ=14)                                                 aTTCCGGGGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTCCCCCGCATGGGCa                                                                                                  <  2:1385373/150‑1 (MQ=255)                                                  ttCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAg                                                                                                 >  1:990230/1‑150 (MQ=255)                                                     cGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAAAGCTCACTCGCTGTCGCTCGAGGCGAACCATACTCCAAGCGTCTCATCCATCCCCGCCTGGGccgcaa                                                                                              >  2:352601/1‑145 (MQ=11)                                                     cGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACACATCGCACCGCGCTCCGCCCTTCGCCGCGTGTACTGCGGCCACTCTCTCTCGCTGGCCCTCGAGGCCAAAAATAGTCGCGGCGTCTCGTCCTTCCCCGCGGGGGccggct                                                                                              >  2:1417017/1‑145 (MQ=2)                                                           ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCGGACGCACGAGCCCAACCATAGTCGAAGCATATCATCCTACCCCGCAAGGGCAGCATATTTgt                                                                                        >  1:1094883/1‑149 (MQ=11)                                                             ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                                      <  2:776432/150‑1 (MQ=39)                                                                aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGGCGCGCGAGCCGAACCATAGTCGAAGCTTCTCATACAACCCCCCATGGGCCGAAGATATGATTGCg                                                                                   >  1:1407282/1‑150 (MQ=11)                                                                       tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGc                                                                            <  1:866001/150‑1 (MQ=40)                                                                         gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                          <  2:573440/150‑1 (MQ=255)                                                                         gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                          <  1:1501590/150‑1 (MQ=255)                                                                         gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                          <  2:232102/150‑1 (MQ=255)                                                                                   aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCAGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAACCTTACCCGCATGGGCAGAATAATTAATTGCGGATTCGGTGGGGAGGTCAg                                                                >  2:589379/1‑150 (MQ=11)                                                                                   aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAAAGCTATATCGCTTACGCTCGAGCCGAACCTTAGTAGAAGATTCTCATCAATCCCCGCATTGGCAGAATATTTTATTTCGGATTCGAGTGGTAGTTCAg                                                                >  2:589426/1‑150 (MQ=11)                                                                                      gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                                             <  2:990230/150‑1 (MQ=255)                                                                                               gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                    <  2:879670/150‑1 (MQ=255)                                                                                               gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                    <  1:946337/150‑1 (MQ=255)                                                                                                   tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                <  1:746728/150‑1 (MQ=255)                                                                                                   tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTGGTGTAAGTgt                                                >  2:871012/1‑149 (MQ=255)                                                                                                     cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtg                                              <  1:140754/150‑1 (MQ=255)                                                                                                          ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                         <  2:342566/150‑1 (MQ=255)                                                                                                          ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                                         <  1:1676844/150‑1 (MQ=255)                                                                                                            cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                       <  1:1265050/150‑1 (MQ=255)                                                                                                            cGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                       <  1:1265062/150‑1 (MQ=255)                                                                                                             gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                      <  1:1212940/150‑1 (MQ=255)                                                                                                             gCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGACGCTCGAGCCGAACCATAGTCGAATCTTCACATACATCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGCGAGGTCAGGGACATTTGAAAGTGATGGTGGTggg                                      >  2:620669/1‑150 (MQ=14)                                                                                                              cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAGGGGCAGAAGATGTGATTGCGGATTCGCTTGAGAGTTCCGGGACTTTTGAAAGTGATTGTGgggggg                                     >  2:853560/1‑150 (MQ=14)                                                                                                                 ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                  <  1:451045/150‑1 (MQ=255)                                                                                                                  tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                 <  1:1016856/150‑1 (MQ=255)                                                                                                                  tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                 <  1:1116406/150‑1 (MQ=255)                                                                                                                         gTTGCCTGCGGCAACGCTCTATCGCAGTCGCTCGAGCCGAACCATAGTCGAAGCATCTCATCCTTCCCCGCATGGGCAGAATATTTGTTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGTAAGGGGGGGGGGGGGGGGACGGGTTAcg                          >  1:1270268/1‑149 (MQ=17)                                                                                                                           tGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGCACCCTAGTCTTATCTTCTCCACCTTCCACTCTTTGTCAGAATACGTGGCTGGTGATTCGCTTGAGCGTTCTGTGACTTATGGAAGTGATGGTGGTGGGGGAAGGATTACTCa                        <  2:574744/150‑1 (MQ=9)                                                                                                                           tGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCTTCCTTCCCCGCTTGGGCAGAATATTTGTTTGCGGATTCGCTTGAGATTTCAGGGACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCa                        <  2:466831/150‑1 (MQ=14)                                                                                                                           tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCTGATTCGCTTGAGTGTTCAGGGACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCa                        <  2:904611/150‑1 (MQ=255)                                                                                                                           tGCCTGCGGCAACGCTCTCTCGCGGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAAACTTCCCCGCAGGGGCAGAATATTTTATTGCGGATTAGCTTGTGAGTTCAGGGACTTTTGAAAGGGATGGTGGTGGGGGAAGGATTAcaaa                        >  1:266749/1‑147 (MQ=14)                                                                                                                               tGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGAATTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAgcgc                    >  2:1493177/1‑150 (MQ=255)                                                                                                                               tGCGGCAACGCTCTCTCGCTGACGCGCGAGCCGAACCTTAGTCGAAGCTTATCATCCTTCCCCCCCTGGGCAGAATAGTTTATTGCGGATTCGCTGGAGAGTTAAGGGACTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTAAgcgc                    >  1:1121047/1‑150 (MQ=14)                                                                                                                                 cGGGAACGCCCCCTCGCTGTCGCTCGAGCCGAACCTCGGCCGAAGCGCGACATCGGTACCCGCTTGGGCAGCCGTCATGATTGGGTAGTCGCTTGAGCGTTCTGGTACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTg                  <  2:792510/150‑1 (MQ=2)                                                                                                                                        gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACCGCTGCGCGCtt           >  1:1285562/1‑150 (MQ=255)                                                                                                                                                 cTGACTCTCTAGCCGAACCTTTGTCTAAGCTTCTGATCCTTCCCCGCTTTGGCAGAATATTTTATTGCGTATTCGCTTGAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCg  <  2:160809/150‑1 (MQ=14)                                                                                                                                                    |                                                                                                                                                   CATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCG  >  NC_000913/1286818‑1287111

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑