BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,820,243:1	+T	7.0%	coding (928/1683 nt)	ligB ←	DNA ligase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,820,243	1	.	T	7.0%	612.5 / 42.7	157	coding (928/1683 nt)	ligB	DNA ligase B
Reads supporting (aligned to +/- strand): ref base . (84/62); new base T (4/7); total (88/69)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 158 total aligned reads displayed.

AATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCA  >  NC_000913/3820094‑3820391

aaTATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTg                                                                                                                                                       >  2:1405188/1‑150 (MQ=255)

    tCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCt                                                                                                                                                  >  1:1414460/1‑150 (MQ=255)

     cccACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTa                                                                                                                                                 >  2:28647/1‑150 (MQ=255)

         cTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAgg                                                                                                                                             >  1:1187011/1‑150 (MQ=255)

           ccTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGtt                                                                                                                                           >  2:1422575/1‑150 (MQ=255)

              gCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCT‑CAGGTTgcc                                                                                                                                        <  1:352039/150‑3 (MQ=255)

               ccAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCT‑CAGGTTgccg                                                                                                                                       <  1:974546/150‑4 (MQ=255)

                cAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGatat                                                                                                                                      >  2:1240033/1‑150 (MQ=255)

                 aGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCT‑CAGGTTgccgaa                                                                                                                                     <  1:1310244/150‑6 (MQ=255)

                     ccTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCa                                                                                                                                 <  1:276702/150‑1 (MQ=255)

                      cTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAg                                                                                                                                <  1:491754/150‑1 (MQ=255)

                        gACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACAGCTCTTACACACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGc                                                                                                                              >  2:875408/1‑150 (MQ=255)

                          cGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGcca                                                                                                                            >  2:463635/1‑150 (MQ=255)

                              aCCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAg                                                                                                                        <  2:177826/150‑1 (MQ=255)

                              aCCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAg                                                                                                                        <  2:236698/150‑1 (MQ=255)

                               ccAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCT‑CAGGTTgcc                                                                                                                                        >  2:1074661/1‑131 (MQ=255)

                                cAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGcc                                                                                                                      >  2:27934/1‑150 (MQ=255)

                                 aaTATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCa                                                                                                                     <  2:1068230/150‑1 (MQ=255)

                                 aaTATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCa                                                                                                                     >  2:781799/1‑150 (MQ=255)

                                  atatTCACCCGCTGGGCTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCAc                                                                                                                    >  2:728050/1‑150 (MQ=255)

                                  atatTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCAc                                                                                                                    <  1:10187/150‑1 (MQ=255)

                                      tCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTg                                                                                                                <  2:430399/150‑1 (MQ=255)

                                       cACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCCCCAGCCACTCTGc                                                                                                               >  2:1367174/1‑150 (MQ=255)

                                        aCCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGcc                                                                                                              <  1:1426528/150‑1 (MQ=255)

                                        aCCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGcc                                                                                                              <  2:1187011/150‑1 (MQ=255)

                                          ccGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTg                                                                                                            <  1:1173369/150‑1 (MQ=255)

                                           cGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGc                                                                                                                       <  2:486370/138‑1 (MQ=255)

                                           cGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGc                                                                                                                       >  1:486370/1‑138 (MQ=255)

                                            gCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGc                                                                                                          <  1:126786/150‑1 (MQ=255)

                                            gCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGc                                                                                                          <  1:1410074/150‑1 (MQ=255)

                                                 aCTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGt                                                                                                     >  2:414566/1‑150 (MQ=255)

                                                  cTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTa                                                                                                    >  2:382692/1‑150 (MQ=255)

                                                  cTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTa                                                                                                    <  2:1030431/150‑1 (MQ=255)

                                                  cTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTa                                                                                                    >  1:742708/1‑150 (MQ=255)

                                                   ttttttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAg                                                                                                   <  1:763227/150‑1 (MQ=255)

                                                             tcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCgg                                                                                         >  1:946965/1‑150 (MQ=255)

                                                             tcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCgg                                                                                         <  2:1115203/150‑1 (MQ=255)

                                                              ccAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCggg                                                                                        >  2:1199365/1‑150 (MQ=255)

                                                               cAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGa                                                                                       >  1:595236/1‑150 (MQ=255)

                                                                aGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGAt                                                                                      <  2:497359/150‑1 (MQ=255)

                                                                   aTGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCt                                                                                   <  1:27934/150‑1 (MQ=255)

                                                                    tGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTg                                                                                  >  1:851146/1‑150 (MQ=255)

                                                                    tGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTACCTGGCCCGGTAGCCAATGGCGGGATTCTg                                                                                  >  2:1109399/1‑150 (MQ=255)

                                                                      aCAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGc                                                                                <  1:463635/150‑1 (MQ=255)

                                                                      aCAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGc                                                                                >  2:323542/1‑150 (MQ=255)

                                                                      aCAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGc                                                                                <  2:258218/150‑1 (MQ=255)

                                                                       cAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGct                                                                               >  2:977833/1‑150 (MQ=255)

                                                                         ggTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGctct                                                                             >  1:1313741/1‑150 (MQ=255)

                                                                           tGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCttt                                                                           >  1:1404843/1‑150 (MQ=255)

                                                                            gCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTc                                                                          >  2:85808/1‑150 (MQ=255)

                                                                            gCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTc                                                                          >  2:590413/1‑150 (MQ=255)

                                                                               aGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcgcc                                                                       >  1:1112101/1‑150 (MQ=255)

                                                                                 tGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTgg                                                                                 >  2:560409/1‑138 (MQ=255)

                                                                                 tGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTgg                                                                                 <  1:560409/138‑1 (MQ=255)

                                                                                 tGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcgccgc                                                                     >  1:70150/1‑150 (MQ=255)

                                                                                   aCGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCAGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTc                                                                   >  2:1238343/1‑150 (MQ=255)

                                                                                    cGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCg                                                                  <  1:1332775/150‑1 (MQ=255)

                                                                                    cGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCg                                                                  >  1:1259296/1‑150 (MQ=255)

                                                                                    cGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCg                                                                  <  2:742708/150‑1 (MQ=255)

                                                                                      cAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTa                                                                <  2:855142/150‑1 (MQ=255)

                                                                                           accGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACt                                                           >  1:971378/1‑150 (MQ=255)

                                                                                            ccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTAc                                                               <  1:1281443/145‑1 (MQ=255)

                                                                                            ccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTAc                                                               >  2:1281443/1‑145 (MQ=255)

                                                                                                       ccGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTg                                               <  1:1240033/150‑1 (MQ=255)

                                                                                                          ctctTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTTCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACg                                            >  1:114957/1‑150 (MQ=255)

                                                                                                            ctTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGaa                                          >  2:813625/1‑150 (MQ=255)

                                                                                                              tACCCACCGCACACTGATTTGCCTTCACTTCGGCAAGCTT‑ATCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTTCCTGGCCCGTTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAgg                                        <  2:494200/150‑1 (MQ=255)

                                                                                                              tACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAgg                                        <  2:660127/150‑1 (MQ=255)

                                                                                                               aCCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAggg                                       >  2:1299928/1‑150 (MQ=255)

                                                                                                               aCCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAggg                                       <  1:85808/150‑1 (MQ=255)

                                                                                                                cccACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGt                                      >  1:584392/1‑150 (MQ=255)

                                                                                                                   aCCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAAt                                   >  2:179778/1‑150 (MQ=255)

                                                                                                                   aCCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAAt                                   <  1:831133/150‑1 (MQ=255)

                                                                                                                   aCCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAAt                                   >  2:356601/1‑150 (MQ=255)

                                                                                                                      gCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTcc                                <  1:532825/150‑1 (MQ=255)

                                                                                                                        aaaCTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGc                              <  1:1199365/150‑1 (MQ=255)

                                                                                                                        aaaCTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGc                              >  2:924312/1‑150 (MQ=255)

                                                                                                                          aCTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCtt                            <  1:438918/150‑1 (MQ=255)

                                                                                                                          aCTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCtt                            >  1:1388010/1‑150 (MQ=255)

                                                                                                                           cTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCttt                           >  2:929578/1‑150 (MQ=255)

                                                                                                                             gAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTcc                         >  1:968248/1‑150 (MQ=255)

                                                                                                                             gAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTcc                         <  2:14595/150‑1 (MQ=255)

                                                                                                                              aaTTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTcca                        >  1:773608/1‑150 (MQ=255)

                                                                                                                               aTTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTccac                       <  2:595236/150‑1 (MQ=255)

                                                                                                                                  gCCTTCCCTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCCCCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCccccac                    >  2:206987/1‑150 (MQ=255)

                                                                                                                                  gCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCcaccac                    >  1:658772/1‑150 (MQ=255)

                                                                                                                                   ccTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACt                   >  2:409442/1‑150 (MQ=255)

                                                                                                                                    cTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTc                  >  1:708808/1‑150 (MQ=255)

                                                                                                                                    cTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTc                  >  1:215561/1‑150 (MQ=255)

                                                                                                                                     ttCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCa                 >  2:464265/1‑150 (MQ=255)

                                                                                                                                     ttCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCa                 >  1:213046/1‑150 (MQ=255)

                                                                                                                                      tCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCAt                <  2:575198/150‑1 (MQ=255)

                                                                                                                                       cACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCAtt               >  2:332386/1‑150 (MQ=255)

                                                                                                                                        aCTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTg              >  1:299426/1‑150 (MQ=255)

                                                                                                                                        aCTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTg              <  1:875408/150‑1 (MQ=255)

                                                                                                                                          aatatcAACCTGTAGCT‑CAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCg            >  2:126681/6‑150 (MQ=255)

                                                                                                                                           atatcAACCTGTAGCT‑CAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGt           >  1:96222/5‑150 (MQ=255)

                                                                                                                                             atcAACCTGTAGCT‑CAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAc         >  2:1205384/3‑150 (MQ=255)

                                                                                                                                                   cTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACgcgcgc   <  1:10646/150‑1 (MQ=255)

                                                                                                                                                    tG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCa  <  1:464265/150‑1 (MQ=255)

AATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCA  >  NC_000913/3820094‑3820391

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