BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,025,930	G→A	S42F (TCC→TTC)	b1956 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,025,930	0	G	A	96.0%	423.3 / 2.1	124	S42F (TCC→TTC)	b1956	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (2/3); new base A (58/61); total (60/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 129 total aligned reads displayed.

GAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCC  >  NC_000913/2025839‑2026019

gAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGa                                                                                          <  2:3108074/93‑1 (MQ=255)

 aGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGaaa                                                                                         >  2:4339448/1‑93 (MQ=255)

  gcgcGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGaaaa                                                                                        >  2:2224674/1‑93 (MQ=255)

  gcgcGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGaaaa                                                                                        >  1:2752161/1‑93 (MQ=255)

  gcgcGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGaaaa                                                                                        <  2:1377684/93‑1 (MQ=255)

   cgcgTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAaaaa                                                                                       <  1:2872081/93‑1 (MQ=255)

      gTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAg                                                                                    <  1:3247685/93‑1 (MQ=255)

      gTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAg                                                                                    >  2:3187912/1‑93 (MQ=255)

       ttGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGc                                                                                   >  2:441521/1‑93 (MQ=255)

        tGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCa                                                                                  <  2:378032/93‑1 (MQ=255)

        tGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCa                                                                                  >  1:1693464/1‑93 (MQ=255)

        tGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCa                                                                                  >  1:266637/1‑93 (MQ=255)

          gCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAAt                                                                                <  1:1491794/93‑1 (MQ=255)

           cAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATa                                                                               <  2:3680147/93‑1 (MQ=255)

           cAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATa                                                                               <  1:510987/93‑1 (MQ=255)

           cAACGTTGTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATa                                                                               <  1:130797/93‑1 (MQ=255)

            aaCGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGCGAGCAAGGTGAAAAAAAGCAATAc                                                                              >  1:3265796/1‑93 (MQ=255)

             aCGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATAcc                                                                             <  1:2700816/93‑1 (MQ=255)

               gTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCAc                                                                           >  2:4077860/1‑93 (MQ=255)

                ttttCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACg                                                                          >  1:550444/1‑93 (MQ=255)

                 tttCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGa                                                                         <  2:2462406/93‑1 (MQ=255)

                 tttCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGa                                                                         >  2:1422045/1‑93 (MQ=255)

                 tttCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGa                                                                         <  2:375198/93‑1 (MQ=255)

                     cAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATaaa                                                                     <  1:4072125/93‑1 (MQ=255)

                     cAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATaaa                                                                     <  1:2179950/93‑1 (MQ=255)

                       gatgatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGc                                                                   >  1:248377/1‑93 (MQ=255)

                       gatgatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGc                                                                   <  2:408919/93‑1 (MQ=255)

                        atgatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCa                                                                  <  2:451345/93‑1 (MQ=255)

                         tgatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCag                                                                 >  1:462031/1‑93 (MQ=255)

                         tgatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCag                                                                 >  2:501168/1‑93 (MQ=255)

                          gatTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCaga                                                                >  1:4534981/1‑93 (MQ=255)

                            ttACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAgaga                                                              <  1:2408933/93‑1 (MQ=255)

                            ttACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGAagagt                                                              >  1:825017/1‑92 (MQ=255)

                             tACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGAt                                                             <  1:2667808/93‑1 (MQ=255)

                              aCGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGAtt                                                            <  1:1625087/93‑1 (MQ=255)

                              aCGCTGGCTGGAGATATAGGCATCTACCAGCACCACCACTACACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGAtt                                                            <  2:2598484/93‑1 (MQ=255)

                                gctggctgGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTaa                                                          <  1:780984/93‑1 (MQ=255)

                                gctggctgGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTaa                                                          >  2:2292128/1‑93 (MQ=255)

                                gctggctgGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTaa                                                          <  2:330014/93‑1 (MQ=255)

                                  tggctggAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACg                                                        >  2:1281030/1‑93 (MQ=255)

                                   ggctggAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGa                                                       >  1:2331319/1‑93 (MQ=255)

                                   ggctggAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGa                                                       <  1:4255370/93‑1 (MQ=255)

                                     ctggAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACa                                                     >  2:1493970/1‑93 (MQ=255)

                                      tggAGATATAGGAATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACAt                                                    <  1:1161790/93‑1 (MQ=255)

                                        gagaTATAGGCATCTTCCAGCACCACCACTTCGCGCCAGGTGAGCAAGGTGAAAACAAGCAATACCACGATCAAGAGGAGGTTAACGACATGa                                                  >  1:3503264/1‑93 (MQ=255)

                                        gagaTATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATgc                                                  >  1:3503263/1‑92 (MQ=255)

                                         agaTATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGAc                                                 >  1:3810970/1‑93 (MQ=255)

                                          gaTATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGAcc                                                >  2:2035108/1‑93 (MQ=255)

                                          gaTATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGAcc                                                <  2:2672602/93‑1 (MQ=255)

                                          gaTATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGAcc                                                >  1:1116279/1‑93 (MQ=255)

                                            tataGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAg                                              <  1:457773/93‑1 (MQ=255)

                                            tataGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAg                                              >  1:103488/1‑93 (MQ=255)

                                            tataGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAg                                              <  2:2029995/93‑1 (MQ=255)

                                            tataGGCATCATCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAg                                              <  2:2679864/93‑1 (MQ=255)

                                             ataGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAgg                                             >  1:4624128/1‑93 (MQ=255)

                                             ataGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAgg                                             <  2:735153/93‑1 (MQ=255)

                                                 gCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGcccc                                         <  1:2006781/93‑1 (MQ=255)

                                                  cATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCa                                        >  1:4462362/1‑93 (MQ=255)

                                                  cATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCa                                        <  1:1626487/93‑1 (MQ=255)

                                                  cATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCa                                        <  2:1601466/93‑1 (MQ=255)

                                                   aTCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAg                                       <  2:219334/93‑1 (MQ=255)

                                                   aTCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAg                                       >  2:797474/1‑93 (MQ=255)

                                                    tCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAgg                                      <  2:2257475/93‑1 (MQ=255)

                                                    tCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAgg                                      >  1:1476681/1‑93 (MQ=255)

                                                     cTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAggc                                     <  1:3620971/93‑1 (MQ=255)

                                                     cTTCCAGCACCACCACTGCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAggc                                     <  2:166895/93‑1 (MQ=255)

                                                      ttCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAggcg                                    >  2:3055662/1‑93 (MQ=255)

                                                       tCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAggcgg                                   <  2:3585341/93‑1 (MQ=255)

                                                         cAGCACCACCACTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGgcggcg                                 >  2:444587/1‑93 (MQ=255)

                                                            caccaccacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGc                              <  1:3304649/93‑1 (MQ=255)

                                                               aaccacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCTGAGATTAACGACATGACAAGGCCCCAGGCGGCGTGCGAg                           >  1:313219/2‑93 (MQ=255)

                                                                accacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCATGCGGCGTGCGAGt                          >  2:3092579/1‑93 (MQ=255)

                                                                accacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGt                          >  2:288233/1‑93 (MQ=255)

                                                                  cacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGttt                        >  1:4115626/1‑93 (MQ=255)

                                                                  cacTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGttt                        <  1:776945/93‑1 (MQ=255)

                                                                   acTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGtttt                       >  2:3637219/1‑93 (MQ=255)

                                                                   acTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGtttt                       >  1:3857200/1‑93 (MQ=255)

                                                                   acTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGtttt                       <  1:346209/93‑1 (MQ=255)

                                                                    cTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGttttt                      >  1:963119/1‑93 (MQ=255)

                                                                    cTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGttttt                      <  2:4413586/93‑1 (MQ=255)

                                                                    cTTCACGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGttttt                      >  2:123125/1‑93 (MQ=255)

                                                                         cGCCAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAAc                 <  2:4681120/93‑1 (MQ=255)

                                                                           ccAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCa               >  1:4306029/1‑93 (MQ=255)

                                                                           ccAGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCa               <  2:1026780/93‑1 (MQ=255)

                                                                             aGGTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGc             <  2:391759/93‑1 (MQ=255)

                                                                              ggTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGct            >  1:3872789/1‑93 (MQ=255)

                                                                              ggTGAGCAAGGTGAAAAAAAGCAAAACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGct            <  2:3376635/93‑1 (MQ=255)

                                                                               gTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGctc           <  2:3934862/93‑1 (MQ=255)

                                                                               gTGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGctc           >  1:489389/1‑93 (MQ=255)

                                                                                tGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGctct          <  1:2877142/93‑1 (MQ=255)

                                                                                tGAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGctct          >  2:779013/1‑93 (MQ=255)

                                                                                 gAGCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTg         <  2:4443086/93‑1 (MQ=255)

                                                                                   gCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGt       >  1:4077359/1‑93 (MQ=255)

                                                                                   gCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGt       <  1:926096/93‑1 (MQ=255)

                                                                                   gCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGt       <  2:1009718/93‑1 (MQ=255)

                                                                                   gCAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGt       >  2:2380721/1‑93 (MQ=255)

                                                                                    cAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGtt      <  1:1441806/93‑1 (MQ=255)

                                                                                    cAAGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGtt      >  2:3280056/1‑93 (MQ=255)

                                                                                     aaGGTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGttt     >  2:1884309/1‑93 (MQ=255)

                                                                                        gTGAAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTcc  >  2:2057618/1‑93 (MQ=255)

GAGCGCGTTGGCAACGTTTTCCAGATGATTACGCTGGCTGGAGATATAGGCATCTTCCAGCACCACCACTTCACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCC  >  NC_000913/2025839‑2026019

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