BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,949,436	A→G	L6L (CTA→CTG)	yecE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,436	0	A	G	75.7%	1502.4 / inf	811	L6L (CTA→CTG)	yecE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (125/72); new base G (342/272); total (467/344)
Fisher's exact test for biased strand distribution p-value = 5.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 882 total aligned reads displayed.

GCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCA  >  NC_000913/1949344‑1949527

gcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTa                                                                                             >  2:3279564/1‑93 (MQ=255)

 cacaATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTAc                                                                                            <  1:4976843/93‑1 (MQ=255)

   caATAACAGTATCAACCATATCTACCCGCGCCTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCg                                                                                          >  1:2995530/1‑93 (MQ=255)

   caATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCg                                                                                          >  2:2467904/1‑93 (MQ=255)

    aaTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGc                                                                                         >  1:1387053/1‑93 (MQ=255)

    aaTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGc                                                                                         >  2:6044903/1‑93 (MQ=255)

     aTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCa                                                                                        >  1:5740027/1‑93 (MQ=255)

      tAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCaa                                                                                       <  2:5561797/93‑1 (MQ=255)

      tAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCaa                                                                                       >  2:1278741/1‑93 (MQ=255)

         cAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATgg                                                                                    <  1:3769473/93‑1 (MQ=255)

         cAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATgg                                                                                    <  2:43951/93‑1 (MQ=255)

                 atcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATcc                                                                            <  1:2567369/91‑1 (MQ=255)

                 aCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATcc                                                                            >  2:2099636/1‑93 (MQ=255)

                 aCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATcc                                                                            <  1:4267920/93‑1 (MQ=255)

                  tcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCt                                                                           >  2:5515653/2‑93 (MQ=255)

                   cATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTa                                                                          >  2:1933219/1‑93 (MQ=255)

                   cATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTa                                                                          <  1:3397307/93‑1 (MQ=255)

                   cATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTa                                                                          >  1:1089758/1‑93 (MQ=255)

                    atatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaa                                                                         >  2:841923/1‑93 (MQ=255)

                    atatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTaa                                                                         >  1:365480/1‑93 (MQ=255)

                    atatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTaa                                                                         >  1:5227540/1‑93 (MQ=255)

                     tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                                                        <  1:2314594/93‑1 (MQ=255)

                     tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                                                        >  2:4937792/1‑93 (MQ=255)

                     tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                                                        <  2:3590657/93‑1 (MQ=255)

                     tatCTACCCGCGCATCGCCCGTGCGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                                                        >  2:3731454/1‑93 (MQ=255)

                      atCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAAt                                                                       <  2:1357093/93‑1 (MQ=255)

                      atCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAAt                                                                       <  2:5244378/93‑1 (MQ=255)

                      atCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAAt                                                                       >  1:6122277/1‑93 (MQ=255)

                       tCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATg                                                                      >  2:3061972/1‑93 (MQ=255)

                         tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATggg                                                                    >  2:169576/1‑93 (MQ=255)

                         tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATggg                                                                    <  1:4073668/93‑1 (MQ=255)

                         tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAGTGGTCGCATCCTAAAgggg                                                                    >  2:5302391/1‑93 (MQ=255)

                         tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATggg                                                                    <  2:688554/93‑1 (MQ=255)

                         tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATggg                                                                    >  2:4527604/1‑93 (MQ=255)

                          aCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGt                                                                   >  2:404622/1‑93 (MQ=255)

                          aCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGt                                                                   <  1:4858424/93‑1 (MQ=255)

                           cccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTg                                                                  <  2:3489987/93‑1 (MQ=255)

                           cccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTg                                                                  >  1:3545024/1‑93 (MQ=255)

                           cccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTg                                                                  <  1:1264650/93‑1 (MQ=255)

                            ccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGc                                                                 <  2:3403137/93‑1 (MQ=255)

                                gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                                             >  1:477870/1‑93 (MQ=255)

                                gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATTGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                                             <  1:5408806/93‑1 (MQ=255)

                                gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                                             >  2:175179/1‑93 (MQ=255)

                                  aTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgg                                                           >  1:2089343/1‑93 (MQ=255)

                                   tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                                          >  1:1007517/1‑93 (MQ=255)

                                   tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                                          >  1:4383807/1‑93 (MQ=255)

                                   tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                                          >  1:4185222/1‑93 (MQ=255)

                                   tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                                          >  2:1576704/1‑93 (MQ=255)

                                    cGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgggt                                                         <  2:6098438/93‑2 (MQ=255)

                                     gCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgggta                                                        <  1:2666655/93‑3 (MQ=255)

                                     gCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgggta                                                        <  2:1849759/93‑3 (MQ=255)

                                      cccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTAt                                                       <  1:3489242/93‑1 (MQ=255)

                                      cccGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGAt                                                       >  2:5050194/1‑93 (MQ=255)

                                       ccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATc                                                      >  2:2200978/1‑93 (MQ=255)

                                        cGTGTGCGTCGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATcc                                                     >  1:3847206/1‑92 (MQ=255)

                                        cGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCa                                                     >  2:1490576/1‑93 (MQ=255)

                                         gtgtGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAc                                                    >  1:5862253/1‑93 (MQ=255)

                                          tgtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAcc                                                   <  1:98900/93‑1 (MQ=255)

                                           gtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCa                                                  <  1:280089/93‑1 (MQ=255)

                                           gtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCa                                                  <  2:2307121/93‑1 (MQ=255)

                                           gtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCa                                                  <  1:1946771/93‑1 (MQ=255)

                                           gtgCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCa                                                  >  2:5910595/1‑93 (MQ=255)

                                            tgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAg                                                 >  2:5406953/1‑93 (MQ=255)

                                              cttAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGcc                                               <  2:4865230/91‑1 (MQ=255)

                                                 aGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTg                                            <  2:3680602/93‑1 (MQ=255)

                                                  tcGTGTAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGa                                           <  1:564183/92‑1 (MQ=255)

                                                   cGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGaa                                          <  1:2910244/93‑1 (MQ=255)

                                                     tGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAaga                                        >  1:3150245/1‑93 (MQ=255)

                                                     tGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAaga                                        >  2:3140506/1‑93 (MQ=255)

                                                     tGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAaga                                        >  2:1233117/1‑93 (MQ=255)

                                                     tGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTCAATGGGTGCGGTTGGGGATCACCAGCCTTGAaga                                        >  2:2977387/1‑93 (MQ=255)

                                                       gAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGt                                      >  1:5057973/1‑93 (MQ=255)

                                                        aaGAGATCCTTCACTCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa                                     <  2:4766564/93‑1 (MQ=255)

                                                        aaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa                                     <  1:2501097/93‑1 (MQ=255)

                                                        aaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa                                     <  2:269507/93‑1 (MQ=255)

                                                        aaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa                                     <  2:5646625/93‑1 (MQ=255)

                                                        aaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa                                     <  1:5343183/93‑1 (MQ=255)

                                                        aaGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTa                                     <  1:2391079/93‑1 (MQ=255)

                                                         agagATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTAt                                    <  1:911419/93‑1 (MQ=255)

                                                         agagATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTAt                                    <  1:4544738/93‑1 (MQ=255)

                                                          gagaTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATg                                   >  2:2583323/1‑93 (MQ=255)

                                                           agaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGc                                  >  1:1835162/1‑93 (MQ=255)

                                                           agaTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGc                                  >  1:3926228/1‑93 (MQ=255)

                                                                   caCGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACt                          <  1:5548709/92‑1 (MQ=255)

                                                                   aaCGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACt                          >  2:4204538/1‑93 (MQ=255)

                                                                     cgcgTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACttt                        >  2:3235824/1‑93 (MQ=255)

                                                                      gcgTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTa                       >  1:2208997/1‑93 (MQ=255)

                                                                      gcgTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTa                       >  1:1492519/1‑93 (MQ=255)

                                                                       ctttATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTGGGGTACCACCAGCCTTGAAGAGTATGCCCGCCACTTTaa                      <  2:1156469/91‑1 (MQ=255)

                                                                       ctttATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAAGGGGGGCGGTTGGGTACCCCCAGCCTTGAAGAGTATGCCCGCCACTTTaa                      <  2:2824541/91‑1 (MQ=255)

                                                                       cgTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTaa                      <  1:6035371/93‑1 (MQ=255)

                                                                        gTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAAc                     >  2:5130650/1‑93 (MQ=255)

                                                                         ttATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACt                    >  1:3265914/1‑93 (MQ=255)

                                                                         ttATGATTTACATCGGTCTGCCGAAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAAcg                    >  1:1601697/1‑92 (MQ=255)

                                                                           aTGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGc                  >  2:750226/1‑93 (MQ=255)

                                                                            tGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCg                 >  1:1738404/1‑93 (MQ=255)

                                                                               tttACATCGTTCTGCCGCAGTGGTCGCATCCTAAATGGGTGGGGTTGGGTACCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTga              <  2:5710337/93‑2 (MQ=255)

                                                                                 tACATCGGTCTACCGCAATGGTCGCATCCTAAAGGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAg            <  2:4622827/93‑1 (MQ=255)

                                                                                  aCATCGGTCTCCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTgacgc           <  2:3857591/93‑5 (MQ=255)

                                                                                     tCGGTCTACCGCAATGGTCGCATCCAAAGGGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCa        <  1:74803/93‑1 (MQ=255)

GCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCA  >  NC_000913/1949344‑1949527

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