BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,998,893	A→G	L64P (CTG→CCG)	fliZ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,998,893	0	A	G	85.5%	544.2 / 86.4	235	L64P (CTG→CCG)	fliZ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (11/23); new base G (103/98); total (114/121)
Fisher's exact test for biased strand distribution p-value = 6.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

Only 100 of 239 total aligned reads displayed.

GCGTTTTGCAATGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATGGGTC  >  NC_000913/1998801‑1998982

  gTTTTGCAATGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCggg                                                                                         >  2:4015947/1‑93 (MQ=255)

   ttttGCAATGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGc                                                                                        >  2:5215416/1‑93 (MQ=255)

   ttttGCAATGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGc                                                                                        <  1:3760383/93‑1 (MQ=255)

   ttttGCAATGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGc                                                                                        >  1:4586390/1‑93 (MQ=255)

   ttttGCAATGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGc                                                                                        >  2:364631/1‑93 (MQ=255)

    tttGCAATGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCg                                                                                       <  2:5501571/93‑1 (MQ=255)

    tttGCAATGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCg                                                                                       >  2:7128970/1‑93 (MQ=255)

      tGCAATGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCcgggcggg                                                                                     <  2:6366363/93‑1 (MQ=255)

       gCAATGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGa                                                                                    <  1:1639803/93‑1 (MQ=255)

         aaTGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGaaa                                                                                  >  1:2709037/1‑93 (MQ=255)

         aaTGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGaaa                                                                                  >  1:514090/1‑93 (MQ=255)

           tGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAAtc                                                                                >  1:4417236/1‑93 (MQ=255)

             tAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAAtctc                                                                              <  2:3909596/93‑1 (MQ=255)

             tAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTCGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAAtctc                                                                              <  1:5273453/93‑1 (MQ=255)

               aaTCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGa                                                                            <  2:6411697/93‑1 (MQ=255)

               aaTCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGa                                                                            <  2:2705645/93‑1 (MQ=255)

               aaTCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGa                                                                            <  2:2634196/93‑1 (MQ=255)

                 tCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGAtt                                                                          <  1:1954206/93‑1 (MQ=255)

                  cACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGAttt                                                                         <  1:2590097/93‑1 (MQ=255)

                  cACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGAttt                                                                         >  1:6895553/1‑93 (MQ=255)

                   aCCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGAtttt                                                                        <  1:124687/93‑1 (MQ=255)

                   aCCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGAtttt                                                                        >  2:1960006/1‑93 (MQ=255)

                   aCCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGAtttt                                                                        >  2:6685573/1‑93 (MQ=255)

                   aCCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGCGTCCGGGCGGGAAATCTCGAtttt                                                                        >  1:2794533/1‑93 (MQ=255)

                    ccGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTa                                                                       >  2:7487981/1‑93 (MQ=255)

                    ccGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTa                                                                       <  1:6826116/93‑1 (MQ=255)

                     cGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTAt                                                                      >  2:4996692/1‑93 (MQ=255)

                      gCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTAtt                                                                     <  2:7268836/93‑1 (MQ=255)

                       cAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTc                                                                    >  1:1298707/1‑93 (MQ=255)

                       cAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTc                                                                    <  1:508596/93‑1 (MQ=255)

                        aaaaTCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCa                                                                   <  1:5537328/93‑1 (MQ=255)

                         aaaTCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCAc                                                                  >  2:3161657/1‑93 (MQ=255)

                           aTCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTa                                                                >  1:6300653/1‑93 (MQ=255)

                           aTCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTa                                                                <  2:617266/93‑1 (MQ=255)

                            tCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTAt                                                               <  2:2780936/93‑1 (MQ=255)

                            tCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTAt                                                               <  2:2768254/93‑1 (MQ=255)

                            tCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTAt                                                               <  1:7016229/93‑1 (MQ=255)

                            tCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTAt                                                               <  2:3144149/93‑1 (MQ=255)

                              ggCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTAttt                                                             >  1:6740961/1‑93 (MQ=255)

                              ggCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTAttt                                                             >  1:4284396/1‑93 (MQ=255)

                               gCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTAtttt                                                            <  2:1764666/93‑1 (MQ=255)

                               gCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTAtttt                                                            >  2:5416485/1‑93 (MQ=255)

                                    aTGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGt                                                       >  1:600960/1‑93 (MQ=255)

                                     tGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTc                                                      <  2:31650/93‑1 (MQ=255)

                                        cgccCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTcgcg                                                   <  2:6439007/93‑1 (MQ=255)

                                         gccCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGa                                                  >  2:371964/1‑93 (MQ=255)

                                          cccATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGaa                                                 <  2:1236273/93‑1 (MQ=255)

                                          cccATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGaa                                                 <  1:4203009/93‑1 (MQ=255)

                                            cATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACt                                               <  1:2009129/93‑1 (MQ=255)

                                            cATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACt                                               >  2:1355604/1‑93 (MQ=255)

                                             aTGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTa                                              >  1:841279/1‑93 (MQ=255)

                                               gATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAg                                            >  2:2112941/1‑93 (MQ=255)

                                                aTTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAgg                                           <  2:2999289/93‑1 (MQ=255)

                                                   tttGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTaa                                        <  1:289366/93‑1 (MQ=255)

                                                   tttGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTaa                                        <  2:1614671/93‑1 (MQ=255)

                                                   tttGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTaa                                        <  1:3917987/93‑1 (MQ=255)

                                                     tGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATg                                      <  2:57635/93‑1 (MQ=255)

                                                      gTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGc                                     >  1:1673418/1‑93 (MQ=255)

                                                        tCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGa                                   >  2:377404/1‑93 (MQ=255)

                                                         ccGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGAt                                  <  1:3457935/93‑1 (MQ=255)

                                                          cGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATc                                 >  2:4764620/1‑93 (MQ=255)

                                                          cGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATc                                 <  2:4059041/93‑1 (MQ=255)

                                                           gTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCg                                >  1:2545140/1‑93 (MQ=255)

                                                            tttGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGa                               <  1:104644/93‑1 (MQ=255)

                                                             ttGCCAGCCATTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAg                              >  1:3168089/1‑93 (MQ=255)

                                                               gccagccaTTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTa                            >  2:2947630/1‑93 (MQ=255)

                                                               gccagccaTTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTa                            <  1:5686634/93‑1 (MQ=255)

                                                                ccagccaTTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTaa                           <  2:3616515/93‑1 (MQ=255)

                                                                 cagccaTTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAAt                          <  1:6534424/93‑1 (MQ=255)

                                                                 cagccaTTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAAt                          <  2:735262/93‑1 (MQ=255)

                                                                  agccaTTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAAtt                         <  1:726064/93‑1 (MQ=255)

                                                                    ccaTTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTa                       >  1:6372049/1‑93 (MQ=255)

                                                                     caTTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTAc                      <  2:3509706/93‑1 (MQ=255)

                                                                      aTTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACg                     <  2:2530748/93‑1 (MQ=255)

                                                                      aTTTTCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACg                     <  1:1738318/93‑1 (MQ=255)

                                                                       ttttCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACgg                    <  1:7120910/93‑1 (MQ=255)

                                                                       ttttCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACgg                    <  2:3612916/93‑1 (MQ=255)

                                                                       ttttCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACgg                    >  1:4716277/1‑93 (MQ=255)

                                                                       ttttCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACgg                    <  1:971944/93‑1 (MQ=255)

                                                                        tttCATCAAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGAAATTTACGGc                   >  2:781049/1‑93 (MQ=255)

                                                                          tcatcaAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCaa                 >  2:4628645/1‑93 (MQ=255)

                                                                          tcatcaAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCaa                 <  2:3456232/93‑1 (MQ=255)

                                                                           catcaAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAAt                <  2:5284448/93‑1 (MQ=255)

                                                                             tcaAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATgc              <  1:3544264/93‑1 (MQ=255)

                                                                             tcaAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATgc              >  2:2057855/1‑93 (MQ=255)

                                                                             tcaAGCAGCGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATgc              <  1:2908191/93‑1 (MQ=255)

                                                                             tcaAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATgc              >  2:4894623/1‑93 (MQ=255)

                                                                               aaGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATgcgc            >  1:4838015/1‑93 (MQ=255)

                                                                                     cGTGTCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATg      >  1:3993676/1‑93 (MQ=255)

                                                                                      gtgtCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATgg     >  1:7317608/1‑93 (MQ=255)

                                                                                       tgtCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATggg    >  2:2217723/1‑93 (MQ=255)

                                                                                       tgtCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATggg    <  1:2351689/93‑1 (MQ=255)

                                                                                       tgtCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATggg    >  1:4496953/1‑93 (MQ=255)

                                                                                       tgtCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATggg    <  1:6183620/93‑1 (MQ=255)

                                                                                       tgtCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATggg    <  1:6561356/93‑1 (MQ=255)

                                                                                        gtCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATGGGt   >  1:4051698/1‑93 (MQ=255)

                                                                                         tCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATGGGTc  >  1:4653739/1‑93 (MQ=255)

                                                                                         tCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATGGGTc  >  2:5592913/1‑93 (MQ=255)

                                                                                         tCCGGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATGGGTc  <  1:7463972/93‑1 (MQ=255)

                                                                                         tCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATGGGTc  <  1:3635806/93‑1 (MQ=255)

GCGTTTTGCAATGTAAATCACCGCAAAATCGGCACAATGCCGCCCATGATTTTTGTTCCGTTTGCCAGCCATTTTCATCAAGCAGCGTGTCCAGGCGGGAAATCTCGATTTTATTCACTATTTTGCCGTCGCGAACTAAGGTAATGCGATCGAGTAATTTACGGCAATGCGCGCAATGGGTC  >  NC_000913/1998801‑1998982

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