BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,596,487	T→C	Q133Q (CAA→CAG)	nlpB ←	lipoprotein‑34

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,596,487	0	T	C	99.6%	803.1 / NA	245	Q133Q (CAA→CAG)	nlpB	lipoprotein‑34
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (116/128); minor base G (0/1); total (116/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

GATTTGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGTTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAGC  >  NC_000913/2596397‑2596576

gaTTTGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTg                                                                                         >  1:1183682/1‑93 (MQ=255)

 aTTTGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTgg                                                                                        <  1:1559962/93‑1 (MQ=255)

  tttGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTggg                                                                                       >  2:2487291/1‑93 (MQ=255)

  tttGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTggg                                                                                       >  1:4644688/1‑93 (MQ=255)

  tttGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTggg                                                                                       <  2:4623904/93‑1 (MQ=255)

   ttGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGt                                                                                      <  1:5432751/93‑1 (MQ=255)

     gATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGa                                                                                    >  2:5280504/1‑93 (MQ=255)

     gATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGa                                                                                    <  1:3958506/93‑1 (MQ=255)

      aTAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGAt                                                                                   >  2:1029711/1‑93 (MQ=255)

      aTAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGAt                                                                                   >  2:400211/1‑93 (MQ=255)

       tAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATg                                                                                  <  2:209252/93‑1 (MQ=255)

         aCGACCACGATACTGCTCGTCTTCGTCCAGACGGGTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATggg                                                                                >  1:1016818/1‑92 (MQ=255)

           gACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGgtgt                                                                              <  1:3752741/93‑1 (MQ=255)

            aCCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTa                                                                             >  2:3727137/1‑93 (MQ=255)

            aCCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTa                                                                             >  1:2989664/1‑93 (MQ=255)

             ccACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAg                                                                            <  1:2382617/93‑1 (MQ=255)

             ccACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAg                                                                            <  2:354912/93‑1 (MQ=255)

             ccACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAg                                                                            <  1:4724359/93‑1 (MQ=255)

              cACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGt                                                                           <  2:5229298/93‑1 (MQ=255)

              cACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGt                                                                           >  2:1823162/1‑93 (MQ=255)

               aCGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGtt                                                                          >  1:525274/1‑93 (MQ=255)

               aCGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGtt                                                                          <  2:4757529/93‑1 (MQ=255)

               aCGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGtt                                                                          <  2:6249/93‑1 (MQ=255)

                cGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGttt                                                                         >  1:2733948/1‑93 (MQ=255)

                cGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGttt                                                                         >  1:6663435/1‑93 (MQ=255)

                cGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGttt                                                                         >  1:3624811/1‑93 (MQ=255)

                cGATACTGCTCGTCTTCGTCCAGAAGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGttt                                                                         >  2:4841232/1‑93 (MQ=255)

                 gATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGtttt                                                                        >  1:3692236/1‑93 (MQ=255)

                 gATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGtttt                                                                        <  1:2496500/93‑1 (MQ=255)

                  aTACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGttttt                                                                       >  1:1086922/1‑93 (MQ=255)

                   tACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTc                                                                      <  1:3037872/93‑1 (MQ=255)

                    aCTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGAGATGGTGTAGTTTTTCg                                                                     >  1:4224811/1‑93 (MQ=255)

                      tGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGcc                                                                   >  1:5380780/1‑93 (MQ=255)

                       gCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCt                                                                  >  1:6209798/1‑93 (MQ=255)

                       gCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCt                                                                  <  2:809143/93‑1 (MQ=255)

                        cTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTg                                                                 >  2:5141559/1‑93 (MQ=255)

                         tCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTgc                                                                >  2:2828874/1‑93 (MQ=255)

                          cGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTgca                                                               <  2:4368946/93‑1 (MQ=255)

                          cGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTgca                                                               <  1:2087020/93‑1 (MQ=255)

                           gTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTgcag                                                              <  1:903874/93‑1 (MQ=255)

                            tCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTgcagc                                                             >  2:5425833/1‑93 (MQ=255)

                              ttCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCAc                                                           <  2:1451445/93‑1 (MQ=255)

                              ttCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCAc                                                           >  2:1306198/1‑93 (MQ=255)

                              ttCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCAc                                                           <  1:6439321/93‑1 (MQ=255)

                                cGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGc                                                         >  2:2122991/1‑93 (MQ=255)

                                cGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGc                                                         <  1:5937498/93‑1 (MQ=255)

                                 gTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCt                                                        >  2:5381631/1‑93 (MQ=255)

                                 gTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCt                                                        <  2:310663/93‑1 (MQ=255)

                                  tCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTa                                                       <  1:576414/93‑1 (MQ=255)

                                   ccAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTaa                                                      <  2:3194155/93‑1 (MQ=255)

                                   ccAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTaa                                                      <  1:3254513/93‑1 (MQ=255)

                                    cAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAac                                                     >  2:5346056/1‑93 (MQ=255)

                                    cAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAac                                                     >  1:5152270/1‑93 (MQ=255)

                                     aGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAacc                                                    >  2:6443881/1‑93 (MQ=255)

                                     aGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAacc                                                    >  2:4312539/1‑93 (MQ=255)

                                     aGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAacc                                                    >  1:1212806/1‑93 (MQ=255)

                                      gACGGTTCCATTGTACCCAATCGGTGGTGAGTGTCTGACCAGCATCATCACGCTGGGGGATGGTGTAGGTTTTCGCCTGCAGCACGCTAacca                                                   >  1:4370701/1‑93 (MQ=255)

                                      gACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAacca                                                   <  2:3571346/93‑1 (MQ=255)

                                         ggTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCt                                                <  2:4975096/93‑1 (MQ=255)

                                         ggTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCt                                                >  1:2043500/1‑93 (MQ=255)

                                          gTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTg                                               <  1:253223/93‑1 (MQ=255)

                                          gTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTg                                               <  1:2638690/93‑1 (MQ=255)

                                          gTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTg                                               >  2:5683359/1‑93 (MQ=255)

                                           ttCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCATCTGc                                              >  2:4860595/1‑93 (MQ=255)

                                             ccATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCgg                                            >  2:5507525/1‑93 (MQ=255)

                                              cATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGc                                           >  2:5191813/1‑93 (MQ=255)

                                               aTTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGcc                                          >  2:1142458/1‑93 (MQ=255)

                                                 tGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCcac                                        <  1:4335351/93‑1 (MQ=255)

                                                  gTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCcaca                                       <  2:6217455/93‑1 (MQ=255)

                                                    aCCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACaga                                     <  1:1842074/93‑1 (MQ=255)

                                                    aCCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACaga                                     <  2:4558299/93‑1 (MQ=255)

                                                      ccAATCGGTGGACAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACagaga                                   <  1:6724658/93‑2 (MQ=255)

                                                       cAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTa                                  >  2:2721865/1‑93 (MQ=255)

                                                        aaTCGGTGGTTAGTGTCTTACCAGCATCATCACGCTGGGTGATTGTTTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTAt                                 <  1:5956639/93‑1 (MQ=255)

                                                         aTCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTAtt                                <  1:3406562/93‑1 (MQ=255)

                                                         aTCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTAtt                                <  2:4659987/93‑1 (MQ=255)

                                                             gtggtCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCa                            <  2:1790108/93‑1 (MQ=255)

                                                             gtggtCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCa                            >  1:782515/1‑93 (MQ=255)

                                                                gtCTGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGa                         <  1:3495602/93‑1 (MQ=255)

                                                                gtCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGa                         <  2:5541059/93‑1 (MQ=255)

                                                                gtCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGa                         <  2:3238899/93‑1 (MQ=255)

                                                                 tCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGAc                        <  1:1650079/93‑1 (MQ=255)

                                                                  cAGTGCCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGAcc                       <  1:1716939/93‑1 (MQ=255)

                                                                    gtgtCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCAt                     <  1:2918891/93‑1 (MQ=255)

                                                                       tCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCAtttt                  <  1:2128108/93‑1 (MQ=255)

                                                                       tCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCAtttt                  >  2:4940557/1‑93 (MQ=255)

                                                                        cTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTc                 <  1:320329/93‑1 (MQ=255)

                                                                          gACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTcca               <  2:1197731/93‑1 (MQ=255)

                                                                          gACCAGCAGCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTcca               <  2:2338327/93‑1 (MQ=255)

                                                                            ccAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTccacc             >  2:6476461/1‑93 (MQ=255)

                                                                            ccAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTccacc             >  1:5634035/1‑93 (MQ=255)

                                                                              aGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAg           >  2:5031453/1‑93 (MQ=255)

                                                                               gCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGc          >  2:3850776/1‑93 (MQ=255)

                                                                                 atcatcACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCaa        >  2:2532973/1‑93 (MQ=255)

                                                                                  tcatcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGAATTGCCACGACCATTTTCCACCAGCAAt       >  2:1913027/1‑93 (MQ=255)

                                                                                    atcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGa     <  1:1412653/93‑1 (MQ=255)

                                                                                    atcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGa     <  2:2293530/93‑1 (MQ=255)

                                                                                      caCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAg   >  2:6083405/1‑93 (MQ=255)

                                                                                      caCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAg   >  2:1980459/1‑93 (MQ=255)

                                                                                       aCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAGc  <  2:3395555/93‑1 (MQ=255)

GATTTGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGTTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAGC  >  NC_000913/2596397‑2596576

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