BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,786,415	A→G	T230A (ACG→GCG)	b3614 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,786,415	0	A	G	98.1%	663.9 / 0.4	210	T230A (ACG→GCG)	b3614	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/1); new base G (113/93); total (116/94)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 216 total aligned reads displayed.

GCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAT  >  NC_000913/3786326‑3786504

gCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTgcgc                                                                                        >  2:5633/1‑93 (MQ=255)

gCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTgcgc                                                                                        <  1:1984874/93‑1 (MQ=255)

 cAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCt                                                                                       <  1:1331604/93‑1 (MQ=255)

   aCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTgg                                                                                     <  1:1000389/93‑1 (MQ=255)

   aCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGGGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTgg                                                                                     >  2:474373/1‑93 (MQ=255)

    cGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGt                                                                                    <  1:969365/93‑1 (MQ=255)

    cGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGt                                                                                    <  2:3037657/93‑1 (MQ=255)

      gTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGa                                                                                  >  1:3697808/1‑93 (MQ=255)

        tCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGaaa                                                                                <  1:4129603/93‑1 (MQ=255)

        tCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGaaa                                                                                >  1:1225643/1‑93 (MQ=255)

        tCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGaaa                                                                                <  1:589749/93‑1 (MQ=255)

          gACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGc                                                                              >  2:4101970/1‑93 (MQ=255)

          gACCATTGCCATTGGGCATACTCACCCTGCAACGGTACTCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGc                                                                              >  2:1211250/1‑93 (MQ=255)

           aCCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGcc                                                                             <  2:72413/93‑1 (MQ=255)

           aCCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGcc                                                                             <  2:2943817/93‑1 (MQ=255)

            ccattgccattgGGCGTCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCca                                                                            >  2:1675246/1‑93 (MQ=255)

              attgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagc                                                                          <  2:4434749/93‑1 (MQ=255)

               ttgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagca                                                                         <  2:3882273/93‑1 (MQ=255)

                tgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagcag                                                                        >  1:2297832/1‑93 (MQ=255)

                 tccattgGTCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCagcagc                                                                       <  1:3181321/92‑1 (MQ=255)

                 gccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCagcagc                                                                       >  1:2207711/1‑93 (MQ=255)

                 gccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCagcagc                                                                       >  2:3609174/1‑93 (MQ=255)

                   cattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCtt                                                                     >  2:1750436/1‑93 (MQ=255)

                    attgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTg                                                                    >  1:3843404/1‑93 (MQ=255)

                     ttgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGc                                                                   <  1:3568855/93‑1 (MQ=255)

                     ttgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGc                                                                   <  1:1105046/93‑1 (MQ=255)

                      tgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCt                                                                  <  1:607596/93‑1 (MQ=255)

                       gggCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtg                                                                 >  1:1520179/1‑93 (MQ=255)

                       gggCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtg                                                                 >  2:3240937/1‑93 (MQ=255)

                       gggCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtg                                                                 >  1:1392218/1‑93 (MQ=255)

                        ggCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtga                                                                >  2:929017/1‑93 (MQ=255)

                           aTCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatg                                                             >  2:4551594/1‑93 (MQ=255)

                           aTCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatg                                                             >  1:4261059/1‑93 (MQ=255)

                            tCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatga                                                            >  1:866361/1‑93 (MQ=255)

                              cTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTGTAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAAc                                                          >  1:3987655/1‑93 (MQ=255)

                              cTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAAc                                                          >  1:2449022/1‑93 (MQ=255)

                              cTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAAc                                                          >  1:2395377/1‑93 (MQ=255)

                                cACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCg                                                        >  2:2657404/1‑93 (MQ=255)

                                 ccccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGc                                                       <  2:3815676/92‑1 (MQ=255)

                                 aCCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGc                                                       >  1:4445770/1‑93 (MQ=255)

                                  cccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCa                                                      <  1:4158302/93‑1 (MQ=255)

                                   ccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAg                                                     <  1:4408713/93‑1 (MQ=255)

                                   ccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAg                                                     >  2:811591/1‑93 (MQ=255)

                                   ccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAg                                                     >  1:3673277/1‑93 (MQ=255)

                                   ccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAg                                                     <  2:1427159/93‑1 (MQ=255)

                                    cTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAgg                                                    >  1:2896006/1‑93 (MQ=255)

                                     tGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGt                                                   >  1:190870/1‑93 (MQ=255)

                                      gCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGtt                                                  >  2:3543000/1‑93 (MQ=255)

                                      gCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGtt                                                  >  2:1506820/1‑93 (MQ=255)

                                        aaCGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGa                                                >  1:1841497/1‑93 (MQ=255)

                                         aCGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGAt                                               <  2:2828078/93‑1 (MQ=255)

                                         aCGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGAt                                               >  2:143089/1‑93 (MQ=255)

                                          cGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATa                                              >  2:2280251/1‑93 (MQ=255)

                                           ggTACGCGTGTTGCAACAGATGGTTGATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATAc                                             <  1:103144/93‑1 (MQ=255)

                                            gTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACt                                            >  2:1897130/1‑93 (MQ=255)

                                             tACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACtt                                           >  1:679551/1‑93 (MQ=255)

                                              aCGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTc                                          <  2:2249863/93‑1 (MQ=255)

                                              aCGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTc                                          >  2:4242817/1‑93 (MQ=255)

                                              aCGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTc                                          <  2:4383524/93‑1 (MQ=255)

                                               cgcgTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCt                                         >  2:1310379/1‑93 (MQ=255)

                                               cgcgTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCt                                         <  2:2674051/93‑1 (MQ=255)

                                                gcgTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTa                                        <  1:531480/93‑1 (MQ=255)

                                                gcgTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTa                                        >  1:295710/1‑93 (MQ=255)

                                                 cgTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTac                                       <  1:4527110/93‑1 (MQ=255)

                                                  gtgtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTaca                                      <  2:2865190/93‑1 (MQ=255)

                                                   tgtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTacac                                     <  1:4547041/93‑1 (MQ=255)

                                                    gtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAcc                                    >  2:4300415/1‑93 (MQ=255)

                                                        cAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgccg                                <  1:4694287/93‑1 (MQ=255)

                                                        cAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgccg                                <  2:3331896/93‑1 (MQ=255)

                                                        cAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgccg                                >  1:1300178/1‑93 (MQ=255)

                                                         aaCAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGa                               >  2:3186127/1‑93 (MQ=255)

                                                         aaCAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGa                               >  2:1339949/1‑93 (MQ=255)

                                                           cAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaa                             <  1:3006341/93‑1 (MQ=255)

                                                             gATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaaaa                           >  1:4705023/1‑93 (MQ=255)

                                                              aTGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAAc                          >  1:4314943/1‑93 (MQ=255)

                                                              aTGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAAc                          <  2:1952187/93‑1 (MQ=255)

                                                              aTGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAAc                          <  2:3133794/93‑1 (MQ=255)

                                                              aTGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAAc                          >  1:773654/1‑93 (MQ=255)

                                                              aTGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTGCTACACCGCCGAAAAAc                          >  1:4223033/1‑93 (MQ=255)

                                                               tGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACg                         >  1:1437249/1‑93 (MQ=255)

                                                               tGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACg                         >  2:2737660/1‑93 (MQ=255)

                                                                ggTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGc                        >  1:3962687/1‑93 (MQ=255)

                                                                 gTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCt                       >  2:651513/1‑93 (MQ=255)

                                                                 gTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCt                       <  2:1988635/93‑1 (MQ=255)

                                                                 gTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCt                       >  2:3200524/1‑93 (MQ=255)

                                                                  tttATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtg                      <  1:1544189/93‑1 (MQ=255)

                                                                    tataACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtgtg                    >  1:2377523/1‑93 (MQ=255)

                                                                     ataACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGc                   <  1:3324382/93‑1 (MQ=255)

                                                                      taACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGcc                  <  1:624991/93‑1 (MQ=255)

                                                                       aaCCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCt                 <  2:3950863/93‑1 (MQ=255)

                                                                       aaCCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCt                 >  1:4595582/1‑93 (MQ=255)

                                                                       aaCCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCt                 <  1:2704069/93‑1 (MQ=255)

                                                                         ccTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGa               >  2:2118542/1‑93 (MQ=255)

                                                                          cTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGAt              <  2:879174/93‑1 (MQ=255)

                                                                                ccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcg        >  1:2968484/1‑93 (MQ=255)

                                                                                ccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcg        >  1:3053410/1‑93 (MQ=255)

                                                                                   gACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGt     >  2:3838851/1‑93 (MQ=255)

                                                                                    aCATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTa    >  2:2485963/1‑93 (MQ=255)

                                                                                     cATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACAGCCGAAAAACGCTGTGCCTGATGCACCGCGTaa   <  2:2197124/93‑1 (MQ=255)

                                                                                      aTTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAt  <  2:1791075/93‑1 (MQ=255)

GCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAT  >  NC_000913/3786326‑3786504

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