BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,724,581	T→C	Q137R (CAG→CGG)	yiaA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,724,581	0	T	C	99.7%	928.8 / NA	289	Q137R (CAG→CGG)	yiaA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (150/138); minor base G (0/1); total (150/139)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

AGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCA  >  NC_000913/3724491‑3724670

  aaaCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCgtgt                                                                                       >  2:5898535/1‑93 (MQ=255)

  aaaCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCgtgt                                                                                       <  2:4835383/93‑1 (MQ=255)

   aaCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGtt                                                                                      >  1:5828502/1‑93 (MQ=255)

   aaCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGtt                                                                                      >  1:5751669/1‑93 (MQ=255)

   aaCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGtt                                                                                      <  1:2576504/93‑1 (MQ=255)

    aCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGttt                                                                                     <  2:4174739/93‑1 (MQ=255)

    aCAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGttt                                                                                     <  1:264406/93‑1 (MQ=255)

      aGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCt                                                                                   <  2:6290805/93‑1 (MQ=255)

      aGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCt                                                                                   <  2:5813645/93‑1 (MQ=255)

       gATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCtt                                                                                  <  1:3251288/93‑1 (MQ=255)

         ttAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCtttt                                                                                <  1:4691293/93‑1 (MQ=255)

         ttAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCtttt                                                                                <  1:3533849/93‑1 (MQ=255)

            gTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggt                                                                             >  2:46854/1‑93 (MQ=255)

            gTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggt                                                                             >  1:1822906/1‑93 (MQ=255)

             tGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggtg                                                                            <  1:1966179/93‑1 (MQ=255)

             tGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggtg                                                                            <  2:5737591/93‑1 (MQ=255)

              gAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggtgg                                                                           >  1:5634520/1‑93 (MQ=255)

              gAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTtggtgg                                                                           >  1:3298676/1‑93 (MQ=255)

               aGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTggg                                                                          >  1:2938444/1‑93 (MQ=255)

                gTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGt                                                                         >  1:587188/1‑93 (MQ=255)

                 ttAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGtt                                                                        <  2:4511859/93‑1 (MQ=255)

                  tAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGttt                                                                       <  2:873912/93‑1 (MQ=255)

                   aaGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGtttt                                                                      >  1:838351/1‑92 (MQ=255)

                   aaGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTa                                                                      <  2:3857708/93‑1 (MQ=255)

                    aGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTAt                                                                     <  2:3151010/93‑1 (MQ=255)

                     gTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTAtt                                                                    <  2:2191646/93‑1 (MQ=255)

                     gTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTAtt                                                                    <  1:4765411/93‑1 (MQ=255)

                      tAAAAGCCCGGTGACATTGGACTGAGCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCGGTGTTTCTTTTGGTGGGTTTATTc                                                                   <  2:4419278/93‑1 (MQ=255)

                       aaaaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTcc                                                                  <  1:539580/93‑1 (MQ=255)

                       aaaaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTcc                                                                  <  1:5065644/93‑1 (MQ=255)

                        aaaGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCg                                                                 <  1:1136514/93‑1 (MQ=255)

                           gCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCa                                                              >  2:5445071/1‑93 (MQ=255)

                           gCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGATTATTCAGGCa                                                              >  2:3975473/1‑93 (MQ=255)

                            cccGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCAt                                                             >  2:1745120/1‑93 (MQ=255)

                             ccGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATc                                                            >  2:5053619/1‑93 (MQ=255)

                              cGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCa                                                           >  2:274719/1‑93 (MQ=255)

                                gTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACg                                                         >  1:3940931/1‑93 (MQ=255)

                                gTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACg                                                         >  1:1459334/1‑93 (MQ=255)

                                 tCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGCATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGa                                                        >  1:4973057/1‑93 (MQ=255)

                                    caTTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAata                                                     <  2:5836893/93‑1 (MQ=255)

                                      ttGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATAtt                                                   >  2:2321394/1‑93 (MQ=255)

                                         gACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATAttttt                                                <  2:296038/93‑1 (MQ=255)

                                         gACTCACCGTGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATAttttt                                                <  2:3722397/93‑1 (MQ=255)

                                          aCTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTc                                               <  2:4384080/93‑1 (MQ=255)

                                            tGACCGTGCTTACGTGAGTTATTCGCTGTATTCTTCCTCGCTAACCCGTGTTTTTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTg                                             <  2:4054866/93‑1 (MQ=255)

                                            tGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTg                                             >  1:683342/1‑93 (MQ=255)

                                            tGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTg                                             <  2:5601048/93‑1 (MQ=255)

                                             gACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGc                                            >  1:3686699/1‑93 (MQ=255)

                                             gACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGc                                            <  2:1950264/93‑1 (MQ=255)

                                              aCCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCa                                           >  1:6024589/1‑93 (MQ=255)

                                              aCCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCa                                           <  2:616363/93‑1 (MQ=255)

                                               ccGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCAc                                          <  1:4136738/93‑1 (MQ=255)

                                                cGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCAcc                                         <  2:462967/93‑1 (MQ=255)

                                                cGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCAcc                                         >  2:1127122/1‑93 (MQ=255)

                                                    cTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTa                                     >  2:5186706/1‑93 (MQ=255)

                                                    cTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTa                                     <  1:4389020/93‑1 (MQ=255)

                                                    cTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTa                                     <  1:1037626/93‑1 (MQ=255)

                                                    cTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTa                                     <  2:3652113/93‑1 (MQ=255)

                                                      tACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACt                                   >  2:681888/1‑93 (MQ=255)

                                                       aCGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACtt                                  >  2:1209025/1‑92 (MQ=255)

                                                       aCGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTg                                  >  1:4509617/1‑93 (MQ=255)

                                                       aCGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTg                                  <  2:81409/93‑1 (MQ=255)

                                                           gAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCAcc                              >  2:4903118/1‑93 (MQ=255)

                                                           gAGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCAcc                              >  2:5639273/1‑93 (MQ=255)

                                                            aGTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCa                             >  2:2524891/1‑93 (MQ=255)

                                                             gTTATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCaa                            >  1:4881477/1‑93 (MQ=255)

                                                              ttATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCaaa                           <  1:52241/93‑1 (MQ=255)

                                                              ttATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCaaa                           >  1:5524727/1‑93 (MQ=255)

                                                              ttATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCaaa                           >  1:5287470/1‑93 (MQ=255)

                                                               tATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAAc                          <  2:3068495/93‑1 (MQ=255)

                                                               tATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAAc                          >  2:2876409/1‑93 (MQ=255)

                                                               tATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAAc                          <  1:3038036/93‑1 (MQ=255)

                                                               tATTCGCTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAAc                          <  1:2044573/93‑1 (MQ=255)

                                                                     cTGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCtt                    <  2:4992390/93‑1 (MQ=255)

                                                                      tGTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTAGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTa                   <  1:3855309/93‑1 (MQ=255)

                                                                       gTATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaa                  <  1:2341046/93‑1 (MQ=255)

                                                                        tATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaag                 >  2:4309519/1‑93 (MQ=255)

                                                                        tATTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaag                 >  2:3358595/1‑93 (MQ=255)

                                                                         aTTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaaga                >  2:3478027/1‑93 (MQ=255)

                                                                         aTTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaaga                <  1:4083918/93‑1 (MQ=255)

                                                                         aTTCTTCCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaaga                <  2:1449248/93‑1 (MQ=255)

                                                                          ttcttcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTaagaa               <  2:4559426/93‑1 (MQ=255)

                                                                            cttcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAagaaga             <  2:2423813/93‑1 (MQ=255)

                                                                            cttcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAagaaga             >  1:743594/1‑93 (MQ=255)

                                                                            cttcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAagaaga             <  1:164596/93‑1 (MQ=255)

                                                                             ttcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAgaagaa            <  1:1902813/93‑1 (MQ=255)

                                                                              tcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGaaa           <  1:962191/93‑1 (MQ=255)

                                                                              tcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGaaa           <  1:2786930/93‑1 (MQ=255)

                                                                              tcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGaaa           <  2:6227087/93‑1 (MQ=255)

                                                                              tcCTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGaaa           >  2:2591370/1‑93 (MQ=255)

                                                                               ccTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAg          >  2:2018618/1‑93 (MQ=255)

                                                                                cTGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGc         <  1:1290938/93‑1 (MQ=255)

                                                                                 tGGGTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGcc        <  1:5174276/93‑1 (MQ=255)

                                                                                  gggTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCa       >  2:607278/1‑93 (MQ=255)

                                                                                  gggTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCa       <  1:6199407/93‑1 (MQ=255)

                                                                                  gggTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCa       >  1:2183545/1‑93 (MQ=255)

                                                                                    gTAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGt     <  2:3693835/93‑1 (MQ=255)

                                                                                     tAACCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTc    <  1:4300240/93‑1 (MQ=255)

                                                                                      aaCCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTcc   >  2:5255231/1‑93 (MQ=255)

                                                                                       aCCCGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCa  >  1:4473251/1‑93 (MQ=255)

AGAAACAGATTAGTGAGTTAAGTAAAAGCCCGGTCACATTGGACTGACCGGGCTTACGTGAGTTATTCGCTGTATTCTTCCTGGGTAACCTGTGTTTCTTTTGGTGGGTTTATTCCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCA  >  NC_000913/3724491‑3724670

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