BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,456,991	A→G	intergenic (+106/‑91)	pmbA → / → b4236	PmbA/cytochrome b(562)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,456,991	0	A	G	81.1%	640.8 / 147.5	307	intergenic (+106/‑91)	pmbA/b4236	PmbA/cytochrome b(562)
Reads supporting (aligned to +/- strand): ref base A (20/38); new base G (124/125); total (144/163)
Fisher's exact test for biased strand distribution p-value = 4.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 307 total aligned reads displayed.

CTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTGCAGTCTCCTC  >  NC_000913/4456901‑4457080

cTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCg                                                                                         <  1:3164577/93‑1 (MQ=255)

  gCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTg                                                                                       >  1:2563468/1‑93 (MQ=255)

   ccGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGa                                                                                      >  2:4035929/1‑93 (MQ=255)

      cgcTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGAccc                                                                                   <  1:186576/93‑1 (MQ=255)

       gcTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCg                                                                                  >  1:1297539/1‑93 (MQ=255)

       gcTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCg                                                                                  >  2:3960794/1‑93 (MQ=255)

        cTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCgg                                                                                 >  2:2592253/1‑93 (MQ=255)

         tataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGt                                                                                >  2:1902187/1‑93 (MQ=255)

          ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTc                                                                               >  1:882181/1‑93 (MQ=255)

          ataAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTc                                                                               >  2:6009959/1‑93 (MQ=255)

           taAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCg                                                                              >  2:144904/1‑93 (MQ=255)

           taAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCg                                                                              >  2:1937496/1‑93 (MQ=255)

             aaCTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGaa                                                                            <  1:1038323/93‑1 (MQ=255)

              aCTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAAt                                                                           >  2:2095288/1‑93 (MQ=255)

              aCTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAAt                                                                           >  2:3777286/1‑93 (MQ=255)

               cTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATg                                                                          >  2:668616/1‑93 (MQ=255)

               cTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATg                                                                          >  2:2971364/1‑93 (MQ=255)

                tAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGa                                                                         <  2:3356110/93‑1 (MQ=255)

                tAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGa                                                                         >  1:3286338/1‑93 (MQ=255)

                 aaaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAg                                                                        <  1:1478972/93‑1 (MQ=255)

                  aaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGc                                                                       >  1:6568513/1‑93 (MQ=255)

                  aaaaaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGc                                                                       >  1:3007383/1‑93 (MQ=255)

                     aaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTa                                                                    <  1:7307333/93‑1 (MQ=255)

                     aaTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTa                                                                    <  2:437535/93‑1 (MQ=255)

                      aTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAg                                                                   >  1:4836829/1‑93 (MQ=255)

                      aTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAg                                                                   <  2:5388706/93‑1 (MQ=255)

                      aTTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAg                                                                   >  1:4020531/1‑93 (MQ=255)

                       ttACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAgt                                                                  >  1:4540148/1‑93 (MQ=255)

                       ttACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAgt                                                                  >  2:6289287/1‑93 (MQ=255)

                       ttACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAgt                                                                  >  2:1535569/1‑93 (MQ=255)

                         acacAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAgtcg                                                                <  1:4209360/93‑1 (MQ=255)

                          cacaAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAgtcgt                                                               >  1:2304743/1‑93 (MQ=255)

                            caAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGtcgtcg                                                             >  1:797385/1‑93 (MQ=255)

                              aaaTCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGAt                                                           >  1:3441306/1‑93 (MQ=255)

                               aaTCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATg                                                          <  2:660965/93‑1 (MQ=255)

                                aTCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGa                                                         >  1:2779358/1‑93 (MQ=255)

                                   aTTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGc                                                      <  1:1741496/93‑1 (MQ=255)

                                   aTTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGc                                                      <  1:2154347/93‑1 (MQ=255)

                                   aTTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGc                                                      >  1:5758588/1‑93 (MQ=255)

                                    ttCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCa                                                     <  1:1139264/93‑1 (MQ=255)

                                    ttCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCa                                                     >  1:636689/1‑93 (MQ=255)

                                     tCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAg                                                    <  1:987818/93‑1 (MQ=255)

                                       gCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAgtg                                                  <  2:3210644/93‑1 (MQ=255)

                                       gCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAgtg                                                  >  1:2457683/1‑93 (MQ=255)

                                        cACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAgtgt                                                 <  2:3304810/93‑1 (MQ=255)

                                          ntGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTNa                                               >  2:5210814/2‑93 (MQ=255)

                                           tGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGaa                                              >  1:3878240/1‑93 (MQ=255)

                                           tGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGaa                                              <  1:4670439/93‑1 (MQ=255)

                                           tGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGaa                                              >  2:1685652/1‑93 (MQ=255)

                                           tGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGaa                                              >  1:5334473/1‑93 (MQ=255)

                                            gCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAg                                             <  2:3918510/93‑1 (MQ=255)

                                            gCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAg                                             >  2:5887603/1‑93 (MQ=255)

                                             cATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAgg                                            >  2:5692361/1‑93 (MQ=255)

                                              aTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGa                                           <  2:2864400/93‑1 (MQ=255)

                                              aTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGa                                           >  2:252510/1‑93 (MQ=255)

                                              aTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGa                                           >  2:6147136/1‑93 (MQ=255)

                                              aTCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGa                                           >  2:1975423/1‑93 (MQ=255)

                                               tCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGAt                                          <  1:3534180/93‑1 (MQ=255)

                                               tCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGAt                                          <  2:494280/93‑1 (MQ=255)

                                               tCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGAt                                          <  1:4804603/93‑1 (MQ=255)

                                               tCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGAt                                          <  1:4136095/93‑1 (MQ=255)

                                                 gagaCGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGCTGa                                        >  2:2394739/1‑93 (MQ=255)

                                                  agaCGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGaa                                       >  2:4743662/1‑93 (MQ=255)

                                                   gaCGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAg                                      <  2:769012/93‑1 (MQ=255)

                                                     cGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAgtg                                    <  2:6049088/93‑1 (MQ=255)

                                                       gAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGta                                  <  1:5156269/93‑1 (MQ=255)

                                                        aaaGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGtaa                                 >  2:7143489/1‑93 (MQ=255)

                                                         aaGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGtaaa                                <  1:5051445/93‑1 (MQ=255)

                                                         aaGTGAATGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGtaaa                                <  2:6434028/93‑1 (MQ=255)

                                                         aaGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGtaaa                                >  1:2966024/1‑93 (MQ=255)

                                                              aaTGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATaaaa                           <  1:5977817/93‑1 (MQ=255)

                                                               aTGAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATaaaaa                          <  1:3760249/93‑1 (MQ=255)

                                                               aTGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATaaaaa                          <  1:779657/93‑1 (MQ=255)

                                                                 gAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAg                        <  1:1823146/93‑1 (MQ=255)

                                                                 gAGGGCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAg                        >  1:7124673/1‑93 (MQ=255)

                                                                   gggCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGa                      <  1:6882066/93‑1 (MQ=255)

                                                                   gggCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGa                      >  2:2234106/1‑93 (MQ=255)

                                                                   gggCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGa                      <  1:5268706/93‑1 (MQ=255)

                                                                    ggCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGaa                     <  2:5848909/93‑1 (MQ=255)

                                                                    ggCCGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGaa                     <  2:5136693/93‑1 (MQ=255)

                                                                      ccGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACt                   <  1:1580773/93‑1 (MQ=255)

                                                                       cGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTa                  >  2:4846399/1‑93 (MQ=255)

                                                                       cGGGGAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTa                  <  2:1723680/93‑1 (MQ=255)

                                                                        ggggAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTAt                 >  2:5102752/1‑93 (MQ=255)

                                                                        ggggAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTAt                 <  1:4274086/93‑1 (MQ=255)

                                                                        ggggAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTAt                 <  2:3520626/93‑1 (MQ=255)

                                                                          ggAGCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTc               <  2:1018026/93‑1 (MQ=255)

                                                                          ggAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTc               >  1:2001487/1‑93 (MQ=255)

                                                                           gAGCGTACTCGCAGTGCGTGACCCGGTCAAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCt              <  2:424601/93‑1 (MQ=255)

                                                                             gCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTg            <  1:5534435/93‑1 (MQ=255)

                                                                             gCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTg            >  1:6809528/1‑93 (MQ=255)

                                                                             gCGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTg            >  2:6818852/1‑93 (MQ=255)

                                                                              cGTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTGc           <  1:7053296/93‑1 (MQ=255)

                                                                              cGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTGc           <  2:1126679/93‑1 (MQ=255)

                                                                               gTACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTGCa          >  1:6091004/1‑93 (MQ=255)

                                                                                tACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTGCAg         >  1:2701376/1‑93 (MQ=255)

                                                                                tACTCGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTGCAg         >  1:5636214/1‑93 (MQ=255)

                                                                                    cGCAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTGCAGtctc     >  1:6595394/1‑93 (MQ=255)

                                                                                      cAGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTGCAGTctcct   <  2:3629263/93‑1 (MQ=255)

                                                                                       aGTGCGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTGCAGTctcctc  >  2:5530356/1‑93 (MQ=255)

CTGCCGCGCTATAAACTAAAAAATTACACAAAATCATTCGCACTGCATCGAGACGGAAAGTGAATGAGGGCCGGGGAGCGTACTCGCAGTACGTGACCCGGTCGAATGAGCGTAGTCGTCGATGAGGCAGTGTGAAGGATGAAGTGTAAATAAAAAAGGAACTATTCTTGCAGTCTCCTC  >  NC_000913/4456901‑4457080

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