BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,785,180	Δ1 bp	coding (413/495 nt)	b2655 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,785,180	0	T	.	99.4%	819.6 / ‑1.8	181	coding (413/495 nt)	b2655	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); major base . (103/76); minor base A (1/0); total (104/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 = 9.99e-01

Only 100 of 185 total aligned reads displayed.

ACCTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCGTCCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAATCCAA  >  NC_000913/2785092‑2785267

aCCTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGAc                                                                                    >  1:85500/1‑93 (MQ=255)

  cTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACgg                                                                                  >  1:2668245/1‑93 (MQ=255)

   tACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGa                                                                                 >  2:1875277/1‑93 (MQ=255)

    aCCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGAc                                                                                >  2:3506893/1‑93 (MQ=255)

     ccAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGAcc                                                                               >  2:2642944/1‑93 (MQ=255)

       agcagATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAg                                                                             >  1:4637190/1‑93 (MQ=255)

       agcagATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAg                                                                             >  2:2823964/1‑93 (MQ=255)

        gcagATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGCCCAGa                                                                            >  2:1262150/1‑93 (MQ=255)

          agATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCGGCG‑CCGACGGACCAGaaa                                                                          >  2:727328/1‑93 (MQ=255)

          agATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGaaa                                                                          >  2:1336476/1‑93 (MQ=255)

          agATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGaaa                                                                          >  1:2240285/1‑93 (MQ=255)

          agATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGaaa                                                                          <  2:336029/93‑1 (MQ=255)

           gATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAAc                                                                         <  2:2755581/93‑1 (MQ=255)

           gATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAAc                                                                         >  2:614084/1‑93 (MQ=255)

              cAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTcc                                                                      <  1:3849451/93‑1 (MQ=255)

              cAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTcc                                                                      >  1:289279/1‑93 (MQ=255)

               agagAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTccc                                                                     >  2:3121045/1‑93 (MQ=255)

                gagaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCg                                                                    >  2:1643456/1‑93 (MQ=255)

                gagaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCg                                                                    <  2:1739397/93‑1 (MQ=255)

                gagaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCg                                                                    >  1:3507183/1‑93 (MQ=255)

                 agaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCga                                                                   <  2:1187585/93‑1 (MQ=255)

                 agaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCga                                                                   <  1:4078622/93‑1 (MQ=255)

                  gaAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaa                                                                  >  1:2227792/1‑93 (MQ=255)

                   aaaaTATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaag                                                                 >  1:1121583/1‑93 (MQ=255)

                   aaaaTATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaag                                                                 >  1:1221504/1‑93 (MQ=255)

                    aaaTATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaaga                                                                <  2:4460035/93‑1 (MQ=255)

                     aaTATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaa                                                               >  2:2699584/1‑93 (MQ=255)

                     aaTATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaa                                                               <  2:2300954/93‑1 (MQ=255)

                      atatGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaag                                                              >  2:927640/1‑93 (MQ=255)

                      atatGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaag                                                              >  2:4201307/1‑93 (MQ=255)

                      atatGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaag                                                              <  2:3448028/93‑1 (MQ=255)

                       tatGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaaga                                                             >  2:1426881/1‑93 (MQ=255)

                        atGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaagaa                                                            <  2:2483119/93‑1 (MQ=255)

                        atGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCgaagaagaa                                                            >  2:1134054/1‑93 (MQ=255)

                         tGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGaaa                                                           <  2:1487938/93‑1 (MQ=255)

                         tGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGaaa                                                           <  2:4234255/93‑1 (MQ=255)

                           tCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAg                                                         >  2:2837165/1‑93 (MQ=255)

                            cTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAgg                                                        <  1:1310612/93‑1 (MQ=255)

                            cTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAgg                                                        <  2:4368751/93‑1 (MQ=255)

                            cTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAgg                                                        <  1:97597/93‑1 (MQ=255)

                             tGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCGGAAACTCCCGAAGAAGAAAAggg                                                       >  2:2644755/1‑93 (MQ=255)

                             tGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAggg                                                       >  2:3096501/1‑93 (MQ=255)

                             tGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAggg                                                       >  2:2243809/1‑93 (MQ=255)

                              gAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCTGAAGAAGAAAAGGGt                                                      <  1:240503/93‑1 (MQ=255)

                                aTTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTa                                                    >  2:1368613/1‑93 (MQ=255)

                                 tttGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTAc                                                   >  2:4737241/1‑93 (MQ=255)

                                 tttGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGGGACAATTTTCCGC‑TCCGACGGTTCAGAAACTCCCGAATAAGAAAAgggcaat                                                   >  1:1698325/1‑89 (MQ=255)

                                  ttGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACt                                                  <  1:268313/93‑1 (MQ=255)

                                    gCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTgg                                                <  1:2724550/93‑1 (MQ=255)

                                      aCGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGtt                                              >  1:2280891/1‑93 (MQ=255)

                                       cGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGttt                                             >  2:452356/1‑93 (MQ=255)

                                         ttGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCa                                           <  2:296743/93‑1 (MQ=255)

                                          tGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCAt                                          <  2:3572057/93‑1 (MQ=255)

                                            ccATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGc                                        <  1:1990791/93‑1 (MQ=255)

                                            ccATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGc                                        <  2:3827163/93‑1 (MQ=255)

                                             cATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCt                                       >  1:1963105/1‑93 (MQ=255)

                                              aTTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGATACTGGTTTCATGCtt                                      >  2:1992403/1‑93 (MQ=255)

                                              aTTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑ACGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCtt                                      >  2:2518657/1‑93 (MQ=255)

                                               ttATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTa                                     <  1:3501747/93‑1 (MQ=255)

                                                 aTGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATc                                   <  2:3592085/93‑1 (MQ=255)

                                                      tatataCAAATCCCTCCAGTGACAACTTTCCGCGT‑CGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAccc                              <  2:3939743/93‑1 (MQ=255)

                                                      tatataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAcccc                              <  2:1156173/93‑1 (MQ=255)

                                                      tatataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAccc                              >  1:558612/1‑93 (MQ=255)

                                                      tatataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAAccc                              >  1:3222121/1‑93 (MQ=255)

                                                       atataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCg                             <  1:3080217/93‑1 (MQ=255)

                                                        tataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGa                            >  2:237868/1‑93 (MQ=255)

                                                        tataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCACAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGa                            >  1:1322337/1‑93 (MQ=255)

                                                         ataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGaa                           <  2:433662/93‑1 (MQ=255)

                                                         ataCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGaa                           >  1:2479555/1‑93 (MQ=255)

                                                          taCAAATCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAg                          >  2:801931/1‑93 (MQ=255)

                                                             aaaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATc                       >  1:1140233/1‑93 (MQ=255)

                                                             aaaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATc                       <  1:1491426/93‑1 (MQ=255)

                                                             aaaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATc                       <  1:944886/93‑1 (MQ=255)

                                                             aaaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATc                       <  1:3107883/93‑1 (MQ=255)

                                                             aaaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATc                       <  1:678140/93‑1 (MQ=255)

                                                             aaaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTCTCATGCTTATCAACCCGAAGATc                       >  2:3237018/1‑93 (MQ=255)

                                                              aaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCt                      <  2:774515/93‑1 (MQ=255)

                                                              aaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCt                      <  1:2735155/93‑1 (MQ=255)

                                                              aaTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCt                      >  1:394181/1‑93 (MQ=255)

                                                               aTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCtt                     <  2:3034904/93‑1 (MQ=255)

                                                               aTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCtt                     >  2:1297786/1‑93 (MQ=255)

                                                               aTCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCtt                     <  1:3122323/93‑1 (MQ=255)

                                                                tCCCTCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCttt                    >  1:1582262/1‑93 (MQ=255)

                                                                ccccTCCAGTGACAACTTTCCGCG‑GCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCttt                    <  1:1676651/92‑1 (MQ=255)

                                                                accccccAGGGACAACTTTCCGCG‑GCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTCACTGGTTTCATGCTTATCAACCCGAAGATCttt                    <  1:1462916/92‑1 (MQ=255)

                                                                    tCCAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAc                >  2:4690492/1‑93 (MQ=255)

                                                                      cAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAcca              <  2:3992045/93‑1 (MQ=255)

                                                                      cAGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAcca              <  1:2123432/93‑1 (MQ=255)

                                                                       aGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAccat             >  1:2197197/1‑93 (MQ=255)

                                                                       aGTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAccat             >  2:4632528/1‑93 (MQ=255)

                                                                        gTGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAccatc            >  1:4436135/1‑93 (MQ=255)

                                                                         tGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAccatcc           >  2:1502963/1‑93 (MQ=255)

                                                                         tGACAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAccatcc           >  1:3383354/1‑93 (MQ=255)

                                                                           aCAACTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTAccatccat         >  1:3009002/1‑93 (MQ=255)

                                                                             aaCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATaa       <  1:1025595/93‑1 (MQ=255)

                                                                              aCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAAt      >  2:3764381/1‑93 (MQ=255)

                                                                              aCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAAt      >  1:3698543/1‑93 (MQ=255)

                                                                              aCTTTCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAAt      >  1:3333915/1‑93 (MQ=255)

                                                                                tttCCGCG‑CAGACGGACCAGAACCTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAATcc    <  2:1057320/93‑1 (MQ=255)

                                                                                  tCCGCG‑CCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAATCCaa  <  1:3255367/93‑1 (MQ=255)

ACCTACCAGCAGATCAGAGAAAATATGTCTGAATTTGCACGTTGCCATTATGGGTATATACAAATCCCTCCAGTGACAACTTTCCGCGTCCGACGGACCAGAAACTCCCGAAGAAGAAAAGGGTTACTGGTTTCATGCTTATCAACCCGAAGATCTTTGTACCATCCATAATCCAA  >  NC_000913/2785092‑2785267

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