BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	477,096	A→G	L251P (CTC→CCC)	ylaB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	477,096	0	A	G	99.1%	709.1 / ‑4.7	221	L251P (CTC→CCC)	ylaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (98/121); total (99/122)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 227 total aligned reads displayed.

TAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGT  >  NC_000913/477007‑477185

tAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCg                                                                                        >  2:3675666/1‑93 (MQ=255)

tAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCg                                                                                        >  1:2665858/1‑93 (MQ=255)

 aGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCgg                                                                                       >  2:2917223/1‑93 (MQ=255)

  gCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCgga                                                                                      >  2:332516/1‑92 (MQ=255)

  gCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGt                                                                                      >  1:4623322/1‑93 (MQ=255)

   cATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTa                                                                                     >  2:5314143/1‑93 (MQ=255)

    ctCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTaa                                                                                    <  2:2737872/92‑1 (MQ=255)

    aTCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTaa                                                                                    >  2:5135285/1‑93 (MQ=255)

     tCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAc                                                                                   >  2:1315932/1‑93 (MQ=255)

     tCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAc                                                                                   >  1:2903630/1‑93 (MQ=255)

     tCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAc                                                                                   >  1:2770687/1‑93 (MQ=255)

     tCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAc                                                                                   >  1:2082367/1‑93 (MQ=255)

     tCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAc                                                                                   >  2:4074116/1‑93 (MQ=255)

         gCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCaaa                                                                               <  1:1081346/93‑1 (MQ=255)

          cAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAAt                                                                              <  1:4975944/93‑1 (MQ=255)

           aGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATa                                                                             <  1:1861937/93‑1 (MQ=255)

           aGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATa                                                                             <  2:781845/93‑1 (MQ=255)

            gCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATaa                                                                            >  1:5193343/1‑93 (MQ=255)

               gatgatGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGa                                                                         <  2:6038842/93‑1 (MQ=255)

               gatgatGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGa                                                                         <  1:1653362/93‑1 (MQ=255)

                atgatgCGGTGACTGAATGCGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGAc                                                                        <  1:1215512/93‑1 (MQ=255)

                 tgatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACt                                                                       <  2:6493953/93‑1 (MQ=255)

                  gatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACtt                                                                      >  1:645788/1‑93 (MQ=255)

                   atgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTg                                                                     <  2:6132299/93‑1 (MQ=255)

                   atgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTg                                                                     <  1:3871812/93‑1 (MQ=255)

                    tgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGc                                                                    <  2:382339/93‑1 (MQ=255)

                      cGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCg                                                                  >  2:710248/1‑93 (MQ=255)

                        gTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGAt                                                                >  1:5863490/1‑93 (MQ=255)

                          gACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGc                                                              <  2:101946/93‑1 (MQ=255)

                          gACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGc                                                              <  1:4023477/93‑1 (MQ=255)

                           aCTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGcc                                                             <  2:2789725/93‑1 (MQ=255)

                            cTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCa                                                            <  1:647315/93‑1 (MQ=255)

                             tGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAg                                                           >  1:579943/1‑93 (MQ=255)

                             tGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAg                                                           >  1:3424832/1‑93 (MQ=255)

                             tGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAg                                                           >  2:5378676/1‑93 (MQ=255)

                              gAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGc                                                          >  2:3576269/1‑93 (MQ=255)

                              gAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGc                                                          <  2:1153274/93‑1 (MQ=255)

                                aTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCt                                                        >  2:816457/1‑93 (MQ=255)

                                 tACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCtt                                                       <  1:679128/93‑1 (MQ=255)

                                  aCGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCttt                                                      >  1:5562260/1‑93 (MQ=255)

                                  aCGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCttt                                                      <  2:5049078/93‑1 (MQ=255)

                                  aCGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCttt                                                      >  1:1620782/1‑93 (MQ=255)

                                   cGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTc                                                     <  2:6426900/93‑1 (MQ=255)

                                   cGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTc                                                     <  1:5311898/93‑1 (MQ=255)

                                    ggCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCt                                                    <  1:2290803/93‑1 (MQ=255)

                                      cgcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGc                                                  <  2:5432921/93‑1 (MQ=255)

                                       gcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCa                                                 >  2:888237/1‑93 (MQ=255)

                                        cAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCaa                                                <  1:3996311/93‑1 (MQ=255)

                                         aaaaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAAc                                               <  2:4562300/93‑1 (MQ=255)

                                         aaaaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAAc                                               >  1:1182487/1‑93 (MQ=255)

                                          aaaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACa                                              >  2:6626158/1‑93 (MQ=255)

                                           aaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAt                                             <  1:6174372/93‑1 (MQ=255)

                                           aaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAt                                             <  1:5634059/93‑1 (MQ=255)

                                           aaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAt                                             <  2:4533399/93‑1 (MQ=255)

                                           aaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAt                                             >  1:4475726/1‑93 (MQ=255)

                                            aTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAtt                                            >  2:1992583/1‑93 (MQ=255)

                                            aTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAtt                                            <  2:4243421/93‑1 (MQ=255)

                                            aTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAtt                                            >  2:1863296/1‑93 (MQ=255)

                                             tACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAttt                                           <  1:1012581/93‑1 (MQ=255)

                                               cGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTc                                         <  2:4702842/93‑1 (MQ=255)

                                               cGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTc                                         <  2:3291345/93‑1 (MQ=255)

                                                  agcaCAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGtt                                      <  2:5964311/93‑1 (MQ=255)

                                                   gcaCAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTg                                     <  1:6259187/93‑1 (MQ=255)

                                                   gcaCAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTg                                     >  2:2582897/1‑93 (MQ=255)

                                                     acaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGaa                                   >  2:1373252/1‑93 (MQ=255)

                                                     acaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGaa                                   >  1:6197392/1‑93 (MQ=255)

                                                     acaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGaa                                   <  1:5391984/93‑1 (MQ=255)

                                                       aaaCATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGc                                 <  2:2311755/93‑1 (MQ=255)

                                                        aaCATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGcc                                >  1:5589390/1‑93 (MQ=255)

                                                          cATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCa                              <  1:824551/93‑1 (MQ=255)

                                                           aTCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAt                             <  1:5992398/93‑1 (MQ=255)

                                                            tCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtg                            >  1:5358165/1‑93 (MQ=255)

                                                             cGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtgt                           <  1:1714031/93‑1 (MQ=255)

                                                             cGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtgt                           <  2:2504928/93‑1 (MQ=255)

                                                             cGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtgt                           <  1:4690595/93‑1 (MQ=255)

                                                               cTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtga                         <  1:3623387/93‑1 (MQ=255)

                                                                 tccAGCAGGCCAATCACCAACCCGGGCCGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatg                       <  2:3170270/92‑1 (MQ=255)

                                                                 gCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatg                       <  2:5791269/93‑1 (MQ=255)

                                                                  ccAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatga                      <  1:3193860/93‑1 (MQ=255)

                                                                  ccAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatga                      <  1:5163890/93‑1 (MQ=255)

                                                                  ccAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatga                      <  2:3642254/93‑1 (MQ=255)

                                                                  ccAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatga                      <  2:3469546/93‑1 (MQ=255)

                                                                    agcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATc                    >  2:1280931/1‑93 (MQ=255)

                                                                    agcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATc                    >  2:3377809/1‑93 (MQ=255)

                                                                    agcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATc                    >  2:6033760/1‑93 (MQ=255)

                                                                      cagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGa                  >  2:1490247/1‑93 (MQ=255)

                                                                       agGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGaa                 >  1:2371180/1‑93 (MQ=255)

                                                                       agGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGaa                 >  2:1686890/1‑93 (MQ=255)

                                                                        ggCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGaaa                <  1:2428348/93‑1 (MQ=255)

                                                                        ggCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGaaa                >  1:6440100/1‑93 (MQ=255)

                                                                        ggCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGaaa                >  1:1801810/1‑93 (MQ=255)

                                                                        ggCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGaaa                >  2:4853731/1‑93 (MQ=255)

                                                                          ccAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAata              <  1:41753/93‑1 (MQ=255)

                                                                           cAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAatat             <  1:3693730/93‑1 (MQ=255)

                                                                           cAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAatat             <  1:1894540/93‑1 (MQ=255)

                                                                             aTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTc           >  1:1561534/1‑93 (MQ=255)

                                                                             aTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTc           <  2:1300352/93‑1 (MQ=255)

                                                                                 ccACCCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtct       <  2:5631821/93‑1 (MQ=255)

                                                                                  cAACCCGGGCGTTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtctc      <  1:6283447/93‑1 (MQ=255)

                                                                                   aaCCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTcc     <  2:4188229/93‑1 (MQ=255)

TAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGT  >  NC_000913/477007‑477185

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: `‑`