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.3%	462.8 / ‑4.6	150	L251P (CTC→CCC)	ylaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (73/76); total (73/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

Only 100 of 158 total aligned reads displayed.

ATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGT  >  NC_000913/477006‑477185

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

   gCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGt                                                                                      <  2:860092/93‑1 (MQ=255)

     aTCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTaa                                                                                    <  1:2546948/93‑1 (MQ=255)

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

      tCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAc                                                                                   <  1:3065150/93‑1 (MQ=255)

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

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

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

      tCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACACGGGCGGTAAc                                                                                   >  1:3114443/1‑93 (MQ=255)

        cTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCa                                                                                 >  2:3969848/1‑93 (MQ=255)

        cTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCa                                                                                 >  1:2798653/1‑93 (MQ=255)

        cTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCa                                                                                 >  2:700351/1‑93 (MQ=255)

          gCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCaaa                                                                               >  2:352425/1‑93 (MQ=255)

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

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

           cAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAAt                                                                              >  1:2901330/1‑93 (MQ=255)

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

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

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

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

              ccGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATaaa                                                                           >  1:2710968/1‑93 (MQ=255)

              ccGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATaaa                                                                           <  1:3790522/93‑1 (MQ=255)

                gatgatGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGa                                                                         >  2:136424/1‑93 (MQ=255)

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

                 atgatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCCACCCGGGCGGTAACCAAATAAAGAc                                                                        >  2:1534754/1‑93 (MQ=255)

                 atgatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGAc                                                                        <  2:743937/93‑1 (MQ=255)

                  tgatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACt                                                                       >  1:2914267/1‑93 (MQ=255)

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

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

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

                   gatgCGGGGACTGAATACGACGCAAAATACGCAGCACAAACATAGCTGCCAGCAGGACAATCACCAACCCAGGCGGTAAACAAATAAAGAGtt                                                                      >  2:2032748/1‑93 (MQ=255)

                     tgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGc                                                                    >  2:561141/1‑93 (MQ=255)

                      gCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGcc                                                                   <  1:2653878/93‑1 (MQ=255)

                      gCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGcc                                                                   <  1:2536599/93‑1 (MQ=255)

                      gCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGcc                                                                   <  1:1787041/93‑1 (MQ=255)

                        ggTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGa                                                                 <  2:925607/93‑1 (MQ=255)

                         gTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGAt                                                                <  1:1820241/93‑1 (MQ=255)

                          tGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATg                                                               >  2:3991141/1‑93 (MQ=255)

                            aCTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGcc                                                             >  1:200637/1‑93 (MQ=255)

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

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

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

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

                                 aTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCt                                                        <  1:682842/93‑1 (MQ=255)

                                 aTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCt                                                        <  1:2099635/93‑1 (MQ=255)

                                 aTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCt                                                        <  1:1392527/93‑1 (MQ=255)

                                 aTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCt                                                        <  1:4196661/93‑1 (MQ=255)

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

                                      gcgcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTg                                                   >  1:1447410/1‑93 (MQ=255)

                                      gcgcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTg                                                   >  1:3672365/1‑93 (MQ=255)

                                      gcgcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTg                                                   >  2:1473718/1‑93 (MQ=255)

                                        gcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCAGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCa                                                 >  1:4522455/1‑93 (MQ=255)

                                         cAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCaa                                                >  1:104975/1‑93 (MQ=255)

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

                                           aaaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACa                                              <  2:2476289/93‑1 (MQ=255)

                                           aaaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACa                                              <  1:775610/93‑1 (MQ=255)

                                           aaaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAAACCGGGCGGTAACCAAATAAAGACTTGCCGATGACAGCCTTTCTGCAACa                                              >  2:4508329/1‑93 (MQ=255)

                                              tACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAttt                                           >  2:2719815/1‑93 (MQ=255)

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

                                               aCGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAtttt                                          <  1:2662747/93‑1 (MQ=255)

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

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

                                                 gcagcaCAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCg                                        <  1:4005995/93‑1 (MQ=255)

                                                 gcagcaCAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCg                                        <  1:87545/93‑1 (MQ=255)

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

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

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

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

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

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

                                                        aaaCATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGc                                 >  2:3194766/1‑93 (MQ=255)

                                                        aaaCATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGc                                 >  1:1646600/1‑93 (MQ=255)

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

                                                            aTCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCag                             >  1:1381952/1‑92 (MQ=255)

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

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

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

                                                               gCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtgtg                          <  1:831188/93‑1 (MQ=255)

                                                                cTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtga                         >  1:1677276/1‑93 (MQ=255)

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

                                                                cTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtga                         >  1:1502525/1‑93 (MQ=255)

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

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

                                                                   ccAGCAGTCCAATCACCACCCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatga                      <  2:3310797/93‑1 (MQ=255)

                                                                       cagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCCACCTTTTCGTTCAAGCCCATGTGATGATCGa                  >  2:3342844/1‑93 (MQ=255)

                                                                       cagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGACCATGTGATGATCGa                  >  1:995251/1‑93 (MQ=255)

                                                                        agGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGaa                 <  1:631697/93‑1 (MQ=255)

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

                                                                            cAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAatat             >  1:2278201/1‑93 (MQ=255)

                                                                             aaTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATAtt            >  2:3188643/1‑93 (MQ=255)

                                                                             aaTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATAtt            <  1:4464982/93‑1 (MQ=255)

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

                                                                                 aCCATCCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCATCATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtc        <  1:4472927/93‑1 (MQ=255)

                                                                                 aCCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtc        <  1:2622215/93‑1 (MQ=255)

                                                                                 aCCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtc        >  1:1130849/1‑93 (MQ=255)

                                                                                 aCCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtc        >  2:864124/1‑93 (MQ=255)

                                                                                  ccAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtct       >  2:228222/1‑93 (MQ=255)

                                                                                    aaCCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTcc     >  2:3637588/1‑93 (MQ=255)

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

                                                                                     aCCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCg    >  1:1793087/1‑93 (MQ=255)

ATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGT  >  NC_000913/477006‑477185

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