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	97.6%	392.6 / ‑0.2	127	Q133Q (CAA→CAG)	nlpB	lipoprotein‑34
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (59/65); total (61/66)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

Only 100 of 133 total aligned reads displayed.

GATTTGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGTTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAGC  >  NC_000913/2596397‑2596576

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

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

gaTTTGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTg                                                                                         <  2:3810230/93‑1 (MQ=255)

   ttGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCCTCACGCTGGGt                                                                                      >  1:3680853/1‑93 (MQ=255)

   ttGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTgggg                                                                                      >  1:1797995/1‑92 (MQ=255)

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

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

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

      aTAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGAt                                                                                   <  1:1231715/93‑1 (MQ=255)

      aTAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGAt                                                                                   <  1:2981877/93‑1 (MQ=255)

         aCGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTTATGgt                                                                                >  1:4002841/1‑93 (MQ=255)

         aCGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGgt                                                                                <  1:1775405/93‑1 (MQ=255)

           gACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGgtgt                                                                              >  1:420121/1‑93 (MQ=255)

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

           gACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAATGTCTGACCAGCATCATCACGCTGGGTGATGgtgt                                                                              <  1:2549413/93‑1 (MQ=255)

            gccACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTa                                                                             <  1:3857489/92‑1 (MQ=255)

            aCCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTa                                                                             <  1:3857490/93‑1 (MQ=255)

            aCCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTa                                                                             <  1:3170983/93‑1 (MQ=255)

            aCCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTa                                                                             <  1:933232/93‑1 (MQ=255)

            aCCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTa                                                                             <  1:3170982/93‑1 (MQ=255)

            aCCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTa                                                                             <  2:3157850/93‑1 (MQ=255)

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

             ccACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAg                                                                            >  2:1758410/1‑93 (MQ=255)

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

               aCGATACTGCTCGTCTTCGTCCAGACGGTTCCATTTTACCCACTCGGTGGTCAGTGTCTGACCCGCATCATCACGCTGGGTGATGGGGTAGtt                                                                          >  2:3844712/1‑93 (MQ=255)

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

                  aTACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGttttt                                                                       <  1:4726203/93‑1 (MQ=255)

                   tACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTc                                                                      >  2:2827201/1‑93 (MQ=255)

                     cTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGc                                                                    >  1:3670001/1‑93 (MQ=255)

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

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

                      tGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGcc                                                                   <  1:2706899/93‑1 (MQ=255)

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

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

                         tCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTgc                                                                <  2:2828099/93‑1 (MQ=255)

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

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

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

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

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

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

                                        cGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAaccacc                                                 <  1:4738019/93‑1 (MQ=255)

                                        cGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAaccacc                                                 <  1:1266453/93‑1 (MQ=255)

                                        cGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTCTCGCCTGCAGCACGCTAaccacc                                                 >  2:3392426/1‑93 (MQ=255)

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

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

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

                                           ttCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGc                                              <  1:2774568/93‑1 (MQ=255)

                                           ttCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGc                                              >  1:1117510/1‑93 (MQ=255)

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

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

                                                ttGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCca                                         >  2:4626201/1‑93 (MQ=255)

                                                   tACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACag                                      >  1:2907866/1‑93 (MQ=255)

                                                   tACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGCGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACag                                      >  2:1371874/1‑93 (MQ=255)

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

                                                     cccAATCGGTGGTTAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACagag                                    <  2:3429441/93‑1 (MQ=255)

                                                      ccAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGt                                   <  2:3062802/93‑1 (MQ=255)

                                                      ccAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGACGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGt                                   >  1:2045717/1‑93 (MQ=255)

                                                       cAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCGCCTGCGGCCACAGAGTa                                  <  2:830967/93‑1 (MQ=255)

                                                       cAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTa                                  <  1:1107056/93‑1 (MQ=255)

                                                        aaTCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTAt                                 <  2:513235/93‑1 (MQ=255)

                                                        aaTCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTAt                                 <  1:3632820/93‑1 (MQ=255)

                                                         aTCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCCCagaggctt                                >  2:22092/1‑89 (MQ=255)

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

                                                         aTCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTAtt                                >  2:22093/1‑93 (MQ=255)

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

                                                          tCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTg                               <  2:1849600/93‑1 (MQ=255)

                                                           cGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGc                              <  2:2790426/93‑1 (MQ=255)

                                                            ggtggtCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGcc                             <  1:2912597/93‑1 (MQ=255)

                                                            ggtggtCAGTGTCTGACCAGCATAATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGcc                             <  2:3436241/93‑1 (MQ=255)

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

                                                                 tCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGAc                        >  1:1541063/1‑93 (MQ=255)

                                                                 tCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGAc                        >  2:3529762/1‑93 (MQ=255)

                                                                   aGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCa                      <  1:3675776/93‑1 (MQ=255)

                                                                   aGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCa                      <  1:2237018/93‑1 (MQ=255)

                                                                   aGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCa                      <  1:458623/93‑1 (MQ=255)

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

                                                                         tGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTcc                >  2:250299/1‑93 (MQ=255)

                                                                         tGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTcc                >  2:3019578/1‑93 (MQ=255)

                                                                          gACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTcca               >  1:4730628/1‑93 (MQ=255)

                                                                              aGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAg           <  1:1407854/93‑1 (MQ=255)

                                                                               gCATCATCACGTTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGc          >  1:2086062/1‑93 (MQ=255)

                                                                               gCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGc          <  1:4509162/93‑1 (MQ=255)

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

                                                                                catcatCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCa         >  2:3674514/1‑93 (MQ=255)

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

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

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

                                                                                  tcatcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAAt       <  2:2468844/93‑1 (MQ=255)

                                                                                  tcatcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAAt       >  2:1006117/1‑93 (MQ=255)

                                                                                  tcatcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAAt       <  2:2435404/93‑1 (MQ=255)

                                                                                    atcaCGTTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGa     >  1:1027889/1‑93 (MQ=255)

                                                                                    atcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGa     >  1:1704210/1‑93 (MQ=255)

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

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

                                                                                     tcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGaa    >  2:2018336/1‑93 (MQ=255)

                                                                                      caCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAg   <  2:4091713/93‑1 (MQ=255)

                                                                                       aCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAGc  >  1:1613193/1‑93 (MQ=255)

                                                                                       aCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAGc  >  2:3542195/1‑93 (MQ=255)

                                                                                       aCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAGc  <  1:869284/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` `< 30 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`