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%	822.1 / ‑5.4	255	Q133Q (CAA→CAG)	nlpB	lipoprotein‑34
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (108/146); total (109/146)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 258 total aligned reads displayed.

GATTTGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGTTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAGC  >  NC_000913/2596397‑2596576

 aTTTGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTgg                                                                                        >  1:5661284/1‑93 (MQ=255)

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

 aTTTGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTgg                                                                                        >  1:3411183/1‑93 (MQ=255)

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

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

   ttGAGAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGt                                                                                      <  2:3093868/93‑1 (MQ=255)

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

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

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

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

       tAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATg                                                                                  >  1:1528441/1‑93 (MQ=255)

        aaCGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATgg                                                                                 >  2:7883267/1‑93 (MQ=255)

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

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

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

            aCCACGATACTGCTCGTCTTCGTCCAGACGGTTCCAATGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTa                                                                             <  1:4362487/93‑1 (MQ=255)

                 gATACTGCTCGTCTTCGTCCAGACGGTTCCATGGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGtttt                                                                        <  1:820191/93‑1 (MQ=255)

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

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

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

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

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

                          cGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTgca                                                               >  1:1526372/1‑93 (MQ=255)

                          cGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTgca                                                               >  1:9203980/1‑93 (MQ=255)

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

                             cTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTgcagca                                                            <  2:6903853/93‑1 (MQ=255)

                              ctCGTCCAGACGGTTCCGTTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCAc                                                           <  2:4580601/92‑1 (MQ=255)

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

                                cGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTGTTTCGCCTGCAGCACGc                                                         <  1:4053241/93‑1 (MQ=255)

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

                                  tCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTa                                                       >  2:2708947/1‑93 (MQ=255)

                                    cAGACGGTTCCATTTTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAac                                                     <  2:5620202/93‑1 (MQ=255)

                                    cAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAac                                                     <  1:5544218/93‑1 (MQ=255)

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

                                       actgTTCCCTTTTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAaccac                                                  <  1:3189941/90‑1 (MQ=255)

                                       aCGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAaccac                                                  <  2:5173980/93‑1 (MQ=255)

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

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

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

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

                                            tCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCg                                             <  1:669332/93‑1 (MQ=255)

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

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

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

                                                ttGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCca                                         <  2:5641485/93‑1 (MQ=255)

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

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

                                                     cccAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACagag                                    <  1:9983770/93‑1 (MQ=255)

                                                     cccAATCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACagag                                    <  1:2838375/93‑1 (MQ=255)

                                                     cccAATCGGTGGGCAGTGTCTGACCAGCATCATCACGCTGGATGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACagag                                    <  2:4668552/93‑1 (MQ=255)

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

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

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

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

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

                                                        aaTCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTAt                                 >  1:4661475/1‑93 (MQ=255)

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

                                                         aTCGGTGGTCAGTGTCTGACCAGCATCATCACGCTGGGGGATGGGGTAGTTTTTCGCCTGCACTACGCTAACCACCTCCGCACACACAGTAtc                                >  2:75794/1‑92 (MQ=255)

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

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

                                                            ggtggtCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGcc                             >  1:7788210/1‑93 (MQ=255)

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

                                                              tgggCAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCAc                           <  2:4314419/93‑1 (MQ=255)

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

                                                                 gcAGTGTCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGAc                        <  2:2219080/92‑1 (MQ=255)

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

                                                                    gtgtCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCAt                     >  1:8092491/1‑93 (MQ=255)

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

                                                                     tgtCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCAtt                    >  1:4093473/1‑93 (MQ=255)

                                                                     tgtCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCAtt                    >  2:877167/1‑93 (MQ=255)

                                                                     tgtCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCAtt                    >  1:6299255/1‑93 (MQ=255)

                                                                     tgtCTGACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCAtt                    <  1:7480134/93‑1 (MQ=255)

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

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

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

                                                                          gatcAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGACCACAGAGTATTGCCACGACCATTTTcca               <  1:6757506/90‑1 (MQ=255)

                                                                          gACCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCTCCTGCGGCCACAGAGTATTGCCACGACCATTTTcca               <  2:8446336/93‑1 (MQ=255)

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

                                                                           aCCAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTccac              <  1:8342833/93‑1 (MQ=255)

                                                                             cAGCATCATCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTccacca            <  2:4096617/93‑1 (MQ=255)

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

                                                                               gCATCGTCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGc          <  1:7784692/93‑1 (MQ=255)

                                                                                catcatCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCa         <  2:5432568/93‑1 (MQ=255)

                                                                                catcatCACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCa         <  2:6509815/93‑1 (MQ=255)

                                                                                 atcatcACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCaa        <  2:8018534/93‑1 (MQ=255)

                                                                                 atcatcACGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCaa        <  2:1363012/93‑1 (MQ=255)

                                                                                 atcaccACGCTGGGTGATGGTGTAGTTTTCCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCaa        <  1:8843728/93‑1 (MQ=255)

                                                                                  tcatcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCGCAGAGTATTGCCACGACCATTTTCCACCAGCAAt       >  2:6745581/1‑93 (MQ=255)

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

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

                                                                                   catcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATg      <  2:7431662/93‑1 (MQ=255)

                                                                                   catcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATg      <  2:1142691/93‑1 (MQ=255)

                                                                                    atcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCGCCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGa     <  1:5928819/93‑1 (MQ=255)

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

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

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

                                                                                     tcaCGCTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGaa    <  2:5936185/93‑1 (MQ=255)

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

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

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

GATTTGATAACGACCACGATACTGCTCGTCTTCGTCCAGACGGTTCCATTGTACCCAATCGGTGGTCAGTGTCTGACCAGCATCATCACGTTGGGTGATGGTGTAGTTTTTCGCCTGCAGCACGCTAACCACCTGCGGCCACAGAGTATTGCCACGACCATTTTCCACCAGCAATGAAGC  >  NC_000913/2596397‑2596576

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