BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,679,869	A→G	S62S (AGT→AGC)	ydgC ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,679,869	0	A	G	99.6%	959.8 / ‑4.8	274	S62S (AGT→AGC)	ydgC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (155/118); total (155/119)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 283 total aligned reads displayed.

AGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGG  >  NC_000913/1679779‑1679958

 gAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGCAAACAAAATAGGGAATAATCGACCACATGCTa                                                                                        >  1:6148157/1‑93 (MQ=255)

   aCCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTaaa                                                                                      <  1:533895/93‑1 (MQ=255)

   aCCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTaaa                                                                                      >  1:4218650/1‑93 (MQ=255)

    ccAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTaaaa                                                                                     <  1:7066882/93‑1 (MQ=255)

     caacaaAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTaaaaa                                                                                    <  2:2278242/93‑1 (MQ=255)

      aacaaAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAat                                                                                   <  2:2907987/93‑1 (MQ=255)

      aacaaAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAat                                                                                   >  2:706048/1‑93 (MQ=255)

       acaaAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAatg                                                                                  >  2:5226550/1‑93 (MQ=255)

         aaaaGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAatgat                                                                                >  1:721251/1‑93 (MQ=255)

         aaaaGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAatgat                                                                                >  1:5154051/1‑93 (MQ=255)

          aaaGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAatgatg                                                                               >  2:4201100/1‑93 (MQ=255)

          aaaGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAatgatg                                                                               <  2:7865123/93‑1 (MQ=255)

          aaaGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCAAAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAatgatg                                                                               <  1:4562549/93‑1 (MQ=255)

           aaGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATgg                                                                              <  1:1291366/93‑1 (MQ=255)

           aaGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATgg                                                                              <  2:2582279/93‑1 (MQ=255)

            aGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGt                                                                             >  2:8708873/1‑93 (MQ=255)

            aGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGt                                                                             <  1:7392897/93‑1 (MQ=255)

              cggcggGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTg                                                                           <  2:4684575/93‑1 (MQ=255)

               ggcggGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTgc                                                                          >  1:2824353/1‑93 (MQ=255)

               ggcggGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTgc                                                                          <  2:6581155/93‑1 (MQ=255)

                gcggGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTgcg                                                                         <  2:3337868/93‑1 (MQ=255)

                gcggGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTgcg                                                                         <  2:6747012/93‑1 (MQ=255)

                   ggCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGcgcga                                                                      >  1:2664917/1‑92 (MQ=255)

                   ggCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGt                                                                      >  2:444234/1‑93 (MQ=255)

                   ggCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGt                                                                      >  1:3813026/1‑93 (MQ=255)

                   ggCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGt                                                                      >  2:4051795/1‑93 (MQ=255)

                    gCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTa                                                                     <  2:1089916/93‑1 (MQ=255)

                      aaGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAg                                                                   <  1:826443/93‑1 (MQ=255)

                      aaGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAg                                                                   >  2:1041108/1‑93 (MQ=255)

                       aGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAgg                                                                  <  1:8548014/93‑1 (MQ=255)

                       aGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAgg                                                                  >  2:3978022/1‑93 (MQ=255)

                        gcgcATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGc                                                                 >  2:3068387/1‑93 (MQ=255)

                        gcgcATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAAGAGGGAGTAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGc                                                                 <  1:6386089/93‑1 (MQ=255)

                          gcATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCtt                                                               >  2:1711571/1‑93 (MQ=255)

                          gcATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCtt                                                               >  2:8300602/1‑93 (MQ=255)

                          gcATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCtt                                                               <  1:7662470/93‑1 (MQ=255)

                          gcATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCtt                                                               >  2:5033664/1‑93 (MQ=255)

                            atcatTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCa                                                             >  2:8495019/1‑93 (MQ=255)

                            atcatTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCa                                                             >  1:4396026/1‑93 (MQ=255)

                             tcatTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCaa                                                            >  2:5496066/1‑93 (MQ=255)

                             tcatTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCaa                                                            >  2:2111668/1‑93 (MQ=255)

                             tcatTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCaa                                                            >  1:2313709/1‑93 (MQ=255)

                               cttCGGGTGAAATACGACAGCTACACCAGGTAAACAAAATAGGGAATAATCGCCCACATGCTAAAAATGCTGGTTGCGCGTAAGGCTTAAATg                                                          <  2:3176098/92‑1 (MQ=255)

                               atTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATg                                                          >  1:4630924/1‑93 (MQ=255)

                                  ccGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCcg                                                       <  1:6707934/93‑1 (MQ=255)

                                   cGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCcgc                                                      >  1:2027825/1‑93 (MQ=255)

                                   cGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCcgc                                                      <  2:1866890/93‑1 (MQ=255)

                                   cGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCcgc                                                      >  2:1274328/1‑93 (MQ=255)

                                   cGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCcgc                                                      >  1:8439125/1‑93 (MQ=255)

                                      tGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGtt                                                   <  2:8263133/93‑1 (MQ=255)

                                      tGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGtt                                                   <  1:1990104/93‑1 (MQ=255)

                                      tGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGtt                                                   <  1:1454563/93‑1 (MQ=255)

                                       gAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTc                                                  <  1:4657051/93‑1 (MQ=255)

                                       gAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTc                                                  <  2:346239/93‑1 (MQ=255)

                                         aaTACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGc                                                <  1:2717550/93‑1 (MQ=255)

                                           tACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTg                                              >  2:1389406/1‑93 (MQ=255)

                                            aCCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTgg                                             <  1:3274907/93‑1 (MQ=255)

                                              cacaGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGca                                           >  1:4168367/1‑93 (MQ=255)

                                               acaGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGcaa                                          >  1:5591749/1‑93 (MQ=255)

                                                caGCGACGCCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGcaac                                         <  2:6261172/93‑1 (MQ=255)

                                                   cGACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAat                                      <  2:5897725/93‑1 (MQ=255)

                                                    gACACCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAata                                     <  1:2885965/93‑1 (MQ=255)

                                                     acacCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAatat                                    <  2:2065251/93‑1 (MQ=255)

                                                     acacCAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAatat                                    >  2:7326871/1‑93 (MQ=255)

                                                       acCAGGTAAACAAAATAGGGAATAATCAACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATaa                                  >  2:6716142/1‑93 (MQ=255)

                                                        ccAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAAt                                 >  1:6843419/1‑93 (MQ=255)

                                                         cAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATg                                <  2:8771985/93‑1 (MQ=255)

                                                         cAGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATg                                >  1:8774460/1‑93 (MQ=255)

                                                          aGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATgc                               <  1:2117906/93‑1 (MQ=255)

                                                          aGGTAAACAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGT‑AGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATgcg                              >  2:8106712/1‑93 (MQ=255)

                                                              aaaCAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGa                           >  1:1013359/1‑93 (MQ=255)

                                                               aaCAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGAt                          <  2:2810000/93‑1 (MQ=255)

                                                               aaCAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGAt                          >  1:3061059/1‑93 (MQ=255)

                                                               aaCAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGAt                          >  1:3383143/1‑93 (MQ=255)

                                                               aaCAAAATAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGAt                          >  1:4612399/1‑93 (MQ=255)

                                                                  aaaaTAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAg                       >  2:7490062/1‑93 (MQ=255)

                                                                  aaaaTAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAg                       <  1:5424312/93‑1 (MQ=255)

                                                                   aaaTAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAgc                      <  1:4775233/93‑1 (MQ=255)

                                                                    aaTAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAgcg                     >  1:6084332/1‑93 (MQ=255)

                                                                      tAGGGAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCa                   >  1:7592147/1‑93 (MQ=255)

                                                                        gggAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCaaa                 >  2:436296/1‑93 (MQ=255)

                                                                         ggAATAATCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAg                <  2:694409/93‑1 (MQ=255)

                                                                          gAATAATCGACCACATGCTAAAAATGATGGTTGCACGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAgg               <  1:5646843/93‑1 (MQ=255)

                                                                           aataatCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGt              >  2:3712865/1‑93 (MQ=255)

                                                                           aataatCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGt              <  2:3909981/93‑1 (MQ=255)

                                                                             taatCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCg            >  1:6526937/1‑93 (MQ=255)

                                                                              aatCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCgg           >  1:5934418/1‑93 (MQ=255)

                                                                              aatCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCgg           >  1:4864988/1‑93 (MQ=255)

                                                                               atCGACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCggg          >  2:8408433/1‑93 (MQ=255)

                                                                                  gACCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGaaa       <  1:4054247/93‑1 (MQ=255)

                                                                                   aCCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGaaaa      >  2:2962735/1‑93 (MQ=255)

                                                                                   aCCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGaaaa      >  2:439693/1‑93 (MQ=255)

                                                                                   aCCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGaaaa      >  2:610730/1‑93 (MQ=255)

                                                                                   aCCACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGaaaa      <  1:8028793/93‑1 (MQ=255)

                                                                                    ccACATGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAg     <  1:8138088/93‑1 (MQ=255)

                                                                                     cacaTGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGt    >  2:6549219/1‑93 (MQ=255)

                                                                                     cacaTGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGt    <  2:1171899/93‑1 (MQ=255)

                                                                                     cacaTGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGt    <  1:1624076/93‑1 (MQ=255)

                                                                                     cacaTGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGt    <  1:8184323/93‑1 (MQ=255)

                                                                                      acaTGCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTg   <  1:8489471/93‑1 (MQ=255)

AGAACCAACAAAAGCGGCGGGCAAGCGCATCATTCCGGTGAAATACCACAGCGACACCAGGTAAACAAAATAGGGAATAATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGG  >  NC_000913/1679779‑1679958

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