BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,926,900	T→G	5.1%	D647A (GAC→GCC)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,926,900	0	T	G	5.1%	436.1 / 9.5	151	D647A (GAC→GCC)	ptrB	protease II
Reads supporting (aligned to +/- strand): ref base T (63/80); new base G (6/2); total (69/82)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

Only 100 of 153 total aligned reads displayed.

TAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATAC  >  NC_000913/1926751‑1927049

tAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                                                                                                                                                       >  1:896562/1‑150 (MQ=255)

 aaCACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                                                                                                                      <  1:1690550/150‑1 (MQ=255)

 aaCACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTc                                                                                                                                                      <  1:591114/150‑1 (MQ=255)

  acacTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGACTGAGTCCATGTCg                                                                                                                                                     >  1:52708/1‑150 (MQ=255)

  acacTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGAGCCAGCGCGACCAGAAAAGAATATTCCATCGCTACGACTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGACTTAGTCCAAGTCg                                                                                                                                                     >  2:629322/1‑150 (MQ=255)

  acacTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTACCCTGCGCCAGCGCGACCAGAAAAGAATATTACATCGCTACGACTTCGTACGATTTAAAGCGACAAGATTTACCGCCATGTCCTGAGTCCAAGTct                                                                                                                                                     >  2:637754/1‑149 (MQ=255)

   cacTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATGGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCTTGGCCTGAGTCCATGTCgg                                                                                                                                                    <  1:1103817/150‑1 (MQ=255)

   cacTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                                                                    <  1:76789/150‑1 (MQ=255)

   cacTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                                                                    <  2:1147651/150‑1 (MQ=255)

   cacTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTAAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCgg                                                                                                                                                    <  2:131989/150‑1 (MQ=255)

     cTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTa                                                                                                                                                  <  1:1803540/150‑1 (MQ=255)

       gAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACa                                                                                                                                                >  1:428744/1‑150 (MQ=255)

         aaCATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACaga                                                                                                                                              <  1:1983261/150‑1 (MQ=255)

         aaCATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGACATGGACTCAGGCCATGGCGGTaaatc                                                                                                                                              <  2:343743/150‑5 (MQ=255)

           cATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTTAGTCCATGTCGGTACAGAGc                                                                                                                                            <  1:380833/150‑1 (MQ=255)

           cATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGc                                                                                                                                            >  1:577191/1‑150 (MQ=255)

            aTTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCa                                                                                                                                           <  2:1345801/150‑1 (MQ=255)

              tATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCAGATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGCTACAGAGCatt                                                                                                                                         >  1:99904/1‑148 (MQ=255)

                 cTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaa                                                                                                                                      >  1:741907/1‑150 (MQ=255)

                  tGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                                                                                     <  2:1961258/150‑1 (MQ=255)

                  tGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaaa                                                                                                                                     >  1:1872680/1‑150 (MQ=255)

                      aaaTACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATg                                                                                                                                 <  2:1004916/150‑1 (MQ=255)

                       aaTACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                                                                                                                <  1:705833/150‑1 (MQ=255)

                       aaTACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATgg                                                                                                                                <  1:1254704/150‑1 (MQ=255)

                        aTACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGt                                                                                                                               >  2:1791390/1‑150 (MQ=255)

                           cTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATATAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGGAATAAAAGAAGGTCAt                                                                                                                            >  2:482847/1‑150 (MQ=255)

                              aGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCgg                                                                                                                         <  1:517476/150‑1 (MQ=255)

                              aGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGAGTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCgg                                                                                                                         >  1:1071777/1‑150 (MQ=255)

                                 ccGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGttt                                                                                                                      >  1:1362774/1‑150 (MQ=255)

                                   gCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTc                                                                                                                    <  1:165258/150‑1 (MQ=255)

                                    aaGGGGTTGCTGGTAATGTTCCCTGCGCAATAGCGACCAGAAAAGCATATTCCATCGCTTGGCCTTCGTACGATTTAAAGCGACCAGTTTTAGCGCCATGGCCTGAGTCCATGTCGGTACAGAGCTATATAAGATGGTCATCGGTTTTCa                                                                                                                   <  2:1780077/149‑1 (MQ=255)

                                      ggCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGAATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGc                                                                                                                 <  1:1752374/150‑1 (MQ=255)

                                         gTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATAATCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcg                                                                                                              <  1:1299453/150‑1 (MQ=255)

                                          tAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTAGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTcgc                                                                                                             >  2:419626/1‑150 (MQ=255)

                                                    tGTTCCCTTCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGc                                                                                                   <  2:946450/150‑1 (MQ=255)

                                                    tGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATGCCATCGCTACGCCTTCGTACGAGTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGc                                                                                                   <  1:306729/150‑1 (MQ=255)

                                                        cccTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGAcc                                                                                               >  2:1611784/1‑150 (MQ=255)

                                                           tGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAt                                                                                            >  2:1302767/1‑150 (MQ=255)

                                                              gcCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCAtttt                                                                                         <  2:797072/150‑1 (MQ=255)

                                                                cAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGc                                                                                       <  2:1169394/150‑1 (MQ=255)

                                                                 aGCGCGACCAGAAATGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGAACAGATTTACCGCCATGGCCTGAGTCCATGTCGTTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGcc                                                                                      <  2:1931745/150‑1 (MQ=255)

                                                                   cgcgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                                                                                    <  2:1871140/150‑1 (MQ=255)

                                                                   cgcgACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCgg                                                                                    >  2:143500/1‑150 (MQ=255)

                                                                      gACCAGAAAAGCATATTCCCTCGCTACGCATTCTTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGGCCATGTCGGTAAAGAGCAATAAAAGAAGGTCATCGGTGTTAAGCTCGCGAAATTTAGCGCCCCATGTTGACGGGTc                                                                                 >  2:316946/1‑150 (MQ=255)

                                                                      gACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGAATTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTc                                                                                 >  2:575365/1‑150 (MQ=255)

                                                                       aCCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCGTCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTcc                                                                                >  2:1720039/1‑150 (MQ=255)

                                                                       aCCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTcc                                                                                <  1:294542/150‑1 (MQ=255)

                                                                        ccAGAAAAGCATATTCCATCGCTCCGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc                                                                               <  1:1760461/150‑1 (MQ=255)

                                                                        ccAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTccc                                                                               <  1:2015064/150‑1 (MQ=255)

                                                                           gAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAat                                                                            <  2:1578963/150‑1 (MQ=255)

                                                                               aGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCACAATTTAGCGACCCATTTTGCCGGTTCCCAATATTg                                                                        <  2:1561501/150‑1 (MQ=255)

                                                                                gCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc                                                                       <  1:1032516/150‑1 (MQ=255)

                                                                                gCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGc                                                                       <  1:1605534/150‑1 (MQ=255)

                                                                                 cATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCa                                                                      >  1:500839/1‑150 (MQ=255)

                                                                                  atatTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCAc                                                                     >  1:181551/1‑150 (MQ=255)

                                                                                     ttCCATCGCTACGCCTTCGTACTATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGATCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTg                                                                  <  2:1090522/150‑1 (MQ=255)

                                                                                     ttCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTg                                                                  <  2:1747780/150‑1 (MQ=255)

                                                                                       ccATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTgag                                                                >  1:404815/1‑150 (MQ=255)

                                                                                       ccATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTgag                                                                >  1:1187238/1‑150 (MQ=255)

                                                                                       ccATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTgag                                                                >  1:1317695/1‑150 (MQ=255)

                                                                                            gCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCg                                                           >  2:466634/1‑150 (MQ=255)

                                                                                            gCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGACTCCATGTAGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCACGCAATTTAGAGACCCAATTTGACGGTTCCCAATATTGCACCTGAGAATCg                                                           >  1:1263770/1‑150 (MQ=255)

                                                                                             cTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCTTGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGt                                                          <  1:1572343/150‑1 (MQ=255)

                                                                                               aCGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGc                                                        <  1:1366390/150‑1 (MQ=255)

                                                                                               aCGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGCCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTATCAGCTCGCGCAATCTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGc                                                        <  2:1400383/150‑1 (MQ=255)

                                                                                                cGCCTTCGTACGATTTAAAGCGACCAGATTTACCGGGATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCa                                                       <  1:1861230/150‑1 (MQ=255)

                                                                                                cGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCTATAAAAGATGGTCATCTGTTTTCAGCTCGAGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCa                                                       <  1:154174/150‑1 (MQ=255)

                                                                                                cGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCa                                                       >  1:1488879/1‑150 (MQ=255)

                                                                                                   cTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAAc                                                    <  1:1302767/150‑1 (MQ=255)

                                                                                                      cGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCgg                                                 <  2:448405/150‑1 (MQ=255)

                                                                                                      cGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCgg                                                 >  1:1781197/1‑150 (MQ=255)

                                                                                                         aCGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCg                                              <  1:624775/150‑1 (MQ=255)

                                                                                                         aCGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCg                                              <  1:703548/150‑1 (MQ=255)

                                                                                                              ttAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTAcc                                         <  1:1471250/150‑1 (MQ=255)

                                                                                                                  aGCGACCACATTTACCGCCAAGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCAGCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTa                                     <  1:1231576/150‑1 (MQ=255)

                                                                                                                   gCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTaa                                    <  1:74581/150‑1 (MQ=255)

                                                                                                                    cGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTaaa                                   >  1:881041/1‑150 (MQ=255)

                                                                                                                     gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGCGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAAt                                  >  1:663628/1‑150 (MQ=255)

                                                                                                                     gACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCACAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAAt                                  >  1:570204/1‑150 (MQ=255)

                                                                                                                      aCCAGATTTACAGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCATCTCGCGCATTTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATg                                 <  1:629322/150‑1 (MQ=255)

                                                                                                                       ccAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGc                                >  1:1583602/1‑150 (MQ=255)

                                                                                                                        cacaTTTACCCCCATGGCCTGAGTCCATGTAGGAACAGAGCCATAAAACATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTGGCCGGTTCCCAATATTGCACCAGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCg                               <  2:1254959/147‑1 (MQ=255)

                                                                                                                        cAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCg                               <  1:1855771/150‑1 (MQ=255)

                                                                                                                         aGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCgg                              <  1:1553827/150‑1 (MQ=255)

                                                                                                                           aTTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGAt                            <  2:1915135/150‑1 (MQ=255)

                                                                                                                            tttACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATa                           <  2:741907/150‑1 (MQ=255)

                                                                                                                             gtACCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATaa                          >  2:1856087/2‑150 (MQ=255)

                                                                                                                              tACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGCCCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATaac                         >  1:1738258/1‑149 (MQ=255)

                                                                                                                               aCCGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGc                        >  2:1674362/1‑150 (MQ=255)

                                                                                                                                ccGACATGGACTCAGGCCATGGCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGGTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGcc                       >  2:144767/1‑150 (MQ=255)

                                                                                                                                  gCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCtg                     >  2:1093604/1‑150 (MQ=255)

                                                                                                                                  gCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCtg                     >  1:1696936/1‑150 (MQ=255)

                                                                                                                                       ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCgg                >  2:1389372/1‑150 (MQ=255)

                                                                                                                                       ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCgg                >  2:1397207/1‑150 (MQ=255)

                                                                                                                                       ggCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCgg                <  2:570769/150‑1 (MQ=255)

                                                                                                                                        gCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGt               <  1:1898871/150‑1 (MQ=255)

                                                                                                                                         gcTGAGTCCTTGTCGGTACAGAGCTATATATGATGGTCATCGGTATTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCTAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTg              <  1:449694/149‑1 (MQ=255)

                                                                                                                                         ccTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATTCGGATAAGCCAGTGCGGTg              >  1:791394/1‑150 (MQ=255)

                                                                                                                                                 caagTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTc      <  1:1791390/147‑1 (MQ=255)

                                                                                                                                                 cATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTc      >  1:1627843/1‑150 (MQ=255)

TAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAAAGATGGTCATCGGTTTTCAGCTCGCGCAATTTAGCGACCCATTTTGCCGGTTCCCAATATTGCACCTGAGAATCGTGCAAACCGGTCGTTACCAGTAAATGCGGATAAGCCTGTGCGGTGACGTTGTCATAC  >  NC_000913/1926751‑1927049

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