BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	991,040	A→G	N399S (AAC→AGC)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	991,040	0	A	G	99.6%	816.1 / ‑5.0	254	N399S (AAC→AGC)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (112/141); total (112/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

Only 100 of 261 total aligned reads displayed.

CATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGT  >  NC_000913/990951‑991129

cATTGAGATGAACAACTTCTACGCCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCtt                                                                                        <  1:1495208/93‑1 (MQ=255)

 aTTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTc                                                                                       >  2:3505914/1‑93 (MQ=255)

 aTTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTc                                                                                       <  1:5319669/93‑1 (MQ=255)

 aTTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCACGATCCACACCCTGCTTGGCGAAGAAAGCTTc                                                                                       >  1:4831251/1‑93 (MQ=255)

  ttGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTcc                                                                                      <  2:4006276/93‑1 (MQ=255)

   tGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGCTTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCa                                                                                     <  2:5133747/93‑1 (MQ=255)

   tGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCa                                                                                     <  1:3813706/93‑1 (MQ=255)

   tGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCa                                                                                     >  1:737332/1‑93 (MQ=255)

    gagaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAg                                                                                    <  2:713252/93‑1 (MQ=255)

     agaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGa                                                                                   <  1:4674482/93‑1 (MQ=255)

       aTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGaaa                                                                                 >  2:18019/1‑93 (MQ=255)

       aTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGaaa                                                                                 >  2:6085757/1‑93 (MQ=255)

        tGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAg                                                                                <  1:1333049/93‑1 (MQ=255)

        tGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAg                                                                                >  1:335798/1‑93 (MQ=255)

         gAACAACTTCTGCACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAgg                                                                               <  1:2760097/93‑1 (MQ=255)

         gAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAgg                                                                               <  2:3599188/93‑1 (MQ=255)

          aacaacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAggg                                                                              <  2:5916445/93‑1 (MQ=255)

          aacaacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAggg                                                                              >  1:3049612/1‑93 (MQ=255)

           acaacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGa                                                                             >  2:1701346/1‑93 (MQ=255)

             aacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATg                                                                           >  1:6131199/1‑93 (MQ=255)

             aacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATg                                                                           >  1:1015067/1‑93 (MQ=255)

              acTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTTGCATGATACACACCCTGCTTGGCGACGAAAGCTTCCAGAAAGGGATGc                                                                          >  2:5660072/1‑93 (MQ=255)

              acTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGc                                                                          <  1:887525/93‑1 (MQ=255)

              acTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGc                                                                          <  2:2935884/93‑1 (MQ=255)

              acTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGc                                                                          >  2:1169731/1‑93 (MQ=255)

                 tCTACGCCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGc                                                                       <  2:2877262/93‑1 (MQ=255)

                 tCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGc                                                                       <  2:3462880/93‑1 (MQ=255)

                 tCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGc                                                                       <  1:5542578/93‑1 (MQ=255)

                   tACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCtt                                                                     >  1:1973701/1‑93 (MQ=255)

                      acCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTAt                                                                  <  2:4934768/93‑1 (MQ=255)

                      acCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTAt                                                                  >  1:311684/1‑93 (MQ=255)

                       cccTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTAtt                                                                 >  2:2180991/1‑93 (MQ=255)

                       cccTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTAtt                                                                 >  2:1585198/1‑93 (MQ=255)

                         cTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTc                                                               >  1:4478496/1‑93 (MQ=255)

                         cTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTc                                                               <  2:5260493/93‑1 (MQ=255)

                         cTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTc                                                               <  2:3662900/93‑1 (MQ=255)

                          tGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCAGACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCg                                                              <  2:6213859/93‑1 (MQ=255)

                          tGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCg                                                              >  2:4575282/1‑93 (MQ=255)

                           gACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGa                                                             >  1:3075621/1‑93 (MQ=255)

                           gACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGa                                                             <  2:5818432/93‑1 (MQ=255)

                            aCCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAg                                                            >  1:5288031/1‑93 (MQ=255)

                            aCCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAg                                                            <  1:99698/93‑1 (MQ=255)

                            aCCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAg                                                            >  1:1131957/1‑93 (MQ=255)

                             ccGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGc                                                           <  2:1222061/93‑1 (MQ=255)

                             ccGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGc                                                           >  1:649470/1‑93 (MQ=255)

                              cGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCATCTTTATTTCGAGCg                                                          <  2:1431122/93‑1 (MQ=255)

                              cGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCg                                                          >  1:5214518/1‑93 (MQ=255)

                                tttACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTc                                                        >  2:2760596/1‑93 (MQ=255)

                                   aCGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCatg                                                     >  2:1152960/1‑93 (MQ=255)

                                     gagaAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCatgat                                                   <  1:4524108/93‑1 (MQ=255)

                                      agaAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCatgatg                                                  >  1:1675523/1‑93 (MQ=255)

                                       gaAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATgg                                                 >  1:4955092/1‑93 (MQ=255)

                                       gaAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATgg                                                 >  2:1994030/1‑93 (MQ=255)

                                       gaAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATgg                                                 >  1:181050/1‑93 (MQ=255)

                                        aaGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGt                                                <  1:6045619/93‑1 (MQ=255)

                                         aGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTa                                               <  2:2079533/93‑1 (MQ=255)

                                         aGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTa                                               >  1:4305758/1‑93 (MQ=255)

                                         aGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTa                                               <  2:4829429/93‑1 (MQ=255)

                                               cgGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGccg                                         >  2:2709582/1‑91 (MQ=255)

                                                ggAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGc                                        >  1:165629/1‑93 (MQ=255)

                                                ggAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGc                                        <  2:3735774/93‑1 (MQ=255)

                                                  aaGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGa                                      <  2:5216855/93‑1 (MQ=255)

                                                  aaGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGa                                      <  2:5321312/93‑1 (MQ=255)

                                                   aGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGAc                                     >  2:1627642/1‑93 (MQ=255)

                                                      gATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCtg                                  >  1:5579329/1‑93 (MQ=255)

                                                       aTTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCtgt                                 >  1:569706/1‑93 (MQ=255)

                                                         tCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTga                               <  2:3043389/93‑1 (MQ=255)

                                                         tCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTga                               <  1:3939122/93‑1 (MQ=255)

                                                         tCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTga                               <  1:3843287/93‑1 (MQ=255)

                                                          cGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgac                              <  2:4202890/93‑1 (MQ=255)

                                                            cATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacga                            <  1:5198970/93‑1 (MQ=255)

                                                             aTGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacgcc                           >  2:435785/1‑91 (MQ=255)

                                                             aTGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacgac                           >  2:1613456/1‑93 (MQ=255)

                                                             aTGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacgac                           <  2:423768/93‑1 (MQ=255)

                                                               gATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACtt                         <  2:3148251/93‑1 (MQ=255)

                                                                aTCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACttt                        >  1:2673120/1‑93 (MQ=255)

                                                                aTCCACACCCTGCTGGGCGAAGAAAGCTTCCAGAAAGAGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACttt                        <  2:2801550/93‑1 (MQ=255)

                                                                 tCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTtg                       >  1:3157964/1‑93 (MQ=255)

                                                                  ccACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTtgt                      <  1:1039807/93‑1 (MQ=255)

                                                                  ccACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTtgt                      <  1:2378259/93‑1 (MQ=255)

                                                                   cacaCCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTtgtg                     >  2:5417599/1‑93 (MQ=255)

                                                                    acacCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGc                    <  2:737332/93‑1 (MQ=255)

                                                                    acacCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGc                    <  2:2859688/93‑1 (MQ=255)

                                                                     cacCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCa                   >  1:4045182/1‑93 (MQ=255)

                                                                     cacCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCa                   >  1:5867607/1‑93 (MQ=255)

                                                                      acCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAg                  <  2:4065201/93‑1 (MQ=255)

                                                                      acCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAg                  >  1:3460624/1‑93 (MQ=255)

                                                                      acCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAg                  >  2:2529611/1‑93 (MQ=255)

                                                                         cTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCg               >  1:11838/1‑93 (MQ=255)

                                                                          tGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGa              <  1:3555199/93‑1 (MQ=255)

                                                                           gCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGAt             <  2:4809481/93‑1 (MQ=255)

                                                                               ggCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGaa         <  2:5474645/93‑1 (MQ=255)

                                                                               ggCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGaa         >  2:5576441/1‑93 (MQ=255)

                                                                               ggCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGaa         <  1:2961004/93‑1 (MQ=255)

                                                                               ggCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGaa         <  2:466128/93‑1 (MQ=255)

                                                                                 cGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGa       <  1:669641/93‑1 (MQ=255)

                                                                                 cGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGa       <  2:4502142/93‑1 (MQ=255)

                                                                                   aagaaAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATg     >  1:2276476/1‑93 (MQ=255)

                                                                                    agaaAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGc    >  2:5087726/1‑93 (MQ=255)

                                                                                      aaagCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGt  >  1:324628/1‑93 (MQ=255)

CATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGT  >  NC_000913/990951‑991129

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