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	99.8%	1666.2 / ‑5.3	511	T230A (ACG→GCG)	b3614	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (283/227); total (283/228)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 523 total aligned reads displayed.

CGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATC  >  NC_000913/3786325‑3786505

  cAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCt                                                                                        >  1:8389788/1‑93 (MQ=255)

   aaCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTg                                                                                       >  1:3915668/1‑93 (MQ=255)

    aCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTgg                                                                                      >  1:4734177/1‑93 (MQ=255)

    aCGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTgg                                                                                      >  2:9097211/1‑93 (MQ=255)

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

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

      ggTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTg                                                                                    >  2:58715/1‑93 (MQ=255)

        gtCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGGACGCGTGTCGCTACAGATGGATTATAAGCTGCCGCGAGACAAAAAGCTGGCGaa                                                                                  >  1:1266629/2‑93 (MQ=255)

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

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

            aCCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGcc                                                                              >  2:3604484/1‑93 (MQ=255)

            aCCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGcc                                                                              >  1:5484801/1‑93 (MQ=255)

             ccattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCca                                                                             <  2:1381131/93‑1 (MQ=255)

             ccattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCca                                                                             >  2:6485395/1‑93 (MQ=255)

             ccattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCca                                                                             >  2:8655222/1‑93 (MQ=255)

             ccattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCca                                                                             >  2:9306469/1‑93 (MQ=255)

              cattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcag                                                                            <  2:2990758/93‑1 (MQ=255)

              cattgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcag                                                                            <  2:4805178/93‑1 (MQ=255)

                ttgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagca                                                                          >  1:5433284/1‑93 (MQ=255)

                ttgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagca                                                                          >  1:2256897/1‑93 (MQ=255)

                ttgccattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCcagca                                                                          >  2:9007087/1‑93 (MQ=255)

                    cattgGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCtt                                                                      <  1:9302770/93‑1 (MQ=255)

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

                         ggCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtga                                                                 <  1:9629922/93‑1 (MQ=255)

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

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

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

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

                             tCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCtgaatga                                                             <  1:177972/93‑1 (MQ=255)

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

                               cTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAAc                                                           <  1:8103865/93‑1 (MQ=255)

                               cTCACCCGGCAACGGTACGCGTGTGGCAACAGATGGTGGAGAACCTGCCGCCAGACAACGCGCGGGCGAAAGCCAGCAGCTGGCTGAATGAAc                                                           >  2:7722044/1‑93 (MQ=255)

                                tCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAAcc                                                          <  1:522467/93‑1 (MQ=255)

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

                                   cccTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCa                                                       >  1:6892023/1‑93 (MQ=255)

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

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

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

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

                                     cTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAgg                                                     <  2:3405758/93‑1 (MQ=255)

                                      tGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGt                                                    <  1:5746754/93‑1 (MQ=255)

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

                                      tGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGt                                                    <  2:2200209/93‑1 (MQ=255)

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

                                        cAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTg                                                  >  2:808283/1‑93 (MQ=255)

                                         aaCGGTTCTGGTGTTGCAACAGATGGTTGATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGa                                                 <  2:4590810/93‑1 (MQ=255)

                                         aaCGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGa                                                 <  2:7910388/93‑1 (MQ=255)

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

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

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

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

                                            ggTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATAc                                              <  1:1158670/93‑1 (MQ=255)

                                            ggTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATAc                                              <  1:3123016/93‑1 (MQ=255)

                                             gTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACt                                             <  2:9913109/93‑1 (MQ=255)

                                              tACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACtt                                            <  1:2457412/93‑1 (MQ=255)

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

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

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

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

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

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

                                                   gtgtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTaca                                       >  1:5623845/1‑93 (MQ=255)

                                                    tgtTGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTacac                                      >  1:9794889/1‑93 (MQ=255)

                                                      ttGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccg                                    <  1:8917509/93‑1 (MQ=255)

                                                       tGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgc                                   >  1:1449150/1‑93 (MQ=255)

                                                       tGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgc                                   <  2:7273041/93‑1 (MQ=255)

                                                       tGCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgc                                   <  2:6014827/93‑1 (MQ=255)

                                                        gCAACAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAccgcc                                  <  1:7424668/93‑1 (MQ=255)

                                                        gCAACAGATGGTTAATAACCTGCCGCCAGACATAGCGCTGGTGAAAGCCAGCAGATTGCTGAATAAACCGCAGGTTGAAACTTCAAAAccgcc                                  <  2:345988/93‑1 (MQ=255)

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

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

                                                           aCAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaa                               >  1:4992426/1‑93 (MQ=255)

                                                            cAGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaa                              >  1:4113717/1‑93 (MQ=255)

                                                             aGATGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGaaaa                             >  1:6771823/1‑93 (MQ=255)

                                                                tGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACg                          <  1:7717144/93‑1 (MQ=255)

                                                                tGGTTTATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAACGCCGAAAAACg                          >  1:3671900/1‑93 (MQ=255)

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

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

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

                                                                   tttATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtg                       >  2:7040036/1‑93 (MQ=255)

                                                                    ttATAACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtgt                      <  2:2510464/93‑1 (MQ=255)

                                                                    ttATAACCGGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCtgt                      <  1:4212615/93‑1 (MQ=255)

                                                                      ataACCTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGc                    >  1:7632376/1‑93 (MQ=255)

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

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

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

                                                                           cTGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGAt               >  1:3794966/1‑93 (MQ=255)

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

                                                                            tGCCGCCAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATg              <  2:6896855/93‑1 (MQ=255)

                                                                             gccgccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGc             >  2:8350494/1‑93 (MQ=255)

                                                                             gccgccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGc             <  1:6043462/93‑1 (MQ=255)

                                                                                gccAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCAcc          >  1:5049895/1‑93 (MQ=255)

                                                                                  cAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcgc        >  2:4206847/1‑93 (MQ=255)

                                                                                  cAGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcgc        <  2:9390469/93‑1 (MQ=255)

                                                                                   aGACATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcgcg       >  1:7955107/1‑93 (MQ=255)

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

                                                                                      cATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTaa    <  2:3668144/93‑1 (MQ=255)

                                                                                      cATTGCGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTaa    <  1:1577154/93‑1 (MQ=255)

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

                                                                                        ttACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtc  <  1:9935756/93‑1 (MQ=255)

CGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATC  >  NC_000913/3786325‑3786505

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