BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,335,654	T→A	5.3%	intergenic (‑64/+40)	eptA ← / ← adiC	lipid A phosphoethanolamine transferase/arginine:agmatine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,335,654	0	T	A	5.3%	439.3 / 9.3	127	intergenic (‑64/+40)	eptA/adiC	lipid A phosphoethanolamine transferase/arginine:agmatine antiporter
Reads supporting (aligned to +/- strand): ref base T (57/63); new base A (4/3); total (61/66)
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

Only 100 of 128 total aligned reads displayed.

AGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCAC  >  NC_000913/4335508‑4335799

aGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGtac                                                                                                                                                 <  4:897716/149‑3 (MQ=255)

 ggCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCt                                                                                                                                                <  3:303708/149‑1 (MQ=255)

    aGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTaa                                                                                                                                             <  3:685176/149‑1 (MQ=255)

        agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGtt                                                                                                                                         >  3:366653/1‑149 (MQ=255)

        agagATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGt                                                                                                                                          <  1:537206/148‑1 (MQ=255)

           gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTgg                                                                                                                                       >  2:150790/1‑148 (MQ=255)

           gaTATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGc                                                                                                                                      <  3:188338/149‑1 (MQ=255)

           gaTATAAAAAGCGGCCAACAATAACCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCCAACAAAGTGCTCATATGCTCGCAAACCAAAATTAAGCCAACCTTAAGTTCTTAAGGGTGGc                                                                                                                                      <  4:315014/149‑1 (MQ=255)

             tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCtt                                                                                                                                    >  3:230391/1‑149 (MQ=255)

             tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCt                                                                                                                                     <  1:406575/148‑1 (MQ=255)

             tataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCt                                                                                                                                     <  1:167817/148‑1 (MQ=255)

              ataAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCtt                                                                                                                                    >  1:270469/1‑148 (MQ=255)

                aaaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCtttt                                                                                                                                  <  2:806359/148‑1 (MQ=255)

                 aaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTa                                                                                                                                 <  1:905068/148‑1 (MQ=255)

                 aaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTa                                                                                                                                 >  2:911489/1‑148 (MQ=255)

                 aaaaGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTAt                                                                                                                                >  4:7648/1‑149 (MQ=255)

                     gCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGttt                                                                                                                            <  3:321282/149‑1 (MQ=255)

                     gCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGttt                                                                                                                            <  3:689596/149‑1 (MQ=255)

                        gCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGc                                                                                                                          <  2:579789/148‑1 (MQ=255)

                         ccAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAgg                                                                                                                                           >  1:148908/1‑130 (MQ=255)

                         ccAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGAACTTAAgg                                                                                                                                           <  2:148908/130‑1 (MQ=255)

                                tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTaa                                                                                                                  <  1:467205/148‑1 (MQ=255)

                                tAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAg                                                                                                                 <  3:336574/149‑1 (MQ=255)

                                  gCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAgg                                                                                                                <  1:493571/148‑1 (MQ=255)

                                     aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGc                                                                                                                                      <  2:237415/123‑1 (MQ=255)

                                     aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGc                                                                                                                                      >  1:237415/1‑123 (MQ=255)

                                     aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGaa                                                                                                            >  4:878605/1‑149 (MQ=255)

                                     aGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGa                                                                                                             <  1:510085/148‑1 (MQ=255)

                                      ggCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGaa                                                                                                            >  2:893257/1‑148 (MQ=255)

                                      ggCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGaa                                                                                                            >  2:622819/1‑148 (MQ=255)

                                       gCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAAt                                                                                                           <  2:30460/148‑1 (MQ=255)

                                          aGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTaa                                                                                                        <  2:21237/148‑1 (MQ=255)

                                          aGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTaa                                                                                                        <  1:55259/148‑1 (MQ=255)

                                              aaTTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCtt                                                                                                    <  2:137293/148‑1 (MQ=255)

                                                ttCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGc                                                                                                 >  3:689241/1‑149 (MQ=255)

                                                  cAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCt                                                                                                >  2:358357/1‑148 (MQ=255)

                                                     agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGc                                                                                                                          >  1:99752/1‑119 (MQ=255)

                                                     agagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGc                                                                                                                          <  2:99752/119‑1 (MQ=255)

                                                      gagGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTAtt                                                                                            >  2:788176/1‑148 (MQ=255)

                                                            cTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCa                                                                                      >  2:87094/1‑148 (MQ=255)

                                                                ttAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTaa                                                                                 <  4:689241/149‑1 (MQ=255)

                                                                      ggCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGAt                                                                            >  1:273100/1‑148 (MQ=255)

                                                                      ggCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATa                                                                           >  3:41871/1‑149 (MQ=255)

                                                                      ggCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATa                                                                           <  4:472909/149‑1 (MQ=255)

                                                                       gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTaa                                                                                 >  3:513386/1‑142 (MQ=255)

                                                                       gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTaa                                                                                 <  4:513386/142‑1 (MQ=255)

                                                                       gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATa                                                                           >  1:452902/1‑148 (MQ=255)

                                                                       gcgcTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATAc                                                                          <  3:165187/149‑1 (MQ=255)

                                                                         gcTTCAACATCACGTTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAATCAAAACTTAGCCAACCTTACGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACg                                                                         <  1:695760/148‑1 (MQ=255)

                                                                           ttCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGt                                                                      <  3:146784/149‑1 (MQ=255)

                                                                             cAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGttt                                                                    <  3:919916/149‑1 (MQ=255)

                                                                               aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGtttt                                                                   <  1:781586/148‑1 (MQ=255)

                                                                               aCATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGt                                                                      >  1:616004/1‑145 (MQ=255)

                                                                                cATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                         >  4:700437/1‑125 (MQ=255)

                                                                                cATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGt                                                                                         <  3:700437/125‑1 (MQ=255)

                                                                                 aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                                                                 <  2:286846/148‑1 (MQ=255)

                                                                                 aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                                                                 <  2:273100/148‑1 (MQ=255)

                                                                                 aTCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTa                                                                 <  1:224851/148‑1 (MQ=255)

                                                                                     cGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTa                                                             >  1:484810/1‑148 (MQ=255)

                                                                                       gtgtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGcc                                                          <  4:633760/149‑1 (MQ=255)

                                                                                        tgtTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGCTGGCTTTTATGTTTGCTGGCTTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGcc                                                          >  1:598117/1‑148 (MQ=255)

                                                                                           ttCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGtt                                                      <  3:750881/149‑1 (MQ=255)

                                                                                               aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTTGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                  <  4:303788/149‑1 (MQ=255)

                                                                                               aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                  >  3:109084/1‑149 (MQ=255)

                                                                                               aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                  >  3:400561/1‑149 (MQ=255)

                                                                                               aTCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTaa                                                  <  4:22293/149‑1 (MQ=255)

                                                                                                 cGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATg                                                               <  2:846089/134‑1 (MQ=255)

                                                                                                 cGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATg                                                               >  1:846089/1‑134 (MQ=255)

                                                                                                 cGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAAt                                                 >  2:668874/1‑148 (MQ=255)

                                                                                                   aaCAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCa                                              <  4:639971/149‑1 (MQ=255)

                                                                                                     cAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAg                                             >  2:725647/1‑148 (MQ=255)

                                                                                                      aaaGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAgcg                                           >  3:301533/1‑149 (MQ=255)

                                                                                                          tGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCAt                                        >  2:396315/1‑148 (MQ=255)

                                                                                                               atatGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATgg                                  >  3:821377/1‑149 (MQ=255)

                                                                                                                  tGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCg                                <  1:725647/148‑1 (MQ=255)

                                                                                                                     tCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtga                            >  3:462449/1‑149 (MQ=255)

                                                                                                                      tgcaaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGCTTTTATGTAGCCGGTTGCAATTCAGCGCATACATGGCGGtga                            >  2:923057/2‑148 (MQ=255)

                                                                                                                      cGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtga                            <  1:160474/148‑1 (MQ=255)

                                                                                                                      cGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtga                            >  2:138925/1‑148 (MQ=255)

                                                                                                                      cGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtga                            <  2:500567/148‑1 (MQ=255)

                                                                                                                      cGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtga                            >  2:799164/1‑148 (MQ=255)

                                                                                                                         agagcaaaATTAAGCCAACCTAAAGTTCTTGAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGAc                        <  3:777290/147‑1 (MQ=255)

                                                                                                                         aaagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGtgatga                         <  1:893257/148‑1 (MQ=255)

                                                                                                                          aagcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGAcc                       <  3:298105/149‑1 (MQ=255)

                                                                                                                            gcaaaATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCa                      >  1:317271/1‑148 (MQ=255)

                                                                                                                               aaaTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAg                  >  4:329446/1‑149 (MQ=255)

                                                                                                                                aaTTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAgg                 <  4:301533/149‑1 (MQ=255)

                                                                                                                                  ttAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGt               <  4:501233/149‑1 (MQ=255)

                                                                                                                                   tAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGt               >  1:856210/1‑148 (MQ=255)

                                                                                                                                     aGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTATATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGa             >  1:170737/1‑148 (MQ=255)

                                                                                                                                      gCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACa           >  4:243773/1‑149 (MQ=255)

                                                                                                                                      gCCAACCTTAAGAACTTAAGGTTGGCCGTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATATAAGGGCTACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCGTACATGGCGGTGATGACCATCAGGGTGAc            >  2:519311/1‑148 (MQ=255)

                                                                                                                                         aaCCTTAAGAACTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACaaa         >  1:362083/1‑148 (MQ=255)

                                                                                                                                           ccTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATg       <  2:215478/148‑1 (MQ=255)

                                                                                                                                           ccTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATg       <  2:551923/148‑1 (MQ=255)

                                                                                                                                           ccTTAAGAACTTAAGGTTGGC‑TTAAT‑TTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAc     >  2:892660/1‑148 (MQ=255)

                                                                                                                                            cTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCAt                     <  2:329252/133‑1 (MQ=255)

                                                                                                                                            cTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCAt                     >  1:329252/1‑133 (MQ=255)

                                                                                                                                             ttAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAcc    <  4:230391/149‑1 (MQ=255)

                                                                                                                                             ttAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGAc     >  1:602892/1‑148 (MQ=255)

AGGCAGATAGAGATATAAAAAGCGGCCAACAATAGCCAGGCGAGTAAATTCAAAGAGGGTCTTTTTAGTAGGCGCTTCAACATCACGGTGTTTCCATCGAACAAAGTGCGCATATGCTCGCAAAGCAAAATTAAGCCAACCTTAAGTTCTTAAGGTTGGCTTTTATGTTTGCTGGATTAAGGGGAATTAATCTTTGCTTATTGGTGCATCTAAGGGATACGGGTTTTTATGTAGCCGGTTGTAATTCAGCGCATACATGGCGGTGATGACCATCAGGGTGACAAATGACCAC  >  NC_000913/4335508‑4335799

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 16 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 38 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`