BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,458,870	G→A	T304T (ACG→ACA)	pmbA →	metalloprotease PmbA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,458,870	0	G	A	99.7%	1173.5 / NA	291	T304T (ACG→ACA)	pmbA	metalloprotease PmbA
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (147/143); total (147/144)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

AAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACGCCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCGGCGCGGAAAC  >  NC_000913/4458772‑4458967

aaaTCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACAcc                                                                                                 <  1:7615354/101‑1 (MQ=255)

aaaTCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACAcc                                                                                                 <  1:7615579/101‑1 (MQ=255)

aaaTCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACAcc                                                                                                 >  1:7616016/1‑101 (MQ=255)

 aaTCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCa                                                                                                <  2:7615666/101‑1 (MQ=255)

  aTCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCAt                                                                                               >  1:7616022/1‑101 (MQ=255)

  aTCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCAt                                                                                               >  2:7616020/1‑101 (MQ=255)

     tACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCg                                                                                            >  1:7616026/1‑101 (MQ=255)

     tACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCg                                                                                            <  2:7615319/101‑1 (MQ=255)

     tACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCg                                                                                            <  2:7615635/101‑1 (MQ=255)

      aCCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGa                                                                                           >  2:7616028/1‑101 (MQ=255)

        cTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACa                                                                                         >  1:7616032/1‑101 (MQ=255)

         ttCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAg                                                                                        <  1:7615507/101‑1 (MQ=255)

         ttCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAg                                                                                        >  2:7616033/1‑101 (MQ=255)

          tCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGc                                                                                       <  1:7615530/101‑1 (MQ=255)

            ctgctgGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGa                                                                                     <  2:7615570/101‑1 (MQ=255)

             tgctgGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGaa                                                                                    >  2:7616036/1‑101 (MQ=255)

              gctgGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAg                                                                                   <  2:7615599/101‑1 (MQ=255)

               ctgGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAgg                                                                                  >  1:7616039/1‑101 (MQ=255)

                tgGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGgt                                                                                 >  1:7616044/1‑101 (MQ=255)

                 ggATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAgg                                                                                  >  2:7616046/1‑99 (MQ=255)

                 ggATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGgtg                                                                                <  1:7615608/101‑1 (MQ=255)

                  gATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGgtgt                                                                               <  2:7615729/101‑1 (MQ=255)

                   aTTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGtgtg                                                                              <  1:7615482/101‑1 (MQ=255)

                    ttCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGa                                                                                     >  2:7616051/1‑93 (MQ=255)

                    ttCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTgc                                                                             <  2:7615417/101‑1 (MQ=255)

                       gCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCa                                                                          <  1:7615552/101‑1 (MQ=255)

                        cTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCAc                                                                         >  1:7616055/1‑101 (MQ=255)

                        cTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCAc                                                                         >  2:7616054/1‑101 (MQ=255)

                         tGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACa                                                                                         <  2:7615678/84‑1 (MQ=255)

                          gggTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCg                                                                       >  2:7616061/1‑101 (MQ=255)

                           ggTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGa                                                                      <  2:7615557/101‑1 (MQ=255)

                            gTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAg                                                                     >  1:7616063/1‑101 (MQ=255)

                            gTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAg                                                                     <  1:7615424/101‑1 (MQ=255)

                             tAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCAc                                                                         >  1:7616066/1‑96 (MQ=255)

                                 cAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTcg                                                                >  1:7616071/1‑101 (MQ=255)

                                  aaaTTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTcgc                                                               >  1:7616073/1‑101 (MQ=255)

                                   aaTTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGaa                                                                                    <  2:7615518/79‑1 (MQ=255)

                                   aaTTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTcgcg                                                              <  2:7615537/101‑1 (MQ=255)

                                      tCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGata                                                           >  1:7616077/1‑101 (MQ=255)

                                      tCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGCGTGCGCACCGAGCGTCGCGata                                                           >  2:7616079/1‑101 (MQ=255)

                                       cTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGat                                                            <  2:7615489/99‑1 (MQ=255)

                                        tGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGAtatt                                                         <  2:7615611/101‑1 (MQ=255)

                                         gCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGAtatta                                                        <  1:7615640/101‑1 (MQ=255)

                                           cGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATattatt                                                      >  2:7616083/1‑101 (MQ=255)

                                            ggACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATAttatta                                                     <  1:7615199/101‑1 (MQ=255)

                                            ggACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATAttatta                                                     <  2:7615583/101‑1 (MQ=255)

                                            ggACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATAttatta                                                     >  1:7616086/1‑101 (MQ=255)

                                              aCTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTaaa                                                   >  2:7616092/1‑101 (MQ=255)

                                              aCTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTaaa                                                   >  1:7616090/1‑101 (MQ=255)

                                               cTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAg                                                  >  2:7616096/1‑101 (MQ=255)

                                                tGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGa                                                 <  2:7615769/101‑1 (MQ=255)

                                                tGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGa                                                 >  1:7616098/1‑101 (MQ=255)

                                                  gCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATg                                               <  1:7615683/101‑1 (MQ=255)

                                                   cTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATgg                                              <  1:7615659/101‑1 (MQ=255)

                                                    tGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGc                                             >  1:7616101/1‑101 (MQ=255)

                                                      aCCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCAt                                           <  2:7615906/101‑1 (MQ=255)

                                                       ccATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATc                                          <  2:7615445/101‑1 (MQ=255)

                                                       ccATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATc                                          >  2:7616105/1‑101 (MQ=255)

                                                        cATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATcc                                         >  1:7616106/1‑101 (MQ=255)

                                                          ttGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTg                                       <  2:7615521/101‑1 (MQ=255)

                                                          ttGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTg                                       <  2:7615718/101‑1 (MQ=255)

                                                          ttGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTg                                       >  1:7616110/1‑101 (MQ=255)

                                                          ttGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTg                                       >  2:7616111/1‑101 (MQ=255)

                                                           tGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGa                                      <  1:7615854/101‑1 (MQ=255)

                                                             aaGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACt                                    >  2:7616115/1‑101 (MQ=255)

                                                             aaGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACt                                    <  2:7615654/101‑1 (MQ=255)

                                                                agCATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAg                                 >  2:7616119/1‑101 (MQ=255)

                                                                  cATCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTg                               <  2:7615756/101‑1 (MQ=255)

                                                                   aTCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTgg                              <  1:7615615/101‑1 (MQ=255)

                                                                   aTCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTgg                              >  1:7616121/1‑101 (MQ=255)

                                                                   aTCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTgg                              >  2:7616124/1‑101 (MQ=255)

                                                                    tCCGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGgc                             >  2:7616128/1‑101 (MQ=255)

                                                                     ccGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGgct                            >  1:7616129/1‑101 (MQ=255)

                                                                      cGCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGgctg                           <  1:7615393/101‑1 (MQ=255)

                                                                       gCATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGgctgc                          <  1:7615589/101‑1 (MQ=255)

                                                                        cATCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGgctgct                         <  1:7615750/101‑1 (MQ=255)

                                                                         aTCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGgctgc                          >  2:7616134/1‑99 (MQ=255)

                                                                         aTCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGctgctg                        <  2:7615690/101‑1 (MQ=255)

                                                                         aTCTGCTGAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGctgctg                        <  1:7615974/101‑1 (MQ=255)

                                                                            tgctgAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGAcc                     <  2:7615632/101‑1 (MQ=255)

                                                                            tgctgAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGAcc                     >  1:7616136/1‑101 (MQ=255)

                                                                              ctgAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAg                   >  1:7616139/1‑101 (MQ=255)

                                                                              ctgAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACGCAGTGGCTGCTGACCAg                   <  1:7615646/101‑1 (MQ=255)

                                                                               tgAAAGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGc                  >  2:7616141/1‑101 (MQ=255)

                                                                                  aaGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTAc               <  2:7615825/101‑1 (MQ=255)

                                                                                  aaGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTAc               >  2:7616145/1‑101 (MQ=255)

                                                                                   aGGGCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACt              >  1:7616146/1‑101 (MQ=255)

                                                                                    gggCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTc             >  1:7616150/1‑101 (MQ=255)

                                                                                      gCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTc             <  1:7615765/99‑1 (MQ=255)

                                                                                      gCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCgg           <  2:7615808/101‑1 (MQ=255)

                                                                                      gCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCgg           >  2:7616153/1‑101 (MQ=255)

                                                                                      gCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCgg           <  1:7615902/101‑1 (MQ=255)

                                                                                      gCTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCgg           <  1:7615818/101‑1 (MQ=255)

                                                                                       cTGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCGgc          >  2:7616157/1‑101 (MQ=255)

                                                                                        tGGCGTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCGgcg         >  1:7616158/1‑101 (MQ=255)

                                                                                            gTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCGGCGCGGa     >  2:7616164/1‑101 (MQ=255)

                                                                                            gTCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCGGCGCGGa     >  1:7616163/1‑101 (MQ=255)

                                                                                             tCGACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCGGCGCgg      <  1:7615714/99‑1 (MQ=255)

                                                                                               gACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGAtatta                                                        <  1:7615439/47‑1 (MQ=255)

                                                                                               gACACCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCGGCGCGGAAAc  >  1:7616169/1‑101 (MQ=255)

AAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACGCCATTCGACAGCGAAGGTGTGCGCACCGAGCGTCGCGATATTATTAAAGATGGCATCCTGACTCAGTGGCTGCTGACCAGCTACTCGGCGCGGAAAC  >  NC_000913/4458772‑4458967

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