BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,573,188	G→A	5.1%	pseudogene (775/1503 nt)	yjiT →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,188	0	G	A	5.1%	455.7 / 8.5	133	pseudogene (775/1503 nt)	yjiT	pseudogene
Reads supporting (aligned to +/- strand): ref base G (66/60); new base A (5/2); total (71/62)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

Only 100 of 133 total aligned reads displayed.

TTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCG  >  NC_000913/4573041‑4573333

ttAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                                                                                                                  >  3:404786/1‑149 (MQ=255)

 tAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                                                                                                                  >  2:262316/1‑148 (MQ=255)

 tAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGt                                                                                                                                                  <  1:213120/148‑1 (MQ=255)

   gTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTaga                                                                                                                                               <  4:630868/149‑1 (MQ=255)

          ctCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTaa                                                                                                                                         <  2:428757/148‑1 (MQ=255)

            cacaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAgg                                                                                                                                             <  1:194265/142‑1 (MQ=255)

             acaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCt                                                                                                                                     >  4:181980/1‑149 (MQ=255)

             acaACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCt                                                                                                                                     >  3:283088/1‑149 (MQ=255)

              caACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAACCCt                                                                                                                                     <  2:91801/148‑1 (MQ=255)

                aCCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTa                                                                                                                                          <  1:501837/141‑1 (MQ=255)

                aCCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTa                                                                                                                                          >  2:501837/1‑141 (MQ=255)

                  cTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGAc                                                                                                                                 >  2:525926/1‑148 (MQ=255)

                   taataaATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAGACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTa                                                                                                                               >  3:556355/1‑149 (MQ=255)

                    aataaATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTaga                                                                                                                                               <  1:140590/132‑1 (MQ=255)

                      taaATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaa                                                                                                                            >  4:52157/1‑149 (MQ=255)

                       aaaTCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaa                                                                                                                            >  1:63877/1‑148 (MQ=255)

                        aaTCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaa                                                                                                                          >  4:421352/1‑149 (MQ=255)

                         aTCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaa                                                                                                                          >  2:760488/1‑148 (MQ=255)

                           cacaAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACaaaaaaa                                                                                                                       <  3:181980/149‑1 (MQ=255)

                             caAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAAt                                                                                                                      >  1:99568/1‑148 (MQ=255)

                              aaaGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                    <  3:450096/149‑1 (MQ=255)

                              aaaGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                    <  3:584571/149‑1 (MQ=255)

                              aaaGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATaa                                                                                                                    >  3:627742/1‑149 (MQ=255)

                                aGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAaga                                                                                                                  >  4:346560/1‑149 (MQ=255)

                                 gAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAagag                                                                                                                 <  3:509883/149‑1 (MQ=255)

                                 gAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAaga                                                                                                                  <  1:134316/148‑1 (MQ=255)

                                     cAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTAcg                                                                                                             >  4:177536/1‑149 (MQ=255)

                                         acaacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGtt                                                                                                         <  3:421352/149‑1 (MQ=255)

                                          caacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAAGGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTg                                                                                                        >  3:596113/1‑149 (MQ=255)

                                           aacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGa                                                                                                                                  >  4:340884/1‑122 (MQ=255)

                                           aacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGa                                                                                                                                  <  3:340884/122‑1 (MQ=255)

                                           aacaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGc                                                                                                       >  4:368408/1‑149 (MQ=255)

                                            acaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAag                                                                                                                   <  1:355690/136‑1 (MQ=255)

                                            acaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTg                                                                                                        <  2:884490/147‑1 (MQ=255)

                                            acaGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGc                                                                                                       >  1:147003/1‑148 (MQ=255)

                                                 ttGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGc                                                                                                       >  3:698464/1‑143 (MQ=255)

                                                 ttGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGc                                                                                                       <  4:698464/143‑1 (MQ=255)

                                                 ttGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAAtt                                                                                                  >  1:373317/1‑148 (MQ=255)

                                                 ttGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGATAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGAAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTc                                                                                                 >  4:659864/1‑149 (MQ=255)

                                                   gagaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCt                                                                                                <  1:760488/148‑1 (MQ=255)

                                                   gagaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTa                                                                                               >  4:679039/1‑149 (MQ=255)

                                                     gaAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACCGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTAt                                                                                              <  1:904096/148‑1 (MQ=255)

                                                      aaGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTg                                                                                            <  3:52157/149‑1 (MQ=255)

                                                        ggTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTgg                                                                                           <  2:298470/148‑1 (MQ=255)

                                                          ttCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTc                                                                                         <  1:623083/148‑1 (MQ=255)

                                                              cccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAaga                                                                                                                  <  3:409403/119‑1 (MQ=255)

                                                              cccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAaga                                                                                                                  >  4:409403/1‑119 (MQ=255)

                                                               ccTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaaa                                                                                   >  3:573051/1‑149 (MQ=255)

                                                                cTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaaa                                                                                   >  2:902341/1‑148 (MQ=255)

                                                                cTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGaaaa                                                                                   >  2:575370/1‑148 (MQ=255)

                                                                  aCCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGc                                                                                 <  1:471152/148‑1 (MQ=255)

                                                                     tGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATcc                                                                             >  3:774096/1‑149 (MQ=255)

                                                                       gCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATACGAGTACGCTTTGCCAATTCTATTGGTCAGAAAAGCATCCAg                                                                           >  4:541937/1‑149 (MQ=255)

                                                                         gAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGa                                                                          <  2:413994/148‑1 (MQ=255)

                                                                           aGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAACCATCCAGATGa                                                                       <  3:794172/149‑1 (MQ=255)

                                                                             cGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAg                                                                      <  2:147003/148‑1 (MQ=255)

                                                                                 tttCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTa                                                                  >  2:889266/1‑148 (MQ=255)

                                                                                 tttCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAg                                                                 <  3:765123/149‑1 (MQ=255)

                                                                                    ccAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAggg                                                               >  1:798881/1‑148 (MQ=255)

                                                                                    ccAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGt                                                              >  3:740453/1‑149 (MQ=255)

                                                                                                   gatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTc                                                >  3:878925/1‑148 (MQ=255)

                                                                                                   gatgaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTc                                                <  4:878925/148‑1 (MQ=255)

                                                                                                      gaAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCgg                                             >  1:48113/1‑148 (MQ=255)

                                                                                                       aaaCAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGa                                            >  2:778744/1‑148 (MQ=255)

                                                                                                          cAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGa                                         <  2:179210/148‑1 (MQ=255)

                                                                                                          cAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGa                                         <  2:99568/148‑1 (MQ=255)

                                                                                                               aCGCATTTTTTGAATGGACTGCTATGCTCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttc                                   <  3:621749/149‑1 (MQ=255)

                                                                                                               aCGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttc                                   <  3:177536/149‑1 (MQ=255)

                                                                                                                 gCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcct                                 <  4:774096/149‑1 (MQ=255)

                                                                                                                  cATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcct                                 <  1:902341/148‑1 (MQ=255)

                                                                                                                   aTTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttccttt                               >  3:736948/1‑149 (MQ=255)

                                                                                                                    ttttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcctttc                              <  4:760751/149‑1 (MQ=255)

                                                                                                                    ttttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcctttc                              <  3:346560/149‑1 (MQ=255)

                                                                                                                     tttttGAATGGACTGCTATGCGCAGCCTGTGTAGAGGCTAAGCCTCGACTACAAAAAAATAaga                                                                                                                  <  3:132156/64‑1 (MQ=255)

                                                                                                                     tttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAaga                                                                                                                  >  4:132156/1‑64 (MQ=255)

                                                                                                                     tttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcctttc                              >  2:877974/1‑148 (MQ=255)

                                                                                                                      ttttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTttcctttcct                            <  4:499363/149‑1 (MQ=255)

                                                                                                                        ttGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGt                          <  4:340915/149‑1 (MQ=255)

                                                                                                                            aTGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTACTGTAAc                       >  1:834379/1‑148 (MQ=255)

                                                                                                                                aCTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTa                                                                  <  2:291604/101‑1 (MQ=255)

                                                                                                                                aCTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTa                                                                  >  1:291604/1‑101 (MQ=255)

                                                                                                                                aCTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCg                  <  4:596113/149‑1 (MQ=255)

                                                                                                                                aCTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCg                  >  3:306249/1‑149 (MQ=255)

                                                                                                                                 cTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCg                  <  2:63877/148‑1 (MQ=255)

                                                                                                                                 cTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCg                  <  2:798881/148‑1 (MQ=255)

                                                                                                                                      aTGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAg                    >  2:176351/1‑141 (MQ=255)

                                                                                                                                      aTGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAg                    <  1:176351/141‑1 (MQ=255)

                                                                                                                                          ttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcctttc                              >  3:365388/3‑127 (MQ=255)

                                                                                                                                          ttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcctttc                              <  4:365388/125‑1 (MQ=255)

                                                                                                                                          ttaGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAAc        >  3:248815/3‑149 (MQ=255)

                                                                                                                                          gcAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcc                                  <  2:666578/123‑1 (MQ=255)

                                                                                                                                           taGCCTCTACAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAAc        >  1:882356/2‑148 (MQ=255)

                                                                                                                                             gCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCAt                                                                               <  1:631139/75‑1 (MQ=255)

                                                                                                                                             gCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTa      >  1:765817/1‑148 (MQ=255)

                                                                                                                                               ctctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcc                                  >  4:227185/1‑118 (MQ=255)

                                                                                                                                               ctctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcc                                  <  3:227185/118‑1 (MQ=255)

                                                                                                                                                tctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGtttcctttcc                             >  1:770782/1‑122 (MQ=255)

                                                                                                                                                tctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCa          >  2:384604/1‑141 (MQ=255)

                                                                                                                                                tctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCa          <  1:384604/141‑1 (MQ=255)

                                                                                                                                                tctGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTAcgcg  >  3:871814/1‑149 (MQ=255)

TTAGTACTTACTCACAACCTAATAAATCACAAAGAACCAGTACAACAGCTTGAGAAGGTTCACCCTACCTGGCGAAGCGAGTTTCCAATACCGTTAGATGATGAAACAGGAACGCATTTTTTGAATGGACTGCTATGCGCAGCCTCTGTAGAGGCTAAGCCTCGACTACAAAAAAATAAGAGTACGCGTTGCCAATTCTATTGGTCAGAAAAGCATCCAGATGAGCTTAGGGTGATAGTATCTCTTCCGGACGAAGTTTCCTTTCCTGTAACAAGCGAGCCGTCAACTACGCG  >  NC_000913/4573041‑4573333

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