BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,924,185	A→G	G14G (GGA→GGG)	b1844 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,924,185	0	A	G	99.2%	402.6 / NA	129	G14G (GGA→GGG)	b1844	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (77/51); minor base T (0/1); total (77/52)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

AAGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCTATTAG  >  NC_000913/1924095‑1924274

aaGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCagggggg                                                                                         >  1:1280834/1‑93 (MQ=255)

aaGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCagggggg                                                                                         <  1:1356209/93‑1 (MQ=255)

 aGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCaggggggg                                                                                        >  2:3518950/1‑93 (MQ=255)

  ggCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCTATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAggggggga                                                                                       <  2:4439999/93‑1 (MQ=255)

     gATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCg                                                                                    >  1:3750551/1‑93 (MQ=255)

     gATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCg                                                                                    >  2:4455215/1‑93 (MQ=255)

     gATATCATACCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCg                                                                                    >  2:4545076/1‑93 (MQ=255)

      atatCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGt                                                                                   >  1:3238052/1‑93 (MQ=255)

      atatCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGt                                                                                   >  1:3122100/1‑93 (MQ=255)

       tatCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGtt                                                                                  <  2:1531221/93‑1 (MQ=255)

       tatCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGtt                                                                                  >  1:829481/1‑93 (MQ=255)

       tatCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGtt                                                                                  <  1:3565460/93‑1 (MQ=255)

          cATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTgag                                                                               >  2:3280341/1‑93 (MQ=255)

            ttCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGAt                                                                             <  1:3613464/93‑1 (MQ=255)

             tCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGAtt                                                                            <  1:3170157/93‑1 (MQ=255)

             tCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGAtt                                                                            <  2:2395565/93‑1 (MQ=255)

             tCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGAtt                                                                            <  1:1120702/93‑1 (MQ=255)

              ccATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGGTGAGATTg                                                                           >  2:3261157/1‑93 (MQ=255)

               cATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGc                                                                          >  2:3254477/1‑93 (MQ=255)

                aTTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGcc                                                                         >  2:1379181/1‑93 (MQ=255)

                aTTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGcc                                                                         >  1:592270/1‑93 (MQ=255)

                 ttACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCct                                                                        >  1:848621/1‑93 (MQ=255)

                  tACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCctc                                                                       <  2:200205/93‑1 (MQ=255)

                  tACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCctc                                                                       <  2:2596863/93‑1 (MQ=255)

                   aCGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCctct                                                                      <  1:690547/93‑1 (MQ=255)

                   aCGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCctct                                                                      <  1:2126606/93‑1 (MQ=255)

                   aCGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGCTTGCAGCGGGGGATCGTTGAGATTGCctct                                                                      >  1:1476495/1‑93 (MQ=255)

                     gCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGt                                                                    <  1:248330/93‑1 (MQ=255)

                     gCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGt                                                                    >  1:4610009/1‑93 (MQ=255)

                      cTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGtt                                                                   <  2:1303598/93‑1 (MQ=255)

                       tAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTg                                                                  >  2:2739816/1‑93 (MQ=255)

                        aGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGa                                                                 >  1:3247049/1‑93 (MQ=255)

                         ggCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGAt                                                                >  1:1728428/1‑93 (MQ=255)

                         ggCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGAt                                                                <  2:3371252/93‑1 (MQ=255)

                          gCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATg                                                               <  1:143916/93‑1 (MQ=255)

                          gCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATg                                                               >  2:805875/1‑93 (MQ=255)

                          gCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATg                                                               >  2:1792499/1‑93 (MQ=255)

                          gCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATg                                                               <  1:4320463/93‑1 (MQ=255)

                          gCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATg                                                               >  2:2011278/1‑93 (MQ=255)

                            tttttCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTc                                                             >  1:739730/1‑93 (MQ=255)

                            tttttCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTc                                                             <  1:3992835/93‑1 (MQ=255)

                             ttttCGGCCTGGAGCATGCTATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCa                                                            <  2:1027873/93‑1 (MQ=255)

                              tttCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCAt                                                           <  1:1490815/93‑1 (MQ=255)

                               ttCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCAtt                                                          <  1:1409020/93‑1 (MQ=255)

                               ttCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCAtt                                                          <  2:2347688/93‑1 (MQ=255)

                               ttCGCCATGGAGCATGCCATGTTGCGCATTATCGATACAGAAACTTGCGGTTTGCAGGGTGGGATCGTTGAGATTGCCTCTGTTGATGTCAtt                                                          <  2:2234362/93‑1 (MQ=255)

                                tCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTg                                                         >  2:2725193/1‑93 (MQ=255)

                                tCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTg                                                         <  2:2993027/93‑1 (MQ=255)

                                  ggCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGAc                                                       <  2:4355914/93‑1 (MQ=255)

                                   gCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACg                                                      >  2:3543702/1‑93 (MQ=255)

                                     cTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTTACGGa                                                    >  2:4262095/1‑93 (MQ=255)

                                     cTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGa                                                    >  1:1982508/1‑93 (MQ=255)

                                     cTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGa                                                    >  1:4406811/1‑93 (MQ=255)

                                      tGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGTGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGCCGGaa                                                   >  2:4352525/1‑93 (MQ=255)

                                       ggAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGaaa                                                  <  1:2891783/93‑1 (MQ=255)

                                       ggAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGaaa                                                  >  2:1412541/1‑93 (MQ=255)

                                         aGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTTCCTCTGTTGATGTCATTGACGGaaaaa                                                <  1:3949136/93‑1 (MQ=255)

                                         aGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGaaaaa                                                <  1:1219169/93‑1 (MQ=255)

                                         aGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGaaaaa                                                >  1:1405845/1‑93 (MQ=255)

                                           cATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATc                                              <  2:4231962/93‑1 (MQ=255)

                                            aTGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGGTGAGATTGGCTCTGTTGATGTCATTGACGGAAAAATCg                                             >  2:3400289/1‑93 (MQ=255)

                                             tGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGt                                            >  2:3862578/1‑93 (MQ=255)

                                              gCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTc                                           <  1:1152539/93‑1 (MQ=255)

                                               ccATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCa                                          <  2:1123100/93‑1 (MQ=255)

                                               ccATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCa                                          <  1:4279415/93‑1 (MQ=255)

                                                cATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCaa                                         >  1:2396437/1‑93 (MQ=255)

                                                cATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCaa                                         >  1:2274299/1‑93 (MQ=255)

                                                cATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCaa                                         >  1:3248672/1‑93 (MQ=255)

                                                     tGCGCATTATCTATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCa                                    <  2:1972408/93‑1 (MQ=255)

                                                     tGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCa                                    >  1:2295069/1‑93 (MQ=255)

                                                     tGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCa                                    >  2:1826696/1‑93 (MQ=255)

                                                       cgcATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATg                                  <  2:3405493/93‑1 (MQ=255)

                                                       cgcATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATg                                  >  1:280932/1‑93 (MQ=255)

                                                       cgcATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATg                                  <  1:3955171/93‑1 (MQ=255)

                                                        gcATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGa                                 <  1:4266512/93‑1 (MQ=255)

                                                         cATTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAg                                <  1:3765629/93‑1 (MQ=255)

                                                          aTTATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGc                               >  1:2705122/1‑93 (MQ=255)

                                                           ttATCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGcc                              >  2:1554027/1‑93 (MQ=255)

                                                             aTCGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCAc                            >  1:558025/1‑93 (MQ=255)

                                                               cGATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCt                          >  1:1939829/1‑93 (MQ=255)

                                                                gATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTg                         >  1:1768350/1‑93 (MQ=255)

                                                                gATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTg                         >  2:2949911/1‑93 (MQ=255)

                                                                gATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTg                         <  1:172406/93‑1 (MQ=255)

                                                                gATACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTg                         >  1:3771899/1‑93 (MQ=255)

                                                                  tACAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGt                       >  1:3019628/1‑93 (MQ=255)

                                                                    cAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTgc                     <  2:249569/93‑1 (MQ=255)

                                                                    cAGAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTgc                     <  2:3905314/93‑1 (MQ=255)

                                                                      gAAACCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTgcgc                   >  1:2641018/1‑93 (MQ=255)

                                                                       aaaCCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGcc                  <  2:2262217/93‑1 (MQ=255)

                                                                       aaaCCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGcc                  <  2:4231000/93‑1 (MQ=255)

                                                                         aCCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGcccc                >  1:3097114/1‑93 (MQ=255)

                                                                         aCCTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGcccc                >  2:3882007/1‑93 (MQ=255)

                                                                           cTGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGa              >  2:1360492/1‑93 (MQ=255)

                                                                            tGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGAt             >  1:3025243/1‑93 (MQ=255)

                                                                            tGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGAt             >  1:2426494/1‑93 (MQ=255)

                                                                            tGCGGTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGAt             >  2:1266090/1‑93 (MQ=255)

                                                                                gTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTc         <  1:2642511/93‑1 (MQ=255)

                                                                                gTTTGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTc         >  1:1159000/1‑93 (MQ=255)

                                                                                  ttGCAGGGGGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCt       <  2:1486098/93‑1 (MQ=255)

                                                                                       gggggggaTCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCTATTAg  <  2:885421/93‑1 (MQ=255)

AAGGCGATATCATTCCATTACGCTAGGCTTTTTCGGCCTGGAGCATGCCATGTTGCGCATTATCGATACAGAAACCTGCGGTTTGCAGGGAGGGATCGTTGAGATTGCCTCTGTTGATGTCATTGACGGAAAAATCGTCAACCCCATGAGCCACCTGGTGCGCCCCGATCGTCCTATTAG  >  NC_000913/1924095‑1924274

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