Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,599,755	C→A	5.6%	intergenic (‑97/+174)	livJ ← / ← rpoH	branched‑chain amino acid ABC transporter periplasmic binding protein/RNA polymerase, sigma 32 (sigma H) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,599,755	0	C	A	5.6%	218.4 / 3.1	71	intergenic (‑97/+174)	livJ/rpoH	branched‑chain amino acid ABC transporter periplasmic binding protein/RNA polymerase, sigma 32 (sigma H) factor
Reads supporting (aligned to +/- strand):  ref base C (38/29);  new base A (2/2);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

GAATGCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCA  >  NC_000913/3599608‑3599904                                                                                                                                                    |                                                                                                                                                      gAATGCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCAt                                                                                                                                                     >  2:1865387/1‑150 (MQ=255)      ccAGCGCGATACATCCTTCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActct                                                                                                                                                <  1:772006/150‑1 (MQ=255)      ccAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCCTTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActct                                                                                                                                                >  1:1725937/1‑150 (MQ=255)      ccAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActct                                                                                                                                                <  1:1017297/150‑1 (MQ=255)       cAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTg                                                                                                                                               >  1:674889/1‑150 (MQ=255)        ctcgcgATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATGCTTCTGGTTATTACGGGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGc                                                                                                                                              <  2:1841984/148‑1 (MQ=255)             gATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTACAGGTTATTACGTGTGTTGTGATGTTGTATTTCAGCACTTTATTTCGATTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTaaaa                                                                                                                                         >  1:586102/1‑150 (MQ=255)              ttACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAAc                                                                                                                                        <  2:1890773/149‑1 (MQ=255)               tACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACa                                                                                                                                       <  2:1780552/150‑1 (MQ=255)               tACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACa                                                                                                                                       >  2:1129502/1‑150 (MQ=255)                 cATCCTGCCAGTAACGCTTTACCCTTTATGATCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGAGTATGCTGCTAAAGCacg                                                                                                                                     <  2:1953080/150‑3 (MQ=255)                  aTCCTGCCAGTAACGCTTTACCCTTTAAGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATAc                                                                                                                                    <  1:1102821/150‑1 (MQ=255)                   tCCTTCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATAcc                                                                                                                                   <  1:1796575/150‑1 (MQ=255)                    ccTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATAccc                                                                                                                                  >  2:480804/1‑150 (MQ=255)                     cTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGAGTATGCTGCTAAAGCacgggta                                                                                                                                 <  2:660142/150‑7 (MQ=255)                       gCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAt                                                                                                                               >  2:328962/1‑150 (MQ=255)                        ccAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAtt                                                                                                                              >  1:174689/1‑150 (MQ=255)                          aGTAACTCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAtttt                                                                                                                            <  1:1522783/150‑1 (MQ=255)                          aGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAtttt                                                                                                                            <  1:865298/150‑1 (MQ=255)                           gTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAttttt                                                                                                                           <  2:927182/150‑1 (MQ=255)                                      cccTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTACAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGa                                                                                                                >  1:76458/1‑150 (MQ=255)                                       ccTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGaa                                                                                                               >  1:1642510/1‑150 (MQ=255)                                       ccTTTATGTTCATCCTGAGAATCCACATTCTTCTGGTTTTTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGaa                                                                                                               <  2:1233488/150‑1 (MQ=255)                                        cTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAAt                                                                                                              <  2:548522/150‑1 (MQ=255)                                          ttATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTTTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACAACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAg                                                                                                            <  2:83356/150‑1 (MQ=255)                                          ttATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAg                                                                                                            >  2:193042/1‑150 (MQ=255)                                            aTGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTATATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATAGGAATATCt                                                                                                          >  1:753394/1‑150 (MQ=255)                                              gTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTAt                                                                                                        >  2:557185/1‑150 (MQ=255)                                              gTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTAt                                                                                                        >  1:1766542/1‑150 (MQ=255)                                                    ccTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGAc                                                                                                  <  2:955752/150‑1 (MQ=255)                                                     cTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACa                                                                                                 <  2:1746006/150‑1 (MQ=255)                                                           aTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttat                                                                                           >  2:529046/1‑150 (MQ=255)                                                           aTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttat                                                                                           <  1:1983816/150‑1 (MQ=255)                                                             ccccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttatta                                                                                         <  2:29117/150‑1 (MQ=255)                                                              cccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTaa                                                                                        >  1:834761/1‑150 (MQ=255)                                                              cccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTaa                                                                                        >  2:513387/1‑150 (MQ=255)                                                                      tctGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGc                                                                                >  1:2025711/1‑150 (MQ=255)                                                                         ggTTATTACGTGTGTTGTGATGTTGTTATTCAGCACTTTATTTCGTTTTATGAATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATATCTATTTTGACAGTTTATTAACAATCTGCGTg                                                                             >  1:558171/1‑150 (MQ=255)                                                                            tattaCGTGTGTTGTGATGGTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTgggg                                                                          <  2:720480/150‑1 (MQ=255)                                                                                 cGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTgg                                                                     >  2:1177117/1‑150 (MQ=255)                                                                                 cGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTCCTTTAGCAGCATACTCTGCTAAAACATACCCGATATTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTgt                                                                     >  2:1605632/1‑149 (MQ=255)                                                                                     tgttgtGATGTTGTTTTTCAGCACTTTATTTCGTGTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGtt                                                                 <  2:1487849/150‑1 (MQ=255)                                                                                     tgttggGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGtt                                                                 <  2:1444140/150‑1 (MQ=255)                                                                                        tgtgATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTg                                                              >  1:1606193/1‑150 (MQ=255)                                                                                        tgtgATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTg                                                              >  2:1164816/1‑150 (MQ=255)                                                                                            aTGGTGGTTTTCAGCACTTTATTTCGTTATATGCATGCCGCCCCGGGCTTTCGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATAATGGAATAGCTATTTTGACAGTGTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCgg                                                          <  1:967220/150‑1 (MQ=255)                                                                                             tgttgtTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGa                                                         >  2:356108/1‑150 (MQ=255)                                                                                                tgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAggg                                                      <  2:1766541/150‑1 (MQ=255)                                                                                                 gttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATTCTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGGTTATTTACAATCTGCGTGGGGATTGGCGTTTTGCAGGAgggg                                                     >  2:1357586/1‑150 (MQ=255)                                                                                                 gttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAgggg                                                     >  2:725630/1‑150 (MQ=255)                                                                                                       cAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTg                                               <  2:174689/150‑1 (MQ=255)                                                                                                         gCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGat                                             <  2:1916486/150‑1 (MQ=255)                                                                                                         gCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGat                                             >  1:667507/1‑150 (MQ=255)                                                                                                           aCTTTATTTCGTTTTATGCATGCCTAACCGTTCTTTAGCAGCATACTCTGGTAAAACATAGCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGatta                                           <  2:393031/150‑1 (MQ=255)                                                                                                            cTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGattat                                          >  1:678884/1‑150 (MQ=255)                                                                                                                   tCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACagag                                   >  1:223456/1‑150 (MQ=255)                                                                                                                             cATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCaaaaaa                         >  2:1573591/1‑150 (MQ=255)                                                                                                                              aTGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaac                        <  1:2041370/150‑1 (MQ=255)                                                                                                                              aTGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGGTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaac                        <  2:1666628/150‑1 (MQ=255)                                                                                                                              aTGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAAGAGCTATTTTGACAGTTTATTAACAATCTGCGGGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaac                        <  1:341102/150‑1 (MQ=255)                                                                                                                               tGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaaca                       >  2:659207/1‑150 (MQ=255)                                                                                                                               tGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaaca                       >  1:1165419/1‑150 (MQ=255)                                                                                                                                    aCCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaacaaaaac                  >  1:1400093/1‑150 (MQ=255)                                                                                                                                    aCCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaacaaaaac                  >  1:462519/1‑150 (MQ=255)                                                                                                                                    aCCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCAGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaacaaaaac                  >  2:1424464/1‑150 (MQ=255)                                                                                                                                       cGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAAcccc               >  1:522183/1‑150 (MQ=255)                                                                                                                                        gTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCg              <  1:41763/150‑1 (MQ=255)                                                                                                                                          gttttAGCAGAGTATGCTGCTAAAACATACCCGTTTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCggc            >  1:198555/3‑149 (MQ=255)                                                                                                                                          gttttAGCAGAGTATGCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCggc            >  1:1740976/3‑149 (MQ=255)                                                                                                                                           cTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGAc           >  1:776408/1‑150 (MQ=255)                                                                                                                                            tttAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGAAATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGAct          >  1:1073814/1‑150 (MQ=255)                                                                                                                                               agcagcATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACtctc       <  2:1606192/150‑1 (MQ=255)                                                                                                                                                gcagcaTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCa      >  1:388048/1‑150 (MQ=255)                                                                                                                                                gcagcaTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCa      >  2:638764/1‑150 (MQ=255)                                                                                                                                                    caTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCa  <  2:674889/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GAATGCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCA  >  NC_000913/3599608‑3599904

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