BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,855,431	A→G	D121D (GAT→GAC)	b1773 ←	putative aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,855,431	0	A	G	99.1%	698.1 / ‑4.8	226	D121D (GAT→GAC)	b1773	putative aldolase
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (86/138); total (86/140)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

GTGATAGCCCGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAATCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCATGAGTCCAGG  >  NC_000913/1855341‑1855519

 tGATAGCCCGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCg                                                                                       >  1:2652186/1‑93 (MQ=255)

  gATAGCCCGCCAGCGCTTCCTCAGAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGa                                                                                      <  1:5622343/93‑1 (MQ=255)

   aTAGCCCGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGAc                                                                                     <  1:1173100/93‑1 (MQ=255)

   aTAGCCCGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGAc                                                                                     <  2:4016390/93‑1 (MQ=255)

    tAGCCCGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAGAAAGTCGAcc                                                                                    >  1:2740380/1‑92 (MQ=255)

      gCCCGCCGGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATGCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGc                                                                                  <  1:6857272/93‑1 (MQ=255)

      gCCCGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGc                                                                                  >  2:6427639/1‑93 (MQ=255)

        ccGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCct                                                                                >  1:2176894/1‑93 (MQ=255)

        ccGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCct                                                                                <  1:4782446/93‑1 (MQ=255)

         cGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCctc                                                                               <  2:5765614/93‑1 (MQ=255)

         cGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCctc                                                                               >  1:1046448/1‑93 (MQ=255)

         cGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCctc                                                                               >  2:4431495/1‑93 (MQ=255)

         cGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCctc                                                                               <  1:3944168/93‑1 (MQ=255)

         cGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCctc                                                                               <  2:797162/93‑1 (MQ=255)

          gCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCctct                                                                              <  1:3196822/93‑1 (MQ=255)

            cAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCttt                                                                            <  1:3250591/93‑1 (MQ=255)

            cAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCttt                                                                            <  1:5910928/93‑1 (MQ=255)

             aGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTg                                                                           >  2:4486047/1‑93 (MQ=255)

             aGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTg                                                                           <  2:2389804/93‑1 (MQ=255)

             aGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTg                                                                           >  2:5572124/1‑93 (MQ=255)

             aGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTg                                                                           >  2:2158416/1‑93 (MQ=255)

                 cTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTa                                                                       <  1:6733329/93‑1 (MQ=255)

                  ttCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTaa                                                                      <  1:7395175/93‑1 (MQ=255)

                  ttCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTaa                                                                      >  1:2285964/1‑93 (MQ=255)

                  ttCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTaa                                                                      <  2:6328556/93‑1 (MQ=255)

                    ccTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAAc                                                                    <  1:4242879/93‑1 (MQ=255)

                    cataaTAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAAc                                                                    <  2:2010268/89‑1 (MQ=255)

                     cTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACg                                                                   >  1:6363717/1‑93 (MQ=255)

                      tCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGa                                                                  >  1:95938/1‑93 (MQ=255)

                      tCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGa                                                                  <  2:2934544/93‑1 (MQ=255)

                        aTAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAat                                                                >  1:3004851/1‑93 (MQ=255)

                         tAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAata                                                               <  2:876773/93‑1 (MQ=255)

                         tAAACGGTTTCATTACCAACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAata                                                               <  2:3117923/93‑1 (MQ=255)

                          aaaCGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAatat                                                              <  2:3285155/93‑1 (MQ=255)

                          aaaCGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAatat                                                              <  2:6682844/93‑1 (MQ=255)

                           aaCGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATAtt                                                             >  1:7497419/1‑93 (MQ=255)

                           aaCGGTTCGATTACCGGCATGCCCTAACTCCGCGTCTACCGGAAGACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATAtt                                                             <  1:266013/93‑1 (MQ=255)

                             cGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATAtttt                                                           <  1:4603040/93‑1 (MQ=255)

                              gggttCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTttc                                                          <  1:5156597/90‑1 (MQ=255)

                              ggTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTttc                                                          <  2:1213124/93‑1 (MQ=255)

                               gTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTttct                                                         >  1:1881914/1‑93 (MQ=255)

                                tttCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTttctt                                                        <  1:7571660/93‑1 (MQ=255)

                                 ttCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTttcttc                                                       <  1:7756113/93‑1 (MQ=255)

                                 ttCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTttcttc                                                       <  1:4949307/93‑1 (MQ=255)

                                 ttCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTttcttc                                                       <  1:1476176/93‑1 (MQ=255)

                                   cATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGa                                                     >  1:3343929/1‑93 (MQ=255)

                                    aTTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGaa                                                    >  2:7171620/1‑93 (MQ=255)

                                      tACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAAt                                                  >  2:5715982/1‑93 (MQ=255)

                                        ccGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCa                                                <  2:6215872/93‑1 (MQ=255)

                                        ccGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCa                                                >  2:7537330/1‑93 (MQ=255)

                                          gACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCaaa                                              >  2:4551404/1‑93 (MQ=255)

                                          gACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCaaa                                              <  2:1112502/93‑1 (MQ=255)

                                          gACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCaaa                                              >  1:3471951/1‑93 (MQ=255)

                                           ccctGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAAc                                             <  2:6031687/90‑1 (MQ=255)

                                           aCATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAAc                                             <  2:3673919/93‑1 (MQ=255)

                                           aCATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAAc                                             <  1:6551576/93‑1 (MQ=255)

                                             aTGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAAcgc                                           >  2:904777/1‑93 (MQ=255)

                                             aTGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAAcgc                                           >  2:2697527/1‑93 (MQ=255)

                                               ccccTAATTCCGCCTCTACCGGAATCCCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCgg                                         <  2:2856677/92‑1 (MQ=255)

                                                cccTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGa                                        <  2:6168798/93‑1 (MQ=255)

                                                cccTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGa                                        >  1:6160318/1‑93 (MQ=255)

                                                 ccTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGaa                                       <  1:7290921/93‑1 (MQ=255)

                                                  cTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAg                                      >  1:6182502/1‑93 (MQ=255)

                                                     aTTCCGCCTCTACCGGAATAACCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACACGCAAGCAt                                   >  2:5505898/1‑93 (MQ=255)

                                                      ttCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATc                                  >  1:3134709/1‑93 (MQ=255)

                                                       tCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCa                                 <  1:3453551/93‑1 (MQ=255)

                                                              ctACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTa                          >  2:1649067/1‑93 (MQ=255)

                                                              ctACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTa                          >  2:7759192/1‑93 (MQ=255)

                                                              ctACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTa                          >  2:4239227/1‑93 (MQ=255)

                                                              ctACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTa                          >  2:5059832/1‑93 (MQ=255)

                                                               tACCGGAATACCCAGCGGATGGTAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTaa                         >  2:7493244/1‑93 (MQ=255)

                                                               tACCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTaa                         >  2:1961461/1‑93 (MQ=255)

                                                                tccGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAAt                        <  1:643772/92‑1 (MQ=255)

                                                                aCCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAAt                        <  1:6238153/93‑1 (MQ=255)

                                                                aCCGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAAt                        >  2:5361672/1‑93 (MQ=255)

                                                                  cGGAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGa                      >  2:4777147/1‑93 (MQ=255)

                                                                   ggAATACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGaa                     <  2:5875272/93‑1 (MQ=255)

                                                                      aTACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTc                  <  1:778897/93‑1 (MQ=255)

                                                                      aTACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTc                  <  1:5401119/93‑1 (MQ=255)

                                                                       tACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCa                 <  2:728801/93‑1 (MQ=255)

                                                                       tACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCa                 >  1:5235105/1‑93 (MQ=255)

                                                                       tACCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCa                 <  1:5146387/93‑1 (MQ=255)

                                                                        aCCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCAt                >  1:471632/1‑93 (MQ=255)

                                                                        aCCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCAt                >  1:6936650/1‑93 (MQ=255)

                                                                        aCCCAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCAt                >  1:2832531/1‑93 (MQ=255)

                                                                         cccAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATa               >  1:1510003/1‑93 (MQ=255)

                                                                         cccAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATa               <  2:5035055/93‑1 (MQ=255)

                                                                           cAGCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATAcc             >  2:7303355/1‑93 (MQ=255)

                                                                            aGCGCATGTAGGAAGTCGACAGCCACTTTGGGTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATAGCa            <  2:3809002/93‑1 (MQ=255)

                                                                             gCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCAt           <  1:3988606/93‑1 (MQ=255)

                                                                             gCGGATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCAt           <  2:4891917/93‑1 (MQ=255)

                                                                             gCGGATCGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCAt           <  2:3411646/93‑1 (MQ=255)

                                                                              gggATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCATg          <  1:4368143/92‑1 (MQ=255)

                                                                               ggATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCATGa         <  1:7830225/93‑1 (MQ=255)

                                                                               ggATGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCATGa         <  2:1410796/93‑1 (MQ=255)

                                                                                  tGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCATGAGTc      <  2:4802128/93‑1 (MQ=255)

                                                                                  tGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCATGAGTc      >  2:4619488/1‑93 (MQ=255)

                                                                                  tGGAAGAAGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCATGAGTc      >  1:5675098/1‑93 (MQ=255)

                                                                                    gaagaaGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCATGAGTCCa    <  2:1235627/93‑1 (MQ=255)

                                                                                    gaagaaGTCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCATGAGTCCa    <  2:6889758/93‑1 (MQ=255)

GTGATAGCCCGCCAGCGCTTCCTCATAAACGGTTTCATTACCGACATGCCCTAATTCCGCCTCTACCGGAATACCCAGCGGATGGAAGAAATCGACAGCCTCTTTGGTTAAACGAATATTTTCTTCGAAATCAAACGCGGAAGCATCACGCATTAATGAATTCATACCATGAGTCCAGG  >  NC_000913/1855341‑1855519

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