BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,951,082	0	C	T	66.5%	606.7 / inf	459	intergenic (+46/‑50)	ilvE/ilvD	branched‑chain amino‑acid aminotransferase/dihydroxyacid dehydratase
Reads supporting (aligned to +/- strand): ref base C (68/86); new base T (156/149); total (224/235)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 496 total aligned reads displayed.

CGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGGCC  >  NC_000913/3950991‑3951169

 gAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGatac                                                                                       <  1:4659995/93‑1 (MQ=255)

 gAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGatac                                                                                       <  1:10033444/93‑1 (MQ=255)

 gAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGacac                                                                                       <  1:5785005/93‑1 (MQ=255)

   aaCCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGGACGCACCGTCCCATTTACGAGACAGATACTg                                                                                     >  1:7446246/1‑93 (MQ=255)

     ccGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTggg                                                                                   <  2:2572010/93‑1 (MQ=255)

      cGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGa                                                                                  <  2:5030047/93‑1 (MQ=255)

       gAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAg                                                                                 <  2:1147053/93‑1 (MQ=255)

          gATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGtaa                                                                              >  1:4051760/1‑93 (MQ=255)

              aaTGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGtaaataa                                                                          <  2:1584156/93‑1 (MQ=255)

              aaTGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaaataa                                                                          >  1:9548413/1‑93 (MQ=255)

               aTGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGtaaataaa                                                                         <  2:2251009/93‑1 (MQ=255)

               aTGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaaataaa                                                                         >  2:5673564/1‑93 (MQ=255)

               aTGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACACACCGTCCCATTTACGAGACAGACACTGGGAGtaaataaa                                                                         <  1:3507186/93‑1 (MQ=255)

                tGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAg                                                                        >  2:2117943/1‑93 (MQ=255)

                tGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAg                                                                        >  2:3549843/1‑93 (MQ=255)

                  gggCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTa                                                                      <  1:5216866/93‑1 (MQ=255)

                  gggCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTa                                                                      <  1:1281590/93‑1 (MQ=255)

                    gCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATg                                                                    <  1:3826875/93‑1 (MQ=255)

                     cTGGTTAGATCAAGTTAATCAATAAATACAAAAAAGGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGc                                                                   <  2:772960/93‑1 (MQ=255)

                      tGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGcc                                                                  >  2:9700893/1‑93 (MQ=255)

                       ggTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGAAAATAAAGTATGCCt                                                                 >  1:5641860/1‑93 (MQ=255)

                       ggTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTCAATAAAGTATGCCt                                                                 >  1:3641757/1‑93 (MQ=255)

                        gTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTa                                                                >  1:5315532/1‑93 (MQ=255)

                        gTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTa                                                                >  1:5004085/1‑93 (MQ=255)

                        gTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTa                                                                >  2:3307783/1‑93 (MQ=255)

                        gTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTa                                                                >  1:2875872/1‑93 (MQ=255)

                         ttAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTaa                                                               >  2:1369940/1‑93 (MQ=255)

                         ttAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTaa                                                               <  2:519889/93‑1 (MQ=255)

                          tAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAg                                                              <  2:7549269/93‑1 (MQ=255)

                          tAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAg                                                              >  2:6216352/1‑93 (MQ=255)

                           aGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCATAAga                                                             >  1:2540094/1‑92 (MQ=255)

                            gATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGTCTAAGTa                                                            >  2:6666148/1‑93 (MQ=255)

                             aTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                           <  2:6874855/93‑1 (MQ=255)

                             aTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                           >  1:6433734/1‑93 (MQ=255)

                             aTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                           >  1:3411570/1‑93 (MQ=255)

                             aTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                           <  2:6932994/93‑1 (MQ=255)

                             aTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                           >  1:9785853/1‑93 (MQ=255)

                             aTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCACTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                           >  2:9456122/1‑93 (MQ=255)

                              tCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTAcc                                                          <  2:4565697/93‑1 (MQ=255)

                                aaGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGt                                                        >  2:733563/1‑93 (MQ=255)

                                aaGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGt                                                        >  1:8750820/1‑93 (MQ=255)

                                 aGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGtt                                                       >  2:8789642/1‑93 (MQ=255)

                                 aGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGtt                                                       <  1:1516717/93‑1 (MQ=255)

                                  gTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTc                                                      >  1:313768/1‑93 (MQ=255)

                                  gTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTc                                                      >  2:2440001/1‑93 (MQ=255)

                                  gTTAATCAATAAATACAAAAAATGGGACGGAACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTc                                                      >  1:6844487/1‑93 (MQ=255)

                                   ttAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTcc                                                     <  2:788804/93‑1 (MQ=255)

                                   ttAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTcc                                                     <  1:9082241/93‑1 (MQ=255)

                                   ttAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTcc                                                     <  1:6231467/93‑1 (MQ=255)

                                    tAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCg                                                    <  2:6549862/93‑1 (MQ=255)

                                    tAATCAATAAATACAAAAAATGGGAAGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGCAAATAAAGCATGACTAAGTACCGTTcag                                                    >  2:6419287/1‑91 (MQ=255)

                                      aTCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGcc                                                  >  1:9272553/1‑93 (MQ=255)

                                       tCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGcca                                                 >  2:4280546/1‑93 (MQ=255)

                                        cAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccac                                                <  1:8551749/93‑1 (MQ=255)

                                         aataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccatc                                               >  1:2292588/1‑91 (MQ=255)

                                         aataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccacc                                               >  1:4772709/1‑93 (MQ=255)

                                         aataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccacc                                               >  1:7016502/1‑93 (MQ=255)

                                         aataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTAAGCCTAAGGACCGTTCCGacacc                                               >  2:1971801/1‑93 (MQ=255)

                                          ataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccacca                                              >  2:5745934/1‑93 (MQ=255)

                                            aaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacc                                            >  1:3977894/1‑93 (MQ=255)

                                             aataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacca                                           >  2:1835168/1‑93 (MQ=255)

                                             aataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacca                                           >  2:9584668/1‑93 (MQ=255)

                                             aataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacca                                           >  1:6347423/1‑93 (MQ=255)

                                             aataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccaccacca                                           <  2:2801283/93‑1 (MQ=255)

                                                aCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTc                                        >  1:8100279/1‑93 (MQ=255)

                                                 cAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCa                                       <  1:7845468/93‑1 (MQ=255)

                                                 cAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCa                                       <  1:7072787/93‑1 (MQ=255)

                                                  aaaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCAt                                      >  2:3942186/1‑93 (MQ=255)

                                                  aaaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCAt                                      >  1:6058096/1‑93 (MQ=255)

                                                   aaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATg                                     >  2:8898457/1‑93 (MQ=255)

                                                      aaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTc                                  >  1:4504802/1‑93 (MQ=255)

                                                        tGGGACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGt                                >  1:8872482/1‑93 (MQ=255)

                                                           gACGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAat                             >  2:3745813/1‑93 (MQ=255)

                                                            gcGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAata                            <  1:2719195/92‑1 (MQ=255)

                                                             cGGCACGCACCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAatat                           <  1:7930861/93‑1 (MQ=255)

                                                             cGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAatat                           <  1:5527117/93‑1 (MQ=255)

                                                             cGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAatat                           <  1:757008/93‑1 (MQ=255)

                                                              ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATg                          <  1:234385/93‑1 (MQ=255)

                                                              ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATg                          <  1:9860071/93‑1 (MQ=255)

                                                              ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTACGCCACCACCACTCATGGTCGTAATATg                          <  2:3883948/93‑1 (MQ=255)

                                                                cacgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGc                        >  2:8339078/1‑93 (MQ=255)

                                                                 acgcacCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCg                       <  1:1653996/93‑1 (MQ=255)

                                                                  cgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTACGTACCGTTCTGTCACCACCACTCATGGTCGTAATATGGCgg                      >  2:4810972/1‑93 (MQ=255)

                                                                  cgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACACCAACCCATGGTGGAAATAGGGCgg                      >  1:8318098/1‑93 (MQ=255)

                                                                   gcacCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCggg                     >  1:8419681/1‑93 (MQ=255)

                                                                   gcacCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCggg                     <  1:204536/93‑1 (MQ=255)

                                                                    cacCGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCgggg                    >  1:7730366/1‑92 (MQ=255)

                                                                     acCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTg                   <  1:1847922/93‑1 (MQ=255)

                                                                      ccGTCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGc                  >  2:2521178/1‑93 (MQ=255)

                                                                         tCCCATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCg               <  1:4394484/93‑1 (MQ=255)

                                                                         tCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCg               <  2:8566422/93‑1 (MQ=255)

                                                                           ccATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTg             <  2:5921387/93‑1 (MQ=255)

                                                                            cATTTACGAGACAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTgc            <  2:1062235/93‑1 (MQ=255)

                                                                             aTTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTcgcgcg           >  2:8080714/1‑93 (MQ=255)

                                                                                 aCGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCtgt       <  2:7060211/93‑1 (MQ=255)

                                                                                 aCGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCtgt       <  2:5278082/93‑1 (MQ=255)

                                                                                 aCGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCtgt       <  2:9197650/93‑1 (MQ=255)

                                                                                   gagaCAGATACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGGTCGGGCGCTgggg     >  2:3102750/1‑93 (MQ=255)

                                                                                    agacagacACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTggc    <  1:4227680/93‑2 (MQ=255)

                                                                                      acagataCTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTggcgc  <  2:1771421/93‑4 (MQ=255)

CGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGGCC  >  NC_000913/3950991‑3951169

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