BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,107,299	A→G	V289A (GTG→GCG)	ymdD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,107,299	0	A	G	99.2%	416.4 / NA	133	V289A (GTG→GCG)	ymdD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (79/53); minor base C (0/1); total (79/54)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

TGAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCACATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGT  >  NC_000913/1107209‑1107388

tgAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACcgca                                                                                         >  1:3090587/1‑93 (MQ=255)

tgAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACcgca                                                                                         >  2:3905947/1‑93 (MQ=255)

  aaCCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCAtt                                                                                       <  1:2561736/93‑1 (MQ=255)

  aaCCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCAtt                                                                                       <  1:383350/93‑1 (MQ=255)

  aaCCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCAtt                                                                                       <  2:628837/93‑1 (MQ=255)

     cAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCCAAGGAGAAGACCGCATTCAc                                                                                    >  2:4039066/1‑93 (MQ=255)

      agatagatAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCAcc                                                                                   >  1:756903/1‑93 (MQ=255)

       gatagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCa                                                                                  >  1:4168072/1‑93 (MQ=255)

        atagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCAt                                                                                 >  2:885798/1‑93 (MQ=255)

         tagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATc                                                                                >  2:809744/1‑93 (MQ=255)

         tagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATc                                                                                <  1:1697571/93‑1 (MQ=255)

         tagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATc                                                                                >  2:2280879/1‑93 (MQ=255)

         tagataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATc                                                                                >  2:4150554/1‑93 (MQ=255)

          agataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATcc                                                                               <  2:2165734/93‑1 (MQ=255)

          agataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATcc                                                                               >  1:4607236/1‑93 (MQ=255)

           gataAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCca                                                                              <  2:584619/93‑1 (MQ=255)

             taAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCcaca                                                                            <  1:3586882/93‑1 (MQ=255)

             taAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCcaca                                                                            >  1:2492160/1‑93 (MQ=255)

             taAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATAcaca                                                                            >  1:27599/1‑93 (MQ=255)

              aaaCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAg                                                                           >  1:585706/1‑93 (MQ=255)

               aaCAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGa                                                                          >  1:4126841/1‑93 (MQ=255)

                   gCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGa                                                                      >  2:4051704/1‑93 (MQ=255)

                     gATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAgg                                                                    >  2:1095145/1‑93 (MQ=255)

                      aTGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGCAGACCGCATTCACCATCCACAGACCGAGGa                                                                   >  2:2116003/1‑93 (MQ=255)

                      aTGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGa                                                                   <  2:2609454/93‑1 (MQ=255)

                        gCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGAcc                                                                 >  1:3613925/1‑93 (MQ=255)

                         cGTTAACAAAATAAGTCCCCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCa                                                                <  2:3798899/93‑1 (MQ=255)

                         cGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCa                                                                <  2:2692250/93‑1 (MQ=255)

                          gTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCAt                                                               <  2:3770114/93‑1 (MQ=255)

                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                              <  2:1288939/93‑1 (MQ=255)

                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                              <  1:3451499/93‑1 (MQ=255)

                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                              >  2:1683660/1‑93 (MQ=255)

                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                              >  1:4003391/1‑93 (MQ=255)

                           ttAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATg                                                              <  2:3778495/93‑1 (MQ=255)

                                aaaaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGAt                                                         <  2:3005790/93‑1 (MQ=255)

                                  aaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCa                                                       >  1:2771935/1‑93 (MQ=255)

                                  aaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCa                                                       >  1:1521460/1‑93 (MQ=255)

                                  aaTAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCa                                                       >  1:1945780/1‑93 (MQ=255)

                                    tAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCAcc                                                     >  2:892242/1‑93 (MQ=255)

                                      aGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGa                                                   >  2:2968071/1‑93 (MQ=255)

                                       gTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGAc                                                  >  2:2558088/1‑93 (MQ=255)

                                       gTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGAc                                                  >  2:2606594/1‑93 (MQ=255)

                                       gTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGAc                                                  >  1:1739586/1‑93 (MQ=255)

                                        gcccccGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACt                                                 <  1:2928748/90‑1 (MQ=255)

                                           cccGCGCTGACTGGAAGTTAAGCAAACGGAGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCgg                                              >  2:231270/1‑93 (MQ=255)

                                            ccGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGt                                             >  1:3201658/1‑93 (MQ=255)

                                              gcgcTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGttt                                           >  2:1009853/1‑93 (MQ=255)

                                               cgcTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTc                                          >  1:463404/1‑93 (MQ=255)

                                                 cTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGt                                        >  2:165086/1‑93 (MQ=255)

                                                  tGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTa                                       <  1:2157168/93‑1 (MQ=255)

                                                   gACTGGAAGTTAAGCAAACGGTGGCCCAAGGGGAAGAACGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTAc                                      >  2:558685/1‑93 (MQ=255)

                                                   gACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTAc                                      <  1:4356494/93‑1 (MQ=255)

                                                    aCTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACa                                     <  1:2210694/93‑1 (MQ=255)

                                                    aCTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACa                                     >  1:3745987/1‑93 (MQ=255)

                                                    aCTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACa                                     <  1:2475949/93‑1 (MQ=255)

                                                    aCTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACa                                     >  2:1166143/1‑93 (MQ=255)

                                                     cTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACAt                                    >  2:4426715/1‑93 (MQ=255)

                                                     cTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACAt                                    <  1:1471252/93‑1 (MQ=255)

                                                      tGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                   <  1:2487178/93‑1 (MQ=255)

                                                       ggAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATcc                                  <  1:3518064/93‑1 (MQ=255)

                                                       ggAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATcc                                  >  1:3949383/1‑93 (MQ=255)

                                                          aGTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAgg                               <  1:1009785/93‑1 (MQ=255)

                                                           gTTAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGc                              >  2:2824292/1‑93 (MQ=255)

                                                           cgtAAGCAAGCGGTGGCCAAAGGATAAGACCTCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGc                              <  2:2869681/91‑1 (MQ=255)

                                                             tAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCAt                            <  1:2411803/93‑1 (MQ=255)

                                                             tAAGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCAt                            <  2:2445746/93‑1 (MQ=255)

                                                              aaGCAAACGGTGCCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATc                           <  2:2041736/93‑1 (MQ=255)

                                                               aGCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATcc                          >  1:1963861/1‑92 (MQ=255)

                                                                gCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGc                         >  2:2910839/1‑93 (MQ=255)

                                                                gCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGc                         >  1:4213552/1‑93 (MQ=255)

                                                                gCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGc                         <  2:3002690/93‑1 (MQ=255)

                                                                gCAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGc                         >  1:2186676/1‑93 (MQ=255)

                                                                 cAAACGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGcc                        >  2:480163/1‑93 (MQ=255)

                                                                  aaaCGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCa                       >  2:2384548/1‑93 (MQ=255)

                                                                   aaCGGGGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGGTTCGTACATCCAGGCATCGCCAc                      >  2:1038082/1‑93 (MQ=255)

                                                                    aCGGTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCTAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACt                     >  1:3297960/1‑93 (MQ=255)

                                                                      ggTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGc                   <  1:4390381/93‑1 (MQ=255)

                                                                       gTGGCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGcc                  >  2:1635112/1‑93 (MQ=255)

                                                                          gCCAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATa               <  2:1280236/93‑1 (MQ=255)

                                                                            cAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgc             <  2:1204035/93‑1 (MQ=255)

                                                                            cAAAGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgc             <  1:377999/93‑1 (MQ=255)

                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg            >  1:959964/1‑93 (MQ=255)

                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg            >  1:1446024/1‑93 (MQ=255)

                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg            >  1:2428644/1‑93 (MQ=255)

                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg            >  1:267257/1‑93 (MQ=255)

                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAgcg            >  1:2878439/1‑93 (MQ=255)

                                                                             aaaGGAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGCGATCACCGACTCGGTTTCGTACATCCAGGCATAGCCACTGCCATAgcg            >  1:808969/1‑93 (MQ=255)

                                                                                ggAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTg         <  1:1066791/93‑1 (MQ=255)

                                                                                ggAGAAGACCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTg         <  2:4374208/93‑1 (MQ=255)

                                                                                   gaagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGtt      >  2:4221742/1‑93 (MQ=255)

                                                                                   gaagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGtt      >  1:1931043/1‑93 (MQ=255)

                                                                                    aagaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTg     <  2:865704/93‑1 (MQ=255)

                                                                                     agaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGa    <  1:4106879/93‑1 (MQ=255)

                                                                                     agaCCCCATTCACCATCCACAGACCGAGCACCATGGTTATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGa    <  1:258725/93‑1 (MQ=255)

                                                                                     agaCCACATTCACCGTCCACAAACCGAGGAACATGGTGATCACCGACTAGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGa    <  2:4686263/93‑1 (MQ=255)

                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg   >  2:4436149/1‑93 (MQ=255)

                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg   >  2:4664204/1‑93 (MQ=255)

                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg   <  1:4006754/93‑1 (MQ=255)

                                                                                      gaCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAg   <  1:1615831/93‑1 (MQ=255)

                                                                                       aCCGCATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGt  >  2:2373859/1‑93 (MQ=255)

TGAACCAGATAGATAAACAGCGATGCGTTAACAAAATAAGTCACCCGCGCTGACTGGAAGTTAAGCAAACGGTGGCCAAAGGAGAAGACCACATTCACCATCCACAGACCGAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGT  >  NC_000913/1107209‑1107388

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