BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	37,332	A→G	C165R (TGC→CGC)	caiC ←	putative crotonobetain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	37,332	0	A	G	99.7%	917.3 / ‑5.1	283	C165R (TGC→CGC)	caiC	putative crotonobetain
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (126/156); total (126/157)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 293 total aligned reads displayed.

GGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCC  >  NC_000913/37242‑37421

    aaTCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGt                                                                                     <  1:6079920/93‑1 (MQ=255)

    aaTCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGt                                                                                     <  2:3006083/93‑1 (MQ=255)

     aTCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTg                                                                                    <  1:6028348/93‑1 (MQ=255)

      tCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTgg                                                                                   >  2:4702197/1‑93 (MQ=255)

       caccacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGc                                                                                  >  1:2750143/1‑93 (MQ=255)

        accacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCa                                                                                 >  1:4686236/1‑93 (MQ=255)

        accacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCa                                                                                 >  1:1230758/1‑93 (MQ=255)

          cgcacCTTTCGCTCGGGAGGTGGTGCCGGAGTTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAgg                                                                               <  1:7267/91‑1 (MQ=255)

          cacaCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAgg                                                                               >  2:1277698/1‑93 (MQ=255)

           ccacCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgt                                                                              <  1:1790345/92‑1 (MQ=255)

           ccacCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgt                                                                              <  1:5848537/92‑1 (MQ=255)

           acacCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgt                                                                              >  1:4296782/1‑93 (MQ=255)

           aacccTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgt                                                                              <  1:3022183/90‑1 (MQ=255)

            cacCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgtt                                                                             >  2:571759/1‑93 (MQ=255)

            cacCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgtt                                                                             >  1:930668/1‑93 (MQ=255)

            cacCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgtt                                                                             <  1:578577/93‑1 (MQ=255)

             acCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgttg                                                                            >  2:748987/1‑93 (MQ=255)

              ccTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgttgt                                                                           <  1:298411/93‑1 (MQ=255)

              ccTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgttgt                                                                           <  1:26400/93‑1 (MQ=255)

              ccTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgttgt                                                                           <  2:6501014/93‑1 (MQ=255)

               cTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgttgtt                                                                          >  1:592338/1‑93 (MQ=255)

               cTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGgctgtt                                                                          >  1:3109989/1‑93 (MQ=255)

                 ttCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGa                                                                        <  2:4258539/93‑1 (MQ=255)

                   cggtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGAtt                                                                      <  1:559447/93‑1 (MQ=255)

                    ggtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGAttt                                                                     >  1:5077884/1‑93 (MQ=255)

                     gtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGAtttt                                                                    <  2:880767/93‑1 (MQ=255)

                     gtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGAGGCAGGTTGTTGAtttt                                                                    >  1:4082865/1‑93 (MQ=255)

                         ggAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAg                                                                >  2:2858361/1‑93 (MQ=255)

                         ggAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAg                                                                <  2:1921866/93‑1 (MQ=255)

                         ggAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAg                                                                <  2:3477316/93‑1 (MQ=255)

                                gtgCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTa                                                         <  2:1720254/93‑1 (MQ=255)

                                 tcccGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTaa                                                        <  2:447334/91‑1 (MQ=255)

                                  gCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTaaa                                                       >  2:5557567/1‑93 (MQ=255)

                                   ccGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAAc                                                      <  2:5361055/93‑1 (MQ=255)

                                     ggAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGATGATTTTTCAGTTGAGTAAACGa                                                    >  2:1664305/1‑93 (MQ=255)

                                      gAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGaa                                                   <  2:2422434/93‑1 (MQ=255)

                                      gAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGaa                                                   <  2:2044314/93‑1 (MQ=255)

                                      gAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGaa                                                   <  2:2561063/93‑1 (MQ=255)

                                      gAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGaa                                                   <  2:1374062/93‑1 (MQ=255)

                                       aGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAAc                                                  <  2:2248883/93‑1 (MQ=255)

                                         gTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTc                                                >  1:5400669/1‑93 (MQ=255)

                                         gTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTc                                                <  1:5240983/93‑1 (MQ=255)

                                          tGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCa                                               <  2:6280529/93‑1 (MQ=255)

                                           gAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCAc                                              >  2:3670787/1‑93 (MQ=255)

                                           gAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCAc                                              <  1:4105546/93‑1 (MQ=255)

                                           gAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCAc                                              <  2:6192411/93‑1 (MQ=255)

                                            aaGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATGGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACg                                             >  1:5469319/1‑93 (MQ=255)

                                            aaGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACg                                             >  1:6130046/1‑93 (MQ=255)

                                            aaGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGCGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACg                                             >  2:5976885/1‑93 (MQ=255)

                                              gagaATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGcc                                           >  1:362209/1‑93 (MQ=255)

                                              gagaATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGcc                                           >  1:2875022/1‑93 (MQ=255)

                                              gagaATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGcc                                           >  1:1577344/1‑93 (MQ=255)

                                                gaATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCcat                                         <  1:1023667/93‑1 (MQ=255)

                                                 aaTGTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCcatc                                        <  2:4511314/93‑1 (MQ=255)

                                                  aTTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCcatca                                       <  2:3312307/93‑1 (MQ=255)

                                                   tttCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCcatcat                                      <  2:4569918/93‑1 (MQ=255)

                                                     tCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCAtcatca                                    <  1:2345541/93‑1 (MQ=255)

                                                     tCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCAtcatca                                    <  1:3391498/93‑1 (MQ=255)

                                                      ccgccgTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAg                                   <  2:4826480/93‑1 (MQ=255)

                                                         ccgTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCgg                                <  1:748643/93‑1 (MQ=255)

                                                         ccgTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCgg                                >  2:2652284/1‑93 (MQ=255)

                                                         ccgTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCgg                                <  1:6221109/93‑1 (MQ=255)

                                                          cgTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCggg                               <  1:3460228/93‑1 (MQ=255)

                                                          cgTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCggg                               >  1:2437265/1‑93 (MQ=255)

                                                           gTATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGa                              <  2:5826850/93‑1 (MQ=255)

                                                            tATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGaa                             <  1:3955857/93‑1 (MQ=255)

                                                            tATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGaa                             <  1:4612710/93‑1 (MQ=255)

                                                            tATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGaa                             >  1:1999212/1‑93 (MQ=255)

                                                            tATCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGaa                             <  2:4036739/93‑1 (MQ=255)

                                                              tCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGt                           <  2:96522/93‑1 (MQ=255)

                                                              tCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGt                           <  2:380810/93‑1 (MQ=255)

                                                              tCGTCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGt                           >  1:5181097/1‑93 (MQ=255)

                                                                 tCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGcc                        >  1:3165072/1‑93 (MQ=255)

                                                                 tCAGTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGcc                        <  2:2131748/93‑1 (MQ=255)

                                                                    gTCGATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCcaca                     <  1:1679416/93‑1 (MQ=255)

                                                                       gATAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCt                  >  1:6605186/1‑93 (MQ=255)

                                                                        aTAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTg                 >  2:5622904/1‑93 (MQ=255)

                                                                         tAGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGt                <  1:1409576/93‑1 (MQ=255)

                                                                          aGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTc               >  1:1107966/1‑93 (MQ=255)

                                                                          aGCGGCGGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTc               >  1:4819128/1‑93 (MQ=255)

                                                                           gcggcgGTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCa              >  2:3572319/1‑93 (MQ=255)

                                                                           gcggcgGTGCATAGCGCAAGGGGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCa              >  1:811242/1‑93 (MQ=255)

                                                                            cggcggTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAg             <  2:3752162/93‑1 (MQ=255)

                                                                            cggcggTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAg             >  1:3491918/1‑93 (MQ=255)

                                                                            cggcggTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAg             <  2:70899/93‑1 (MQ=255)

                                                                              gcggTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGc           >  2:6415886/1‑93 (MQ=255)

                                                                              gcggTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGc           <  1:4433599/93‑1 (MQ=255)

                                                                               cggTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCa          <  2:4120740/93‑1 (MQ=255)

                                                                                ggTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCaa         >  2:5084986/1‑93 (MQ=255)

                                                                                 gTGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCaaa        >  2:1857573/1‑93 (MQ=255)

                                                                                  tGCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTCGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAAt       <  2:3883753/93‑1 (MQ=255)

                                                                                   gCATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAAtg      >  1:5005788/1‑93 (MQ=255)

                                                                                    cATAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCGCATCTGTCAGGCAAAtgt     <  1:6364623/93‑1 (MQ=255)

                                                                                     aTAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAAtgtg    <  2:6036836/93‑1 (MQ=255)

                                                                                     aTAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAAtgtg    >  2:5216189/1‑93 (MQ=255)

                                                                                      tAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGc   >  1:4537137/1‑93 (MQ=255)

                                                                                      tAGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGc   <  2:659961/93‑1 (MQ=255)

                                                                                       aGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGcc  >  1:2090869/1‑93 (MQ=255)

                                                                                       aGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGcc  <  2:5475762/93‑1 (MQ=255)

                                                                                       aGCGCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGcc  <  2:3197576/93‑1 (MQ=255)

GGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCC  >  NC_000913/37242‑37421

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